72 FOLGE. 1958. BAND 1. HEFT 68 
hie 
OF MICHIGAN 
DER PHYSICS 
LIBR* 
| | | P H Y S | K | 
Gegründet 1799 durch F. A.C. Gren 2 
Fortgeführt durch L.W. Gilbert, J. C. Poggendorff, 
G. und E. Wiedemann, P. Drude, W. Wien, 
M. Planck, E. Griineisen, F. Möglich 
Der ganzen Reihe 456. Band 5 
Kuratorium: 
W.GENTNER, W.GERLACH, F. HUND, B. KOCKEL, 
M. v. LAUE, W.MEISSNER, W. PAUL, R. 'W. POHL, 
it ee R. ROMPE, R. SEELIGER, W. WEIZEL 
d 
t 
—H.KOPFERMANN 
| Heidelberg 
k 
N 
d 
3, 
r JOHANN 


= 
BR: 
avg 
r 
& 


Annalen der Physik. 7. Folge. Band 1. Heft 6-8. 1958 
Inhalt 
G. Oertel, Berlin: Messungen der Sekundärelektronenemission an Selen. Mit T. 
N.Mizkjewitsch, Moskau: Zu den Invarianzeigenschaften der Lagrange-Funk- 
H. Stolz, Berlin: Zum Phänomen der Entelektrisierung in kleinen leitenden Kör- 
pern in einem ausgedehnten elektromagnetischen Hochfrequenzfeld. Mit 
H. Grunewald und W. Neumann, Potsdam: Über die elektrische Leitfähigkeit 
von Blei(II)-oxyd mit Zusätzen höherwertiger gitterfremder Kationen. Mit 
J. Pachner, Prag: Die Bewegungsgleichungen der unitären Feldtheorie in der 
G. Brederlow, Greifswald: Der Potential- und Feldstärkeverlauf im Kathoden- - 
fallgebiet von Glimmentladungen. Mit 16 Abbildungen .................. 359 
S. Kastner, Berlin-Adlershof: Zur Theorie der Relaxation. I. Die thermodynami- u. 
schen Grundlagen und ihre mathematische Lösung. Mit 4 Abbildungen 377 De 
D. Lyons, Zeuthen-Miersdorf: Das Milnesche Problem bei anisotroper Streuung ar 
(Anwendung auf thermische Neutronen). Mit 3 Abbildungen .............. 400 hi 
M. Ziganescu, Timisoara (Rumänien): Der Hodograf der Bewegung eines ma- be 
teriellen Punktes in der relativistischen Mechanik. Mit 1 Abbildung ..... 424 al 
Br. Für diese Zeitschrift werden grundsätzlich nur Arbeiten angenommen, die vorher hi 
weder im Inland noch im Ausland veröffentlicht worden sind. Mit der Annahme des u 
i Manuskriptes und seiner Veröffentlichung geht das alleinige Recht der Vervielfältigung, m 
‘ rs Verbreitung und Übersetzung einschließlich des Rechts der Veranstaltung von Sonder- di 
ee drucken auf den Verlag über. Von jeder Arbeit liefert der Verlag 75 Sonderabziige 
er. Einsendung der Manuskripte erbeten an Prof. Dr. H. Kopfermann, Heidelberg, Albert- m 
Ueberle-Straße 7, Prof. Dr.G. Richter, Zeuthen-Miersdorf, Platanenallee 6. gr 
Die Redaktion bittet, die Manuskripte druckfertig einzuliefern und in den Korrek- se 
turen den bisherigen Raum des Textes nicht zu iiberschreiten. Die Zeichnungen sind g 
in sorgfältiger Ausführung auf besonderen Blättern beizulegen. 
Bei den Zitaten wolle man die Abkürzungen der Zeitschriftentitel nach dem Ver- 
zeichnis wählen, das jedem Jahrgang der ,,Physikalischen Berichte‘ vorangestellt ist. 
d 
Die Zeitschrift erscheint in Bänden mit 8 Heften, die zwanglos ausgegeben werden. a 
Die Lieferung erfolgt durch den Sortimentsbuchhandel, in Ausnahmefällen unmittel- 
bar vom Verlag. Der Bezugspreis beträgt je Band 34.— (U. S. $ 8.16) zuzügl. Postgebiihren. g 
Die Lieferung erfolgt bis zur Abbestellung, die nur fiir das Ende eines Bandes ausge- 7 
sprochen werden kann. 


f 


‘ 

> 
. ~ 
on 


334 


el- 


3 


ANNALEN DER PHYSIK _ 


7. FOLGE * BAND 1. HEFT 6-8 * 1958 


Inhaltsiibersicht 


Es werden mit einer statischen und einer dynamischen Methode die SEE 
von Selen, ihre Temperaturabhängigkeit und die Energieverteilung der SE 
mit verschiedenen Dotierungen von Brom und Tellur als Parameter unter- 
sucht. Dabei ergeben sich für die kristallinen Proben höhere Ausbeutewerte 
als für die amorphen. Die Beifügung von Brom und Tellur führt ebenfalls 
zu höheren 6-Werten, wobei die Maxima der Ausbeutekurven mit höherer 
Dotierung sich nach größeren Energien hin verschieben. Sowohl für die 
amorphe als auch für die kristalline Form ist die Ausbeute temperaturunab- 
hängig. Die Energieverteilungskurven haben ihre Maxima alle ungefähr 
bei 2,3 eV und fallen im wesentlichen nach höheren Energien zu um so steiler 


ab, je größer die Dotierung ist. 
a 

Es sollen der Koeffizient 6 der Sekundärelektronenemission (SEE) in Ab- 
hängigkeit von der Beschleunigungsspannung der Primärelektronen (PE) 
und seine Temperaturabhängigkeit an sehr reinem Selen und an Selenproben 
mit Brom- bzw. Tellurzusatz gemessen werden. An den gleichen Stoffen soll 
die Energieverteilung der Sekundärelektronen (SE) bestimmt werden. 

Die Abhängigkeit des SEE-Koeffizienten ö von der Energie der PE an 
Selen wurde bereits von Wolff!) und von Gobrecht und Speer?) unter- 
sucht. Eigene Messungen bestätigen und erweitern diese schon vorhandenen 
Ergebnisse. Bei den Untersuchungen wurde auf saubere präzis bekannte Aus- 
gangsmaterialien und, um eine Einlagerung von Sauerstoff zu vermeiden, auf 
sehr gutes Vakuum beim Aufdampfen der zu untersuchenden Schichten 
großer Wert gelegt. 

Er 2. Beschreibung der untersuchten Selenproben 


a Tab. 1 zeigt die Daten der spektralanalytischen Untersuchungen fiir 
die verwendeten Selenproben®). Die mit Se M 1 bezeichnete Reinstselenprobe 
ist sehr sauber, die dotierten Proben zeigen größere Verunreinigungen. 

Beim Aufdampfen des Selens auf den Probenträger, dessen Temperatur 
gleich der Zimmertemperatur ist, ergibt sich die amorphe Modifikation, die 


4 H. Wolff, Ann. Physik 89, 591 (1941). 


 *) H. Gobrecht u. F. Speer, Z. Physik 135, 602 (1953). 
s) Herrn Dr. F. Eckart danke ich für die freundliche Überlassung der Selenproben. 
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Tabelle 1 die \ 


Selenproben mit Brom- bzw. Tellurzusatz anzu 
ss Ergebnisse der spektralanalytischen Untersuchungen Mean kath 
SeM1| $eM2| SeM3 | SeM4 | SeM5 sank 
a Reinst- 0,098% | 0,78%Br |0,78% Br 0,78% Br} SeM6 | SeM7 geru 
=> SeBr,- | 0,04% | 0,12%, | 0,24% | 0,1% Te | 0,2% Te verg 
| Zusatz SeBr, SeBr, | SeBr, span 

Element | Kennzahl | | Kennzahl |Kennzahl |Kennzahl |Kennzahl | Kennzahl| Kennzahl und 


beim Tempern der Schicht in den kristallinen hexagonalen Zustand übergeht. 
Die Untersuchung des SEE-Koeffizienten an amorphen und kristallinen hexa- 
gonalen Selenproben wurde mit Hilfe einer statischen und einer dynamischen 
Methode durchgeführt. 


| | | | mom 
| | | 
| | | | | | 
| | | | mu 


3. Die statische Methode ' 


3.1. Die experimentelle Anordnung. Abb. 1 zeigt das verwendete Stu: 
Versuchsrohr. Es besteht im wesentlichen aus Strahlsystem, Probenträger Nei 
und Kollektor (innen graphitierte Kugel aus Kupfer-Nickel-Blech). Nach der Nac 
mechanischen Fertigstellung des Rohres wird es unter Hochvakuum (10~ Torr) gon 
8 Stunden bei 380° C ausgeheizt*). Das Selen wird aus einem Quarzschiffehen Ver: 
auf einen Probenträger aus Tantal aufgedampft. Vorversuche zeigten, daß cheı 


sich Tantal als Metall für den Probenträger deshalb eignet, weil die Wahr- best 
scheinlichkeit des Abplatzens der Selenschicht gering ist. In Stellung 1 Elel 
(Abb. 1) wird der Träger mit Selen bedampft, in Stellung 2 wird die Messung die 
durchgeführt. Da zunächst die amorphe Modifikation des Selens untersucht Gal 


werden soll, ist das Quarzschiffehen 80 mm vom Träger 7 entfernt, um zu Zus 
vermeiden, daß durch die Einwirkung der Wärmestrahlung das Selen kri- stel 
stallisiert. Es erwies sich als notwendig, während des Aufdampfens von Selen Ten 
pot 
% Bar A. Werner hat den mechanischen Aufbau der Röhren, Herr M. Müller . 
die glastechnischen Arbeiten ausgeführt. Beiden Herren danke ich auch an dieser Stelle 
recht herzlich dafür. erg: 
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die Wolframdrahtkathode des Strahlsystems zu heizen und Anodenspannung 
anzulegen, d.h. Strom zu ziehen, da sonst die Emission der Wolframdraht- 
kathode zu klein wurde und der Kathodenstrom unter 1% des Anfangswertes 
sank. Diese Erscheinung wird vermutlich durch absorbiertes Selen hervor- 
gerufen. Eine Oxydkathode wird durch das Selen auch dann vollkommen 
vergiftet, wenn während des Aufdampfens von Selen Heiz- und Anoden- 
spannung angelegt werden. Abschließend wird das Versuchsrohr gegettert 


d ezogen. 


RE 


— 


Kollektor 
wert 


» 


# Abb. 1. Versuchsrohr für die statische Methode 


Zum Kristallisieren wurden die Proben in den abgezogenen Röhren mehrere 
Stunden bei 75 bis 90° C getempert. Dabei zeigten dünne Schichten mehr 
Neigung abzuplatzen oder aufzureißen als dick aufgedampfte Selenschichten. 
Nach dem Heizen hatten die Schichten die metallisch graue Farbe des hexa- 
gonalen Selens. Zur Kontrolle der Kristallisationsprozesse wurde ein weiteres 
Versuchsrohr verwendet, dessen Träger aus einem quadratischen Glasplätt- 
chen (1cm?) bestand, das außer einem Streifen von 1mm Breite mit Platin 
bestäubt war. Auf dieses Plättehen wurde Selen aufgedampft und an die 
Elektroden eine Spannung von 100 V angelegt. Mit dieser Anordnung wurde 
die Abhängigkeit des Widerstandes von der Temperatur mit einem Multiflex- 
Galvanometer MG 3 (Empfindlichkeit: 2. 10-9 A/mm) gemessen. Im amorphen 
Zustand reicht die Meßempfindlichkeit nicht aus, um Änderungen festzu- 
stellen. Beim Übergang vom amorphen zum kristallinen Zustand bei einer 
Temperatur von 75 bis 90°C sank der Widerstand um mehrere Zehner- 
potenzen. 


‘ Ein Versuchsrohr mit einem Leuchtschirm an Stelle des Probenträgers 
ergab, daß der Elektronenstrahldurchmesser auf dem Probenträger ungefähr 
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1mm war. Das entspricht bei einem Strahlstrom in der Größenordnung von 
1 uA einer Belastung der Probe von ungefähr 0,8 wA/mm?. 

3.2 Meßverfahren. Es wird die übliche Versuchsschaltung verwendet, 
Der Kollektor K ist geerdet und liegt auf Anodenpotential. Die Anoden- 
spannung U, liefert ein 4-kV-Netzgerät, dem ein magnetischer Spannungs- 
konstanthalter vorgeschaltet ist. Die Wehnelt-Spannung ist zwischen 0 
und —30 V regelbar, der Heizstrom beträgt 2,5 A. Alle Zuführungen sind 
durch Potentialringe vor Fehlströmen geschützt. Um die elektronenoptischen 
Verhältnisse konstant zu halten, wird der Strahlstrom durch Variation der 
angelegten Heizspannung geändert. Die Probe P kann wahlweise + oder 
—50 V Vorspannung gegen den Kollektor erhalten. Die Ströme J, (Strom 
über die Probe) und /, (Strom über den Kollektor) wurden mit je einem Multi- 
flexgalvanometer MG 3 gemessen. 

3.3. Ergebnisse. 3.3.1 Definitionen. 

In Übereinstimmung mit sehr vielen Arbeiten über SEE setzen wir die 
Grenze zwischen wirklichen SE und rückdiffundierten Elektronen auf 50 eV 
fest. 

Auftretende Tertiärelektronen, das sind am Kollektor durch SE ausge- 
löste Sekundärelektronen, werden häufig durch ein vorgespanntes Netz 
aufgefangen. Da bei unserer Anordnung ohne Netz das Kugelinnere mit 
einem Aquadagbelag versehen ist (Kohlenstoff weist eine niedrige Ausbeute- 
kurve auf, dmax < 1), können diese bei der Auswertung vernachlässigt werden. 

Die Frage, ob die Anzahl der den Kollektor erreichenden rückdiffun- 
dierten Elektronen mit einer Energie größer als 50 eV in beiden Fällen, d.h. 
wenn die Probe einmal positiv und einmal negativ gegen Anode vorgespannt 
ist, gleich ist, wurde mit Hilfe eines Modells der Anordnung im elektrolytischen 
Trog geklärt. Die Konstruktion der Elektronenbahnen aus den gemessenen 
Potentiallinien ergab, daß bei Gegenspannung alle Elektronen mit einer Energie 
von (50 + 0,1) eV, die mit einem Winkel kleiner als 77° gegen das Lot aus- 
treten, den Kollektor erreichen. Berücksichtigt man dabei noch, daß die 
Anzahl der in einem gegebenen Raumwinkelelement austretenden SE mit 
cos g abnimmt, wenn g der Winkel gegen das Lot ist, so kann obige Frage, 
ohne eine große Ungenauigkeit zu begehen, bejaht werden. 

Für die gemessenen Ströme J, und J, gilt folgender Zusammenhang. Bei 


positiver Vorspannung ist N 
=i,-i, It =i, 4 
und bei negativer Vorspannung von P u eu 


Dabei sind i, der Anteil der reflektierten oder nur einmal unelastisch gestreuten 


PE mit einer Energie größer als 50 eV, isg der Strom, der von den eigent- 
lichen SE herrührt, i, der Strahlstrom der PE.” Daraus ergibt sich für die 


Berechnung von 6: 
« ty — ty I5 


Nach dieser Formel wird im folgenden für die statische Methode die Ausbeute 
berechnet. 

3.3.2 Die Ausbeutekurven der untersuchten Proben. Abb. 2 zeigt die 
Ausbeutekurve von amorphem Reinstselen (Probe SeM 1). Diese Probe 
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von wurde mehrmals reproduzierbar, erst mit steigender, dann mit fallender Primär- 

spannung gemessen. Dabei liegen die Werte für abnehmende Anodenspannung 
let. wenig aber systematisch höher als für zunehmende. Bei nochmaliger Messung 
len- mit steigender Anodenspannung liegt die Ausbeutekurve nahezu mit der- 
1gs- jenigen fir abnehmende Beschleunigungsspannung zusammen. 
n 0 
ind d 
hen 
der 2 
der 
om 
ılti- 
eV 
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ge- 
letz E : 8 ==} omorph 
en. 05+ 
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ais Abb. 2. Mit steigender und fallender Primärspannung gemessene Ausbeutekurven 
nen 
rgie In Abb. 2 sind weiter die Ausbeutekurven des kristallinen Selens (Se M 1), 
Lus- das durch Aufdampfen bei 80 + 5°C gleich den kristallinen Zustand (kri- 
die stallin I) annimmt und der kristallinen Modifikation, bei der zunächst bei 
mit Zimmertemperatur aufgedampft und dann durch Tempern die amorphe in die 
ge, kristalline Form umgewandelt wird (kristallin IT), aufgetragen. Die Maximal- 


werte der kristallinen Form liegen um 0,06 bzw. 0,10 für ö höher als die- 
Bei jenigen der amorphen Modifikation. Der zweite (6 = 1)-Punkt liegt bei kri- 
stallinem Selen bei 1900 V, beim amorphen bei 1600 V. Der Vergleich der 
beiden verschieden gewonnenen kristallinen Zustände zeigt, daß im Maximum 6 
bei Selen (kristallin I) bei 1,46 um 0,05 höher liegt als bei Se (kristallin II). 

In Abb. 3 sind die ö-Werte für die amorphen Modifikationen von Se M 1 
(noch einmal zum Vergleich), Se M 2, Se M 3, SeM 4, Se M 5 aufgetragen. Die 
ten ö-Werte wachsen mit zunehmender Fremdkomponentenbeimischung, und die 
nt- Maxima zeigen eine leichte Verschiebung zu höheren Primärenergien. Abb. 4 
die zeigt die gleiche Tendenz für die Tellurbeimischungen SeM 6, SeM 7 wie 
Abb. 6 für die Bromzusätze. Für die Abweichungen der Meßpunkte von der 
gestrichelten Kurve in Abb. 4 haben wir keine Erklärung. Abb. 5 gibt den 
Verlauf der 6-Werte für die kristalline (IT) Form für die Proben SeM 1, 4, 5, 6, 7 


ute an. Aus experimentellen Gründen ist die Form ,,kristallin II“ gewählt worden, 

da auf diese Art und Weise mit einem Versuchsrohr die amorphe und die 
die kristalline Modifikation gemessen werden konnten. Die Ausbeutekurven fiir 
obe die kristallinen Formen zeigen ähnliche Tendenz wie diejenigen fiir die amorphen 
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Die Verschiebung des Maximums mit zunehmender Dotierung nach rechts 
ist noch etwas stärker ausgeprägt als bei den amorphen Modifikationen. Die 
Maxima liegen 0,04 bis 0,06 für 6 höher als die entsprechenden Werte bei den 
amorphen Körpern. 
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Abb. 4. Ausbeutekurven der mit Te dotierten amorphen Se-Proben 


3.3.3 Die Temperaturabhängigkeit der Ausbeute. Zur Bestimmung der 
Temperaturabhängigkeit der Ausbeute wurde das Versuchsrohr in einem 
doppelwandigen Aluminiumblechbehälter aufgestellt, dessen Zwischenraum 
mit festem Kohlendioxyd angefüllt wurde. Zur Wärmeisolierung wurde der 
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Behälter mit einer 10 cm dicken Schicht aus Zellstoff umgeben. Für die Mes- 
sungen bei höherer als Zimmertemperatur wurde das Versuchsrohr im Aus- 
heizofen untergebracht. Schwierigkeiten ergaben sich hierbei für die Mes- 
sungen bei abnehmender Temperatur, nachdem das Selen in die kristalline 
Form übergegangen war, da sehr häufig auf Grund der Trägheit des verwen- 


100 500 7000 1500 2000 wiv) 
Abb. 5. Ausbeutekurven der kristallinen Se-Proben 


leten Ofens die Schicht nach dem Kristallisieren sublimierte und deshalb die 
Messungen mit fallender Temperatur nicht durchgeführt werden konnten. 
Aus den Messungen ergab sich, daß die SEE des Selens von der Temperatur 
unabhängig ist, gemessen von —78° C bis zu derjenigen Temperatur, bei der der 
Kristallisationsprozeß einsetzt. Dieser Temperaturwert hängt von den Zu- 
sätzen zum reinen Selen ab. h 
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Abb. 6. Temperaturabhängigkeit der Ausbeute 


Es wurden die Proben Se M 1, 2, 4 untersucht, wobei es bei der Probe 2 
leider nicht möglich war, die Sublimation zu verhindern, so daß die Werte 
fiir abnehmende Temperatur nach der Kristallisation fehlen. Die Abb. 6 
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zeigt die gemessenen Temperaturabhängigkeiten der Ausbeute für die jeweils 
zum Maximum der Ausbeutekurve gehörige Anodenspannung. Danach ist 
also keine Temperaturabhängigkeit vorhanden. 

3.3.4 Die Energieverteilung der Ausbeute. Mehrere amorphe und kri- 
stalline Selenproben wurden mit einer Primärelektronenenergie von 100 eV 
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Abb. 7. Gegenspannungskurve von Se M 1 amorph bei U, — 100 V 
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Abb. 8. Spektrale Verteilung der SE bei der nicht dotierten Probe 
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untersucht. Es ist notwendig, in enger Abständen zu messen, da sich die 
Verteilungsfunktion 


% 


: aus der Differentiation der Gegenspannungskurve ergibt. Deshalb wurde die 
_ Gegenspannung U, bei der Messung jeweils um 0,5 V geändert und hierzu 

die am Kollektor noch ankommenden SE (i,) gemessen. 
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Die Abb. 7 zeigt eine Gegenspannungskurve (von amorphem SeM 1). 
In Abb. 8 ist die sich daraus ergebende spektrale Verteilung aufgetragen, 
wobei die Verteilungsfunktion im Maximum auf 1 normiert ist. Außerdem 
zeigt Abb. 8 die Verteilungsfunktionen für die beiden verschieden gewonnenen 
kristallinen Formen von Se M 1. Diese beiden Kurven sind nahezu identisch. 
Nach höheren Energien fällt die Funktion für die kristalline Form jedoch 
steiler ab als diejenige für die amorphe. Abb. 9 enthält die Geschwindig- 
keitsverteilung der SE für die amorphen Formen von SeM 1, 3, 4,5. Mit 
Ausnahme des Verhaltens der Probe Se M 5 scheinen die Energieverteilungs- 
kurven um so steiler abzufallen, je höher die Dotierung ist. 


fleu) & 


Abb. 9. Spektrale Verteilung der SE bei den dotierten Proben 


4. Die dynamische Methode 


4.1 Das Versuchsrohr. Das für die Untersuchungen mit der dynami- 
schen Methode verwendete Versuchsrohr enthält ein Strahlerzeugungs- 
system mit magnetischer Fokussierung und Ablenkung und die Signalplatte 
(Probenträger) mit einem Absaugring für die entstehenden SE. Die Innenwand 
des Versuchsrohres ist metallisiert und liegt auf Anodenpotential, die Strahl- 
kanone hat eine thorierte 
Wolframdrahtkathode. 

4.2 Der Proben- 
träger. Auf einem Glas- 
objektträger wird bei einer 
Temperatur von 300°C 
eine salzsaure Lösung von 
Zinntetrachlorid aufge- 
stäubt. Durch das ent- 
stehende Zinndioxyd er- 
hält das Glas eine sehr 
dünne leitfähige Ober- 
flächenschicht, die sehr Abb. 10. Probenträger mit Maske => 
widerstandsfähig ist. 

Eine Hälfte des Objektträgers wird im Rohr durch eine zwise chen zwei 
Glimmerscheiben liegende Konstantandrahtwendel geheizt. Mit einem am 
Objektträger befestigten Kupfer-Konstantan-Thermoelement kann die Tem- 
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- peratur bestimmt werden. Vor dem Einbringen des Objekttragers in das 
Versuchsrohr werden zu Fichzwecken noch je zwei Quadrate von 5 x 5mm? 
mit Platin bzw. Titan bedampft. Das Versuchsrohr wird bei einer Tem- 
peratur von 380°C 8 Stunden lang evakuiert. Mit einem magnetisch einzu- 
führenden Schiebeofen wird aus einem Quarzschiffehen der Probenträger mit 
Selen bedampft, wobei die mit Platin bzw. Titan versehenen Quadrate 
durch eine Maske abgedeckt werden (Abb. 10). Die untersuchten Schichten 
hatten eine interferometrisch bestimmte Dicke von 400-600 A. 

4.3. Durchführung der Messung. Abb. 11 zeigt die Meßanordnung 
zur fernsehmäßigen Abtastung im Blockschaltbild. Die Selenschicht wird 
_ auf ihrer Vorderseite mit einem Elektronenstrahl abgetastet, der wie in einer 
Fernsehröhre zeilenweise über die Schicht bewegt wird. Das geschriebene 
Raster kann in seiner Größe geändert werden, so daß nicht nur die gesamte 
Signalplatte, sondern auch einzelne Gebiete auf ihr mit entsprechender Ver- 

_ größerung (Dehnung in vertikaler und horizontaler Richtung) elektronisch 
sichtbar gemacht werden können. Der Strahlstrom wird während der Rück- 
laufzeit ausgetastet. Der unmittelbar vor der Signalplatte befindliche Draht- 
ring (Absaugring) ermöglicht eine bessere Absaugmöglichkeit für die von der 
Selenschicht stammenden SE und verhindert damit das Auftreten von Stör- 

signalen. Die optimalen Werte der Vorspannungen von Signalplatte bzw. 
 Absaugring liegen bei —40 bis —50 bzw. bei +15 bis +25 V gegen Anode. 
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Die bei der Abtastung der Signalplatte entsprechend den verschiedenen 
Ausbeutekoeffizienten entstehenden Stromimpulse rufen an dem Arbeits- 
widerstand AR, einen Spannungsabfall hervor, der nach genügender Verstärkung 
in der beim Fernsehen üblichen Weise die Helligkeit der Bildwidergaberöhre 
steuert. Außerdem wird in einem Meßoszillographen auf zwei getrennten 
Oszillographenréhren der Verlauf der Bildmodulationsspannung einer be- 
— waagerechten Zeile sowie innerhalb eines ; schmalen senkrechten 


Abb. 11. Meßanordnung zur fernsehmäßigen Abtastung 
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Streifens sichtbar gemacht. Zur Kennzeichnung der ausgewählten Zeilen 
und Spalten werden dem Kontrollbild zusätzliche Aufhellimpulse hinzugeführt, 
die im Bild einen entsprechenden waagerechten und senkrechten Streifen 
hervorrufen. 

4.4. Interpretation der Meßgrößen. Bei konstanter Verstärkung 
lassen sich mit Hilfe des Meßoszillographen die Amplituden A von ,,Schwarz- 
Weiß-Sprüngen‘“ im Verlauf der Bildsignalspannung messen. Die Eichung 
der Amplituden erfolgt (wobei genügend Linearität des Verstärkers und in der 
Ablenkung des Oszillographen vorausgesetzt wird) direkt in Maßzahlen für 
die 6-Werte. Da näherungsweise die Amplitudenwerte A als lineare Funktion 
der ö-Werte angenommen werden, also auch 6 = L (A) gilt, kann aus zwei 
bekannten Punkten für Platin und Titan: 


P, (Opt, Apt)v; Po (Ori, Arı)v 


für jede beliebige feste Primärspannung V aus den gemessenen Werten Age 
das zugehörige dg. bestimmt werden. 

4.5 Ergebnisse. Mit der dynamischen Meßmethode wurden folgende 
Ergebnisse gefunden. 

4.5.1 Die Ausbeute in Abhängigkeit von der Energie der PE. Es wurden 
die Proben SeM 1 und SeM 5 untersucht. Mit Hilfe der bekannten Ausbeute- 
kurven von Platin und Titan wurden die 6-Werte der Proben zwischen 800 
und 2350 (1750) eV ermittelt. Abb. 12 zeigt den Verlauf der Ausbeutekurven 
für die amorphen und kristallinen Zustandsformen der untersuchten Proben. 
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4.5.2 Die Temperaturabhängigkeit der Ausbeute. Bei der fernsehmäßigen 
Abtastung bleiben unterhalb der Kristallisationstemperatur alle Helligkeits- 
werte konstant. Dieses kann an den Bildsignaloszillogrammen von ausge- 
wählten Zeilen und Spalten gemessen werden. Beim Erreichen der Kristalli- 
sationstemperatur beginnt das Selen von einem Keimpunkt aus, der in der 
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Regel ungefähr in der Mitte der geheizten Fläche liegt, sich langsam in die 
‚kristalline Modifikation umzuwandeln, zunächst über die geheizte Fläche und 
dann auch über den nicht geheizten Teil des Trägers. Die Abb. 13 zeigt eine 
Phase eines solchen Kristalli- 
sationsprozesses. Gleichzeitig lassen 
ch an den Meßoszillographen die 
ıgehörigen Signaländerungen mes- 
on. Mit dieser Methode wurde 
uch festgestellt, daß die SEE des 
morphen Selens von der Tempe- 
itur im Bereich der noch zu dis- 
kutierenden Meßfehlergrenzen zwi- 
chen der Zimmer- und Kristalli- 
itionstemperatur unabhängig ist 
es wurden die Proben SeM 1, 5, 7 
untersucht). Da es nicht immer 
gelang, die Probenheizung so ab- 
Abb. 13. Eine Phase des Kristallisationspro- Zuschalten, daß zwar das Selen 
zesses. 1: Se rg 2: Se kristallin, 3: Titan, vollständig kristallisiert, aber noch 
4: Platin, 5: Aufhellimpulse nicht sublimiert war, wurde das 
; Konstantbleiben von 6 auf dem 
Wert für die kristalline Form bei abnehmender Temperatur nur bei Se M 1 


5. Diskussion 


- 
Die mit der statischen und mit der dynamise ee Methode gewonnenen 
__ Ergebnisse zeigen eine gute Ubereinstimmung (Abb. 2 und 12). Zum Vergleich 

unser er Resultate mit den Ergebnissen der Arbeit von Gobrecht und Speer?) 

schließen wir aus der Bemerkung, daß nach optischen Beobachtungen die 
_ Absorbierte Strahlungswärme zur Umwandlung der Schichten in den metalli- 
ne schen Zustand ausgereicht habe, daß die Verfasser die kristalline hexagonale 
_ Modifiktion des Selens untersucht haben. Weiterhin nehmen wir an, daß dem 
= vo wendeten Selen unsere undotierte Probe Se M 1 näherungsweise entspricht. 

j Wie eine Betrachtung der Ausbeutekurve des Selens nach Gobrecht und 

u eo und unserer 6-Kurve des kristallinen (II) Selens (Abb. 5) aufzeigt, 

stimmen die Werte recht gut überein. Das Maximum liegt bei der Messung 

von Gobrecht und Speer bei 400 (350 eig. Mess.) V mit dmax = 1,38 

a4 eig. Mess.), der zweite (6 = 1)-Punkt bei 1900 (1800 eig. Mess.) V. 


Hachenberg findet bei seinen Messungen (Abb. 14) in Kurve a ähn- 
die he Werte für die kristalline Modifikation des Selens®). Die geringen Ab- 
weichungen von 6 zu höheren Werten sind sicher dadurch zu erklären, daß 
das verwendete Selen nicht gleiche Reinheit hatte. Das Ergebnis einer spek- 
 tralanalytischen Untersuchung liegt leider nicht vor. Analog zu unseren Mes- 
- sungen am amorphen Selen Se M1 (die ö-Werte für abnehmende Primär- 
spannung liegen wenig aber systematisch höher als für die vorangehende 


7 5) Herın Prof. Dr. Hachen berg danke ich für diese private Mitteilung aus einer un- 
veröffentlichten Arbeit. 
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Messung mit steigender Anodenspannung) findet auch O. Hachenberg 

ein Anwachsen der 6-Werte nach Elektronenbeschuß (Abb. 20b und 20c). 
Aus dem Vergleich der Abb. 5, 6, 7, 8 geht hervor, daß einbeitlich die Aus- 

beutewerte für die kristallinen Formen höher liegen als die zugehörigen Werte 

für die amorphen Formen. Möglicherweise ist die größere Beweglichkeit der 

Elektronen im kristallinen Zustand für 

das höhere 6 verantwortlich. 

Mit zunehmender Dotierung — so- t 
wohl von Brom als auch von Tellur — 4 
steigt das Maximum der Ausbeutekurve en. 
an (Abb. 3, 4,5), außerdem verschiebt y : 
sich das Maximum für höhere Dotie- ~~ ¢ 
rungen zu größeren Primärenergien hin. 7 a == 
Der Effekt, daß mit wachsender Dotie- 
rung die ö-Werte von Selen zunehmen, 
wird auch von Gobrecht und Speer?) 
für Beimischungen von Jod und Brom 
(je 0,03 %) festgestellt. Die Interpretation Sth. 14. Messung von 0. Hachenberg 

an Se. a und b: Se auf warme Träger 
von Gobrecht und Speer, nach der „ufgedampft, c: b nach Elektronen- 
die Elektronennachlieferung ins Leitungs- beschuß 
band durch eine höhere Störstellen- 
konzentration begünstigt wird und dieser Einfluß die ausbeuteverringernde 
Wirkung der gleichzeitig wachsenden Ladungsträgerkonzentration überwiegt, 
beschreibt das Verhalten richtig. 

Ein Vergleich der auf verschiedene Weise gewonnenen kristallinen For- 
men I und II (Abb. 2) zeigt, daß die Form I, die durch Aufdampfen des Selens 
auf warme Träger entsteht, größere Ausbeutewerte aufweist. Eine Erklärung 
für dieses Verhalten ist vermutlich darin zu suchen, daß das kristalline Selen 
eine größere Dichte (y = 4,82 gem?) hat als das amorphe (y = 4,26 bis 
4,32 gem?) und beim Tempern der Schicht die Oberfläche schrumpft und 
uneben wird. Nicht glatte Flächen ergeben aber bekanntlich®)’) kleinere 
ö-Werte als glatte Flächen aus gleichem Material. 

Sowohl die Ergebnisse aus der statischen als auch aus der dynamischen 
Methode zeigen, daß für die amorphen und kristallinen Formen der Selen- 
proben die Ausbeute von der Temperatur unabhängig ist. Dieses Verhalten 
stimmt überein mit einem Resultat, das Hachenberg?) in einer theore- 
tischen Untersuchung, in der er die Relaxationsstrecken für die Wechsel- 
wirkungsprozesse bei der SEE abschätzt, erhält. Dabei ergibt sich, daß 
für Körper mit hoher Fehlordnung, dazu gehört zweifellos Selen, die Relaxa- 
tionsstrecke nicht mehr durch die Wechselwirkungen mit den Gitterschwin- 
gungen beeinflußt wird. D.h. für solche Stoffe muß die SEE von der Tem- 
peratur unabhängig sein. 

Aus Abb. 9 geht hervor, daß die Umwandlungstemperatur, bei der der 
amorphe Zustand in den kristallinen übergeht, mit wachsender Dotierung 
steigt. Die zur Umwandlung notwendige Zeit nimmt mit wachsender Do- 


6) H. Bruining, Die Sekundärelektronenemission fester Körper, Springer, Berlin 
1942. 

7) O. Hachenberg u. W. Brauer, Fortschritte der Physik 1, 439 (1954). 

8) O. Hachenberg, Ann. Physik 2, 403 (1948). 
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tierung ab. Eine gleiche katalytische Beeinflussung des Kristallisations- 
prozesses finden de Boer®) bei Zusatz von Br und Te zu Selen, und Eckart 
und Rabenhorst?) bei Zusatz von Br. 

Abb. 9 zeigt, daß die Energieverteilungskurven nach höheren Energien 
zu um so steiler abfallen, je höher die Dotierung ist. Die Annahme, daß die 
Austrittstiefe x, der SE mit der Dotierung wächst, würde dieses Verhalten 
dahingehend klären, daß die Energieverteilung für größere x, stärker homo- 
genisiert ist. Gleichzeitig würde diese Annahme zusammen mit der Be- 
merkung von Decker"), daß für große Primärenergien 6 proportional x, ist, 
das Anwachsen und Verschieben zu höheren Primärenergien der Maxima der 
Ausbeutekurven mit zunehmender Fremdkomponente erklären. Bemerkenswert 
ist weiterhin, daß die Maxima der Energieverteilungskurven alle bei ungefähr 
2,3 eV liegen, wobei allerdings noch zu berücksichtigen ist, daß die Elektronen- 
austrittsarbeiten der untersuchten Proben voneinander verschieden sein 
können. 


Herrn Prof. Dr. Hachenberg danke ich herzlich für die Anregung zu der 
vorliegenden Arbeit. Herrn Dr. Eckart danke ich für zahlreiche fördernde 
Diskussionen. 


9) F. de Boer, Rec. Trav. chim. Pays-Bas 62, 151 (1943). 

10) F. Eckart u. H. Rabenhorst, Ann. Physik 19, 381 (1957), : 

4) A. J. Decker, Physic. Rev. 94, 1179 (1954). = 


Berlin, Deutsche Akademie der Wissenschaften, Heinrich-Hertz-Institut. 


Oo Bei der Redaktion eingegangen am 23. September 1957. 
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Zu den Invarianzeigenschaften der Lagrange-Funktione 


der Felder 
Nikolaj Mizkjewitsch 
Inhaltsiibersicht 


Die vorliegende Arbeit konzentriert sich im wesentlichen auf die folgenden 
drei Probleme: 

1. Untersuchung der analytischen Form der Invarianzbedingung für 
Funktionen von Feldgrößen und ihren ersten beiden Ableitungen. Diese 
Funktionen werden mit den physikalischen Lagrange-Funktionen iden- 
tifiziert. 

2. Klassifizierung der erhalten bleibenden Größen mit Hilfe einer neuen 
Formulierung des Noether-Theorems. 

3. Bildung eines neuen Typs der einheitlichen Feldtheorie für Felder mit 
verschwindender Ruhmasse der Quanten. 


$ 1. Einführung 


Von der Invarianz der Wirkungsintegrale ausgehend kann man laut 
Noetherschem Theorem!) die ganze Information über Erhaltungsgesetze 
der physikalischen Größen folgern. Es ist leicht eine äquivalente Formulierung 
dieses Theorems zu finden, in der nur die Invarianzeigenschaften der La- 
grange-Funktion benutzt sind. Das läßt unmittelbar die differentialen 
(in irgendeinem Punkt des Raums gegebenen) Erhaltungssätze ohne integrale 
Ausdrücke folgern. Insbesondere erhält man automatisch eine allgemein- 
relativistische Verallgemeinerung des Drehimpulses. 

In der allgemeinen Relativitätstheorie sind die Erhaltungssätze von 
Einstein, Hilbert, Klein, Levi-Civitä, Lorentz und andere aus- 
fükrlich untersucht worden. Es ist die Diskussion zwischen Einstein einer- 
seits und Lorentz und Levi-Civitä andererseits über die Erhaltung von 
Energie und Impuls mit Rücksicht auf die Gravitation bekannt?). Es scheint, 
daß dieses Problem nur formal gelöst ist und es möglich ist, ein neues Licht 
in diese Frage zu bringen. Man kann nämlich einige Resultate der speziellen 
Relativitätstheorie auf dieses Gebiet verallgemeinern. 

Außer der natürlichen Klassifizierung der erhaltenen Größen läßt die 
neue Formulierung des Noetherschen Theorems ein tieferes Studium des 
Problems der inneren Einheit der physikalischen Felder zu. Man kann ein- 
deutig die konkreten Formen der Lagrange-Funktionen für eine Klasse 


1) E.Noether, Götting. Nachr. 235 (1918). 
2) W. Pauli, Relativitätstheorie, Moskau 1947. 
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von Feldern angeben, wenn zur Forderung der Invarianz der Lagrange- 
Funktion (in allgemeinem Sinn) noch ein natürliches Einfachheitsprinzip 
für die Lagrange-Funktionen hinzugefügt wird. 


§ 2. Allgemeine Beziehungen 


Wir betrachten allgemein-relativistische Theorien der Felder, die durch 
die Potentiale A, (Wellenfunktionen) charakterisiert sind. Der Index i durch- 
läuft die Gesamtheit aller Wellenfunktionen. Diese können auch nicht-ten- 
sorielle Eigenschaften haben. Wir brauchen nur die Kenntnis ihrer Trans- 
formationsgesetze bei kontinuierlichen Transformationen der Koordinaten, 
d.h. die Änderungen 

0A; = Aj (a) — A; (2), 3 
die die Transformationen bis 


erzeugen, sind gegeben. (Man muß betonen, daß unsere 6A; eine andere 
Natur haben als diejenigen beim Wirkungsprinzip. Zum Beispiel verschwin- 
den sie an der Systemgrenze nicht.) 

In der Definition (1) haben wir berücksichtigt, wie die Komponenten 4, 
in verschiedenen Koordinatensystemen miteinander verbunden sind. Zum 
Beispiel, wenn A, ein Tensor ist, haben wir 

00x” 


ox" 02” Ox 


Die Operation 6 hat einen Mangel, weil sie mit dem Vektor Nabla nicht kommu- 


tiert. Deshalb fiihren wir die Operation 6*: En lf 
= 6A, — Ajo Ox" 

ein, die mit der partiellen Differentiation vertauschbar ist. Pr em 

Wir bezeichnen die Ableitung @4;/@a* mit 4; .. 

Im folgenden betrachten wir nur solche A,, die wie ate ig 

= 


transformiert werden. Diese Eigenschaft haben alle Tensoren (und Pseudo- 
tensoren), Tensordichten und im gewissen Sinn auch Spinoren. 

Jetzt suchen wir Bedingungen der Invarianz einer Funktion L dieser 
Potentiale und ihrer ersten beiden Ableitungen: 


L=L (A;; A;, xß> Juvs Juva: Jur, (6) 


Wir schreiben hier den metrischen Tensor g,, gesondert, weil man ohne ihn 
keine nicht triviale Invariante konstruieren kann. Die Form der Funktion L 
möge in allen Systemen der Koordinaten unveränderlich sein. 

Die Invarianzbedingung für Z (Definition eines Skalars) 
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der aus der Gl. (3) folgt, zur zen 


= 0 (9) 
umgebildet. Wir nehmen hier die lide 
(10) 
an. 


Infinitesimal kann man schreiben 


5g 

6A, A, Juv (sz Fr ax? ö*A, (11) 

Wir die bekannte | 


(Lagrange-Euler-Ableitung). 
von (5) und (11) die Form 


202° 
+5 lye oar + 


Wir führen hier die folgenden A ein: 


re 


v 


Die Gl. (9) bekommt dann bei Benutzung 


a bx? 


— th 
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Weil die Koordinatenvariationen und deren Ableitungen willkürlich sind, 
bekommen wir aus der Gl. (13) die folgenden Gleichungen: 


1 au? 
aire (21) 


Die geschweiften Klammern, die Indizes umfassen, bezeichnen die Summe über 

alle Permutationen, wodurch Symmetrie in allen Indizes erzeugt wird. 

Die Gl. (18) ist stets eine Identität, was man bei Verw er von eee und 

at 

bi -=—A;,+ Aj, ap - Aj apy 


(22) 
du Peli er eg 
Die (19—21) führen zu 
a . 
gout: Jar eu, 0 iy? te 4 (23) 
so daB die Identitat (18) die Form 7 
— gun + Te) = 0 (24) 


bekommt. Man kann einfach zeigen, daß diese Relation mit den Identitäten 
von Bianchi verbunden ist. 

Man kann die Funktion Z ,,Lagrange-Funktion“ nennen, wenn die mit L 
gebildeten Wirkungsintegrale bei Variation die Feldgleichungen geben, die 
mit den experimentellen Fakten übereinstimmen müssen. 

Wir zeigen jetzt, daß die Invarianzbedingung (7) vollständig äquivalent 
der Invarianz der Wirkungsintegrale 

8% =J - (de’) =J Q (x) (dx) = 0 (25) 


ist, wenn das Integrationsgebiet R willkürlich ist. Weil die Funktional- 


determinante der Transformation (2) mit der Genauigkeit erster Ordnung 


taver é (x) 


ist, bekommen wir die Gleichung 
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Wegen der Willkür von R muß daraus die GI. (9) folgen. Die Gin. (7) und (25) 
sind also äquivalent, wenn das Gebiet willkürlich ist. Also kann man alle 
Folgerungen des Noether-Theorems ohne Integralausdrücke bekommen. 
Die Feldgleichungen folgen dabei aus dem Wirkungsprinzip; aber das kann 
nichts verändern. 

Fordern wir, daß außer der Invarianz die Lagrange-Funktion zu par- 
tiellen Differentialgleichungen von nicht höher als der zweiten Ordnung 
führen sollen. Man hält es gewöhnlich für offensichtlich, daß diese Bedingung 
der Forderung, daß LZ nur von den Potentialen und deren ersten Ableitungen 
abhängig sein möge, äquivalent ist. Es ist aber klar, daß Gleichungen der 
zweiten Ordnung auch aus dem J, das linear von den zweiten Ableitungen 
abhängt, folgen, wenn die Koeffizienten bei den zweiten Ableitungen der 
Potentiale nur von den Potentialen und nicht von ihren Ableitungen ab- 
hängen. Dieser Fall ist nämlich realisiert bei der Lagrange-Funktion des 
metrischen Feldes. 


Wir definieren die Lagrange-Funktion für das freie Gravitationsfeld als 
L, = L (0; 0; 0; 9,,; Quv,a; Qur,apß)- (28) 


Fügen wir dazu L, als Lagrange-Funktion des übrigen Feldes, so erhalten 
wir die Lagrange-Funktion des 


Alle drei Größen L,, Z, und L, sind invariant bei den willkürlichen Trans- 
formationen der Koordinaten, so daß alle von uns erhaltenen Gleichungen 
für jedes Z gültig sind. So zerfallen die Größen U3, MP uv. a. jeweils in drei 
Terme: 1,3, U,3, Uys; Mid, mee, mz? u.a. Im weiteren lassen wir die 
Indizes t, g oder f in Gleichungen weg, wenn keine Mißverständnisse auf- 
treten können. 

§ 3. Erhaltungssätze 


In der speziellen Relativitatstheorie werden die differentiellen Erhaltungs- 
sätze immer so geschrieben, daß die vierdimensionale Divergenz der erhaltenen 
Größe verschwindet. Solche Größen sind in der Regel Tensoren. Wenn man 
aber die Verallgemeinerung für den Riemannschen Raum durchführt, 
scheint es, daß der Übergang zur kovarianten Ableitung unbefriedigend ist, 
weil die Gleichungen dann keine unmittelbare Erhaltung bedeuten, wie es 
bei Benutzung des Gauß-Theorems offensichtlich wird. Darum reduzierten 
Einstein und nachfolgende Autoren die Erhaltungssätze wegen der Bilanz 
der Gravitationsenergie auf die Form der gewöhnlichen Divergenz. 

Zugleich braucht die Größe, die auf diese Weise erhalten wird, keine ten- 
sorielle Eigenschaften bei willkürlichen Transformationen der Koordinaten 
haben, weil das Äquivalenzprinzip gültig ist. Lorentz und Levi-Civitä 
versuchten in ihrer Diskussion gegen Einstein den tensoriellen Charakter 
der erhaltenen Größen (Energie, Impuls und Spannung) beizubehalten. In 
diesem Punkt dürfte ihr Fehler liegen. 

Unsere Größen UF, 2° und RF" spielen eine grundlegende Rolle in der 
Feststellung der Erhaltungssätze und besitzen die tensoriellen Eigenschaften 
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nur bei den linearen Koordinatentransformationen (affine Tensoren). Im 
folgenden gebrauchen wir auch noch die nichtkovariante Größe — ,,Radius- schr 
vektor“ x*, deren partielle Ableitungen den unitären Kronecker-Tensor 
bilden. Mit Hilfe dieses „Radius-Vektors‘ kann man die folgenden Divergenz- 
gleichungen als Folgerung der Beziehungen (19—21) erhalten: 


han — fer (31) d 
ax? » ave che er 
Die bedeuten: I 
vi 
Mit Ausnahme der Erhaltungssätze (23), (30), (31) (,,starke Erhaltung‘), syst 


dessen physikalischer Sinn später erklärt wird, ist die Erhaltung der physika- 
lischen Größen nicht nur mit der Invarianz der Lagrange-Funktion, sondern puls 
auch mit der Erfüllung der Euler-Lagrange-Feldgleichungen verbunden 
(„schwache Erhaltung‘). 


auch 

Für das Gravitationsfeld g,, und für alle anderen Felder A, können wir 
entsprechend folgende Gleichungen schreiben: 
tens 


3 
54, 54, u 38 5) Dan 
Aus (14) und (34) folgt die Gleichung 


sofort. Der 
A. Erhaltene Größen von der Art der Energie 
Der Erhaltungssatz fiir Energie und Impuls entspricht der Invarianz der = 
Lagrange-Funktion bei Koordinatentranslationen 
6x* = const. (37) 
so d 
Die Gl. (13) ist dann auf (18) zurückgeführt. Alle anderen Gleichungen sat: 
bleiben wohl in Kraft, weil die Lagrange-Funktion ein allgemeiner Skalar Gel 
ist. Im entgegengesetzten Falle kann die Gl. (23) für alle drei Z bestätigt 
sein, wenn man in (18) die Gin. (34—35) berücksichtigt. Also kann man, oan 
gegründet auf Gl. (24), enh HR. stir 
und mit der gewöhnlichen Bezeichnung der kovarianten Ableitung fa m 
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m 
schreiben. Für 7/‘” resultiert infolge der Gl. 
0) Die Größen T7” und Tf” bleiben also im Sinn der kovarianten Divergenz 


und die Größen Tj", Uys, Uys und Uys im Sinn der partiellen Ableitung er- 
halten. Alle diese Größen treten bei der Betrachtung der Invarianz der 
1) Lagrange-Funktion bei Translationen auf; sie sind deshalb von der Art 
der Energie. 


Die Größe 75” wird der symmetrische (metrische) Energie-Impuls-Tensor 


2) (ohne Gravitationsenergie) genannt®). Wir nennen die Größe TY”, die ganz 

3) ähnlich der vorhergehender definiert ist, Energie-Impuls-Tensor des Gra- 
vitationsfeldes, und die Größe T}” Gesamt-Energie-Impuls-Tensor des Feld- 

‘)s systems, sowie die drei Terme Ug die entsprechenden Spinanteile der Energie 

ns der en des RR: Wir definieren den „kanonischen Energie- Im- 
puls- Quasitensor“ 


auch für alle drei Fälle. 


Die Spinenergie dient also zur Symmetrisierung des kanonischen Quasi- 
tensors. Ähnliche Begriffe sind in der speziell-relativistischen Feldtheorie 
4) benutzt; es ist einfach, sich davon zu überzeugen, wenn man ein Feld be- 

trachtet, dessen Lagrange-Funktion die zweite Ableitung des Potentials 
ebenso wie die beiden Ableitungen des metrischen Tensors nicht enthält. 
Dann haben wir für den kanonischen Quasitensor 


6 
) Der symmetrische Tensor ist in der gewöhnlichen Weise definiert. 
Da der die Gravitation enthaltende Energie-Impuls-Tensor des Feld- 
i systems verschwindet (36), haben wir für die Energiegrößen die Beziehung 
Uys + Upp = = — tra = — (typ + typ) (43) 
7 
so daß der kanonische Quasitensor des totalen Systems auch dem Erhaltungs- 
en satz genügt. In den beiden anderen Fällen findet die Erhaltung nur in solchen 
ar Gebieten statt, wo die Ableitungen des metrischen Tensors verschwinden. 
2 Die erhaltenen Ergebnisse lassen eine qualitativ neue Synthese der Re- 


sultate von Einstein sowie Lorentz und Levi-Civitä zu. Mit den letzteren 
stimmt unser Schluß über das identische Verschwinden des Totaltensors 
8) überein. Wir erhalten aber auch die genauen Erhaltungssätze für andere 
Größen (den kanonischen Quasitensor des Systems und die drei Spinenergie- 
größen), die nicht identisch verschwinden. So haben wir die Forderung Ein- 
steins befriedigt: Einstein erwiderte Lorentz und Levi-Civitä, daß die 


9 
8) D. Iwanenko u. Klass. Feldtheorie, Moskau 1951. 
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Erhaltung der verschwindenden Größe keinen physikalischen Sinn hat. Aber 
das konkrete Resultat Einsteins ist ungenau. Der Einsteinsche ,,Pseudo- 
tensor“ der Energie für das Feldsystem, das die Gravitation enthält, ist gleich 


Ss = top = — Uys, (44) 


d.h. er ist einfach der Spinanteil der Gravitationsenergie mit dem umge- 
kehrten Vorzeichen. Weil Einstein die nicht kovariante Lagrange-Funk- 
tion fiir die Gravitation benutzte, tritt noch ein kleiner Unterschied auf. Wir 
diskutieren diese Frage ausführlich im nächsten Abschnitt. 

Die aus der invarianten Lagrange-Funktion bekommenen Größen ent- 
sprechen der von Lorentz und Klein vorgeschlagenen Klasse der erhaltenen 
Größen, weil sie von zweiten Ableitungen des metrischen Tensors abhängen. 


B. Erhaltene Größen von der Art des Drehimpulses 


Der Transformation (37) folgt die Transformation 
= a? dwg; dwg = const 
als die nächst kompliziertere. Wir nennen sie ,,die verallgemeinerte Drehung“, 
weil sie im speziellen Fall zur gewöhnlichen Drehung führt: 


Jay OO, = = — ,- (46) 
Eigentlich gehört die Transformation (45) zur breiteren Klasse, die auch 
Maßstabsänderung enthält. 
Der Transformation (45) entspricht der genaue Erhaltungssatz (30), was 
man leicht bei der Vertauschung der Gl. (45) mit der Grundgleichung (13) 
prüfen kann. Mit Hilfe der Gl. (43) bekommen wir den Erhaltungssatz 


2 


für den verallgemeinerten Drehimpuls 


Be = — = 4,2 — mee. wh 
Aber fiir die beiden nichttotalen Terme haben wir Ste Erhaltungssatz: 


0 Bo a 


Aus der Symmetrie des Energie-Impuls-Tensors folgt die allgemeine Glei- 
chung 


© (Bo t 1 


die auch mit der Transformation (46) verbunden ist. Im speziell-relativisti- 
schen Grenzfall haben wir den gewöhnlichen Erhaltungssatz für OPER 
impuls: 
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AuBerdem folgt aus der Gl. (13) bei den Transformationen (45—46) die Be- 
ziehung 


4 Gab — = 0, 


ax? 


wobei die Größe 


— pr _ 


Terz 


weit 


(52) 


(53) 


in der speziellen Relativitätstheorie erhalten bleibt. Wir benennen die beiden 


Größen (51) und (53) den „Drehimpuls- Quasitensor‘‘. 


In der speziellen 


Relativitätstheorie stimmen beide so definierten Größen bis auf ein im Fall 
des geschlossenen Systems unwesentliches Divergenzglied überein *). In unserem 
Fall gilt dies nur dann, wenn die Lagrange-Funktion von den zweiten Poten- 
tialableitungen unabhängig ist, d.h. wenn in den Indizes x und ß die Größe 


m2? antisymmetrisch ist. 


Wir betonen, daß keine Vierbeine (Tetrapoden) von Belinfante°) in 
unseren Betrachtungen benutzt wurden. 


C. Multimomente 


Die Koordinatentransformationen komplizierend, kommen wir nach der 
Transformation (45) zu den Transformationen 


= x 005,5; 62%, = const, 


mia 


(54) 


die insbesondere den Übergang zu den gleichmäßig beschleunigten Systemen 
bilden. Ihnen entspricht der Erhaltungssatz (31), und wir können die Größe 


u (33) „Bimoment‘ des Systems nennen. Hier können wir auch ein anderes 
Sun einführen, das &:°" und 7“ ähnlich ist; die Größen [°° und 
Br , wie auch t3, gehören der Klasse der kanonischen Größen an. 


nennen wir NP? 


Es ist wichtig zu betonen, daß mit jedem Schritt der Komplizierung der 
Koordinatentransformationen qualitativ neue Größen auftreten, wenn in 
der Lagrange-Funktion alle Ordnungen der Ableitungen zugelassen sind. 
In unserem Fall enthalten alle Größen, die nach dem Bimoment folgen, keine 


prinzipiell neuen Glieder, sondern sind aus ts, me", N ” und x* gebildet. 


Wir nennen —M}" den ,,verallgemeinerten Spin“ 
entsteht bei Antisymmetrisierung der Indizes t und « in — M2"). Analog 
Wenn höhere 


Ableitungen in der Lagrange-Funktion auftreten, so folgen héhere Multispins. 


Die Gin. (20) und (19) ergeben 


„Bispin“ 


D. Eigenmomente (Spins) 


Us = We + Bz, 


4) W. Pauli, Rev. mod. Physics 18, 203 (1941). 
5) F. Belinfante, Physica 7, 449 Rn 


(der gewöhnliche Spin 


oder „das eigene Bimoment“. 
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(56 
ax? het. PR ) 
im eigentlichen Sinn die Spinenergie bildet, und die Göße 
ENTF y 
02° 


die Bispinenergie ist. Dabei ist die Größe %5; der Spin im eigentlichen Sinn, 
und der Summand 


der ispinarenimpuis. 


3 
$4. Lagrange-Funktionen der Felder 

Die hier gefolgerten Beziehungen (19—21) bilden den analytischen Aus- 
druck der Invarianzbedingung für die Funktion Z, dessen physikalischer Sinn 
sich nur bei Betrachtung der Euler-Lagrange- Gleichungen aufdeckt. 
Weil man mit denselben Größen eine beliebige Zahl von Invarianten kon- 
struieren kann, bilden diese Beziehungen noch keine hinreichende Bedingung 
zur Berechnung der Lagrange-Funktionen der physikalischen Felder. Zu 
diesem Zweck formulieren wir ein einheitliches Prinzip zur Aussonderung 
physikalisch sinnvoller Lagrange-Funktionen der Felder. Erst wählen wir 
L nur aus der Funktionenklasse, für die: a) das Wirkungsprinzip auf par- 
tielle Differentialgleichungen führt (nicht auf algebraische Gleichungen) 
b) L ist eine ganze rationale Funktion der Feldpotentiale (unter ihnen auch 
des metrischen Tensors g,, und g“”), sowie ihrer ersten und zweiten Ab- 
leitungen. 
Wir betrachten das Potential, den metrischen Tensor und deren Ablei- 
tungen als unabhängige Variablen. 


Jetzt formulieren wir das einheitliche Einfachheitsprinzip: 


1. Die Lagrange-Funktion Z möge die minimale Zahl der Variablen 
enthalten. 


2. L möge eine Funktion dieser Variablen von der minimalen Potenz sein. 


8 Wir geben das Massenglied (das Glied, in dem keine Potentialablei- 
tungen auftreten) auf. 


Es scheint, daß die Felder, in dessen Gleichungen es keine Massenglieder 
gibt, die einfachste Felderklasse bilden. Wir benennen solche Felder ,,kine- 
matische Felder“. Zu ihnen gehören elektromagnetische, Gravitations- und 
spinorielle Neutrino-Felder. Die Lagrange-Funktionen für all diese Felder 
kann man aus unserem Prinzip und aus den Invarianzbedingungen (19—21) 
eindeutig berechnen. Für das Neutrino-Feld muß man die Lamé-Metrik 
benutzen®). Diese Frage überschreitet aber die Grenzen dieser Arbeit. 


*) G.Rumer, Erforsch. in der 5-Optik, Moskau 1956. 
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Wir betrachten die Lagrange-Funktion L, des metrischen Feldes. Es 
” ist einfach zu prüfen, daß man nur mit dem metrischen Tensor und seinen 
56) ersten Ableitungen keine Invariante kunstruieren kann. Das ist weit bekannt, 
und wir beschäftigen uns mit diese Frage nicht mehr. 

Wir betrachten die dem Einfachheitsprinzip entsprechende Lagrange- 
Funktion, die linear von den zweiten Ableitungen des metrischen Tensors 


abhängt. Laut Bedingung (21) kann man schreiben 


Weil die Determinante des metrischen Tensors nicht verschwindet, bekommen 
wir 


Es ist leicht von g“” die einfachste Kombination zu bilden, die die Gl. (60) 
befriedigt: 
nn vr,aß 
Hieraus folgen zwei wichtige Beziehungen: 
g g ngen: 
ng = 2k (gu? gr? — gus, (62) 
zu und } 
Neo = — 2kV—g (2g%? — gr? d; — g** H). (63) 
a Hier ist Zy; ein Teil der Lagrange-Funktion, der von den zweiten Ab- 
n) leitungen des metrischen Tensors abhangt. 
% Die Ausdriicke (61) resp. (63) in (16) resp. (20) einsetzend, bekommen wir 
ö zwei Ausdrücke für 75. Der Vergleich gibt uns: 
ie a, 


2 Gra ng —kV-g ( (2 ger ght — gb” — gi") 


= WE V-9 (2 g’’ gt? gt” ger ger 


Um die mit der Einführung der Größe %55 aufkommende Unbestimmheit zu 
beseitigen, symmetrisieren wir diesen Ausdruck in den Indizes r und ß: bey 


(64) 


d Zu den Christoffel-Symbolen übergehend und die Identität 
k beriicksichtigend, bekommen wir ad 
dX, —4k aa (g"’g TH gir) (67) 
V 99 uve g g ‘ 
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. . . 
EL Hieraus kann man einfach den von den zweiten Ableitungen des metrischen 


Tensors unabhängigen Teil der Lagrange-Funktion bestimmen: 


Li = 2k (g"’ — gh) (68) 
Die ganze Lagrange-Funktion der Gravitation ist also gleich 
L, = 2k (gt? gx? — ger’) (9°* + Yunap)- (69) 


Gewöhnlich benutzt man bei der Folgerung der Euler-Lagrange- 
Gleichungen anstatt der skalaren Gravitation-Lagrange-Funktion eine 
nichtkovariante Größe, die man als „Gravitations-Lagrange-Funktion“ 
benennt. Das ist nur möglich infolge des Umstands, daß die Skalarkrümmungs- 
dichte in zwei Summanden dargestellt werden kann: ein Glied hängt von den 
zweiten Ableitungen des metrischen Tensors nicht ab, und das zweite Glied 
hat eine Divergenzform, die zu keinem Beitrag führt, weil die Variationen im 
Wirkungsprinzip an der Grenze verschwinden (wir erwähnen, daß diese 
Variationen eine ganz andere Natur haben, als unsere Ö4,). 
Die erwähnte Zerlegung ist: 


Es ist offensichtlich, daß man, ausgehend von A, keine allgemeinen Erhaltungs- 
sätze für das Gravitationsfeld folgern kann. A und A sind aber invariant bei 
linearen Koordinatentransformationen, und man kann einige erhaltenblei- 
bende Größen ermitteln. 

Unsere Lagrange-Funktion Z, muß also mit der Skalarkrümmung 


x 


_ übereinstimmen. Es ist leicht zu prüfen, daß dies bei k = + } der Fall ist. 

Die Maßeinheiten wählen wir hier so, daß, z.B. die Einsteinschen 

Gleichungen die Form 

Ruy — = —2 

annehmen. 

I Wir betrachten nun, zu welchen Möglichkeiten die Zerlegung (70) führt. 
Mit Hilfe der für A umgeschrieberien Gl. (22) kann man mühelos sehen, daß 

die Bedingung (18) im Fall des Gravitationsfeldes zur Gleichung 


j führt. Wir benutzen hier die natürliche Bezeichnung 1,3 und Wap. Die 
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folgt also mit Hilfe der Identität von Bianchi. Die Größe 11,3 ist gleich 4 Hr 
4 And. -* = A — = — (74) 


was mit der von Einstein definierten Gesamtenergie der Gravitation und der 
anderen Felder genau übereinstimmt. Der physikalische Sinn dieser Größe 
ist also klar: sie ist nur ein Teil der Spinenergie des Gravitationsfeldes. 

Es ist weiter noch leicht zu zeigen, daß sich die Zerlegung in Teile, die von 
A oder A abhängen, auch auf die drehimpulsartigen Größen ausdehnt: 

= = — Tye — tenst.c = — (75) 

0x” 

Das ist eine Folgerung der Invarianz von A bei allen linearen Koordinaten- 
transformationen, und insbesondere der verallgemeinerten Drehung. Es ist 
interessant, daß alle energieartigen Größen, die genau erhalten bleiben, eine 
Divergenzform haben (unter ihnen die von Einstein eingeführte Größe, wie 
es aus Gl. (75) offensichtlich ist), so daß man bei der Integration das Gauß- 
Theorem anwenden kann. 

Es ist natürlich, daß sich die Zerlegung in A und A auf die Energie be- 
EDER und den Bispin nicht betrifft; damit ist auch die Ungleichheit von 
und F verbunden. 

Wir schreiben jetzt die Grundgrößen aus, die das Gravitationsfeld charak- 
terisieren: Den verallgemeinerten Spin: 


Moo = —V—g — 
zu dem der Spin im eigentlichen Sinn den Anteil 
liefert. Der A-Spin ist gleich 
= —V—g Tag + Tale — gr? 65) + (78) 
Der Bispin ist schon in der Gl. (63) angeführt. Die Spinenergie ist gleich 
Ue =— V-g [g®? — g**) + Poel ites . (79) 


Diese Größe ist in die Ausdrücke für Spinenergie (im eigentlichen Sinn) 


und Bispinenergie zerlegt. Wir schreiben hier die Bispinenergie us: 


Wyo = —V-9 [( (g”® 2 g°*) (6; Iwe,aß — 


Der symmetrische (metrische) Energie-Impuls-Tensor und der kanonische 
Quasitensor haben die Form 


(76) 


gr’ g**) inne + Gorse (81) 
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-V-9 [R — (g®® — ge} (Gus,10 (82) 

+ bes + . 

Wir gehen nun zur Elektrodynamik über. Erst bemerken wir, daß wir, wenn 

wir die kovarianten und die partiellen Ableitungen für Größen von gleichem 

Einfachheitsgrad halten, eine wesentliche Mehrdeutigkeit bei der Bestimmung 

der Lagrange-Funktion des elektromagnetischen Feldes erhalten. Wir 

betrachten also laut Einfachheitsprinzip die Lagrange-Funktion, die nur 

von A, , und g,, abhängig ist. Dann bekommen wir nach der G!. (16) 

a2, 


und die GI. (20) führt zur 


imal 


art 


(83) 


Wir benennen die Feldstärke des elektromagnetischen Feldes mit H,, = 


A,u— Au,,. Die Bedingung (84) wird nur durch Funktion von H,, be- 
friedigt, und die einfachste Funktion, die keine freien Indizes hat, ist 


L, = — Hy, (85) 


Das ist die bekannte Lagrange-Funktion des Maxwell-Feldes. 
Wir geben nun die das elektromagnetische Feld charakterisierenden 


Der verallgemeinerte Spin ist: 
Mes = — 4A, (86) 


Dieser stimmt hier mit dem eigentlichen Spin iiberein. Die energieartigen 


x B x u ve 


h aß F J at). he} 
to=V—-g(H Hap 6, + 4A, H**); (88) 
‘oy | =—4 Y- g H**. (89) 


Wir schreiben hier nicht die Drehimpulse und Bimomente auf, da diese 
leicht zu bilden sind. 

Es ist auch sehr einfach, die Lagrange-Funktion des skalaren Feldes ein- 
deutig zu erhalten, die die bekannte Form aufweist. 

Es ist interessant, daß ein Feld, dessen Potential sich wie eine Skalar- 
dichte transformiert, in der speziell-relativistischen Näherung ungeachtet 
des Einfachheitsprinzips nichtlineare Gleichungen hat. Das ist auch der 
Fall für ein Spinordichtefeld. Die Nichtlinearität ist ähnlich derjenigen 
von Heisenberg, Iwanenko, Schiff, Thirring und andere?) benutzten. 

7) W. Heisenberg, Z. Naturforschg. 10a, 425 (1955); D. Iwanenko, D. Kurd- 


ge »laidse, S. Larin, Dokl. Acad. Nauk. SSR 88, 245 (1953); L.I. Schiff, Physic. Rev. 
84,1 (1951); W. Thirring, Z. a 7a, 63 (1952). 
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Es tritt aber noch ein in den Ableitungen nichtlineares Glied auf. Es ist 
natürlich, daß die Koeffizienten bei den nichtlinearen Gliedern eindeutig 
durch die Gravitationskonstante sowie Ruhmasse im verallgemeinerten Fall 
gegeben sind. 

Abschließend möchte ich meinen besten Dank Herrn Professor Iwanenko 
und auch Herrn Professor Rumer für ihr großes Interesse an der Arbeit und 
Herrn Professor Christoff für viele wertvolle Diskussionen aussprechen. 
Auch Herrn Dr. Schmutzer danke ich für kritische Diskussionen und für 
seine freundschaftliche Hilfeleistung bei der Übersetzung. 


Moskau, Physikalische Fakultät der Staatlichen Lomonossow-Universität. 


Bei der Redaktion eingegangen am 20. Juni 157. 
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Zum Phänomen der Entelektrisierung in kleinen leitenden mur 
IB Körpern in einem ausgedehnten elektromagnetischen sung 
| 
as 
Mit 5 Abbildungen ve Res 
Inhaltsiibersicht verl 
7 Der aus der Elektrostatik geläufige Ansatz für das Entelektrisierungsfeld en 
wird auf den Fall eines leitenden Körpers übertragen, dessen Abmessungen SR 
sowie Eindringtiefe des Skineffektes klein sind gegen die Wellenlängen eines en 
elektromagnetischen Hochfrequenzfeldes innerhalb und außerhalb des Körpers. 
Es werden Stromdichte und inneres Feld in Abhängigkeit vom äußeren Feld 
unter wenig spezialisierten Annahmen über den Leitfähigkeitsmechanismus Die 
im Körper berechnet. Die Ergebnisse werden für einen Germaniumkörper Bet 
in einem Temperaturbereich von 0—200° C und eine Wellenlänge von 1,24 em für 
des äußeren Feldes diskutiert. 
wol 
1. Einleitung ent: 
4 Im Rahmen einer einfachen klassischen Theorie der Zyklotronresonanz 
ist der Einfluß der Entelektrisierung auf die Stromdichte, welche sich in 
einem kleinen Körper in einem äußeren elektromagnetischen Hochfrequenz- der 
feld ausbildet, erstmalig von Dresselhaus, Kip und Kittel!)?) untersucht 
worden. Sie setzten dabei voraus, daß das elektromagnetische Feld sowohl 
im Inneren des Körpers als auch außerhalb in einem räumlichen Gebiet von 
der gleichen Größenordnung wie der Körper als näherungsweise homogen dur 
betrachtet werden kann. Ferner wurde angenommen, daß die Abmessungen 148 
des Körpers klein sind gegen die Eindringtiefe des Skineffektes, und daß 
dem Leitfähigkeitsmechanismus im Körper eine von der Energie unabhängige 
Relaxationszeit entspricht. Die 
Angesichts der Bedeutung, welche den Zyklotronresonanzmessungen sowie Lei 
Messungen der Leitfähigkeit von Halbleitern im Zentimeterwellengebiet zu- 
kommt, erschien es lohnend, in einer ersten Arbeit?) den Ansatz von Dressel- 
haus, Kip und Kittel für das Entelektrisierungsfeld aufzugreifen und das De 
Problem unter allgemeineren Voraussetzungen über den Leitfähigkeits- mit 
mechanismus zu betrachten. In der vorliegenden Arbeit wollen wir diese 
Untersuchungen fortsetzen und präzisieren. Da 
1) G. Dresselhaus, F. Kip u. C. Kittel, Physic. Rev. 98, 369 (1955). a 


a 2) G. Dresselhaus, F. Kip u. C. Kittel, Physic. Rev. 100, 618 (1955). 
er 3) H. Stolz, Ann. Physik (6), 19, 394 (1957). 


en 
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2. Die Entelektrisierung 


Zur Messung von Zyklotronresonanzen !)2)*) und der dielektrischen Eigen- 
schaften von festen Körpern im Zentimeterwellengebiet®) ist eine Anordnung 
verbreitet, bei welcher die zu untersuchende Substanz in Form eines kleinen 
Probekörpers in einen Hohlraumresonator eingebracht und dessen Verstim- 
mung und Dämpfung beobachtet wird. Wir setzen voraus, daß die Abmes- 
sungen des Probekörpers klein sind gegen die Eindringtiefe des Skineffektes 


sowie gegen A, und Et , wo A, die Wellenlänge des Feldes im leeren Resonator 
e’ 


und e’ die Dielektrizitätskonstante des Probekörpers ist. Dann können wir 
das Feld €,, welches in dem vom Probekörper eingenommenen Gebiet des 
Resonators herrschen würde, wenn man den Probekörper entfernte, als 
näherungsweise homogen annehmen. Den Körper denken wir uns aus einer 
verlustfreien Grundsubstanz mit der Dielektrizitätskonstanten e, bestehend, 
in welcher sich Leitungselektronen mit der Gleichstromleitfähigkeit o° be- 
wegen können. Die Leitungselektronen werden einen Beitrag &, zur Dielek- 
trizitätskonstanten leisten. Wir schreiben für die gesamte Dielektrizitäts- 
konstante des Körpers 


=e, + &. (1) 
Die Leitfähigkeit o’ im Wechselfeld können wir in bekannter Weise in die 


Betrachtungen einbeziehen, indem wir eine komplexe Dielektrizitätskonstante 
für das Leitungselektronengas einführen: 


420’ 


wobei w die Kreisfrequenz des Hochfrequenzfeldes ist. In dieser Schreibart 
entspricht der Dielektrizitätskonstanten ¢, der Imaginärteil 


o=0’+io”. (5) 


Wir betrachten einen Probekörper, in welchem sich die Entelektrisierung 
durch einen konstanten skalaren Entelektrisierungsfaktor L beschreiben 
läßt und setzen für das elektrische Feld im Inneren des Körpers 


E=€,—LF. (6) 
Die Polarisation ® besteht: aus einem Beitrag der Grundsubstanz und einem der 
Leitungselektronen : 


(xo + Xe) ) &. u de 


Der letztere ist wegen (2) j im allgemeinen nicht mit € in Phase. Wir schreiben 
mit Dresselhaus, Kip und Kittel 


E=—ent. (8) 

Dabei ist Me die Leitungselektronenkonzentration und r die über alle Elek- 
— mie 
4) R.N. Dexter, H. J. Zeiger u. B. Lax, Physic. Rev. 104, 637 (1956). an 


5) H.B.G. Casimir, Philips Res. Rep. 6, 162 (1951). 
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tronen im em? gemittelte Verschiebung eines Elektrons gegenüber dem Gitter 
unter der Einwirkung des Feldes ©. Mit (6), (7) und (8) erhalten wir 


3. Die Boltzmann-Gleichung 


Unser Ziel ist es, die Stromdichte zu berechnen, welche sich unter der 
Einwirkung des Hochfrequenzfeldes im Inneren des Probekörpers ausbildet. 
Dazu müßten wir das effektive elektromagnetische Feld Egr, Der kennen, 
welches auf ein einzelnes Leitungselektron wirkt. Dieses Feld ist sicher eine 
Funktion des Feldes €,, $,, enthält aber außerdem in höchst komplizierter 
Weise die Wechselwirkung der Elektronen untereinander sowie mit dem 
Gitter der Grundsubstanz. Wir wollen zur Vereinfachung €; mit dem durch (6) 
definierten makroskopischen Feld € identifizieren: 


L&; 
bzw. mit (9) und (11) 


Cor = (E, + Len, t). (12) 


An dieser Stelle muß darauf hingewiesen werden, daß wir in Gl. (8), welche zum 
Ausdruck bringt, daß en, 1 ~ e'®! ist, das Ergebnis der folgenden Rech- 
nung bereits vorausgenommen haben. In dem Ansatz (9) bzw. (12) ist diese 
Vorwegnahme jedoch nicht enthalten. 

Wir nehmen an, daß €, nur eine von Null verschiedene Komponente E, 
in x-Richtung hat und legen für das Folgende die in üblicher Weise lineari- 
sierte Boltzmann -Gleichung zugrunde: 


setzen voraus, daß für die Änderung der Verteilung durch ,,StéBe eine Re- 
laxationszeit 7 (p) existiert. 
Mit dem Ansatz 


ie 


erhalten wir aus (13) wort 
of° 
_ £—en (E, () + Len, 2) op 

—en, (p, t) p, dp (14) 


ergib 


Wir 


Die 


tion 
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schr 
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ergibt sich durch nochmaliges Differenzieren nach t: is ie 
1% . op 
Wir setzen nun an: shite Ha 


h (p, th m eet, 


‘0 
m* h(p,t Pr dp 
Die Auflösung dieser Gleichung nach h (p, t) ist elementar und liefert Me ” we 
24 
T 1 
Mit (14) und (15) erhalten wir daraus die Stromdichte Zn PE 
1 T of 3 
„eo: Li T op Px - nE, (t) 
1+ 4+— p,d 
m* w l+iworöop, 
(17) 


An (17) knüpfen wir folgende Feststellungen: Die Stromdichte j ist propor- 
tional dem in (10) definierten Feld €; = 7 &,: 


j = E, (18) 


wobei man o, unmittelbar aus (17) entnimmt. Damit ist die in (8) ausge- 
sprochene Zeitabhängigkeit j ei”! bestätigt und die Möglichkeit der Be- 
schreibung der Verhältnisse im Probekörper durch die komplexe Dielektri- 
zitätskonstante &’+ie” bzw. Leitfähigkeit o’ + io’ erwiesen. Diese ist 
offenbar wesentlich an die Linearisierung der Boltzmann-Gleichung und 
die Voraussetzung der Existenz einer Relaxationszeit gebunden. 
Ignorieren wir die p-Abhängigkeit von r in den beiden Integralen in (17), 
so erhalten wir 
e? 
mit 
(1-5), 


2 
ww“ 


Der Ausdruck (19) führt auf die von Groschwitz und Siebertz®) angege- 
nen Ausdrücke für o’ und x; und läßt sich aus den von Dresselhaus, Kip 


6) E. Groschwitz u. R. Siebertz, Z. Naturforschg. 11a, 482 (1956). 
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und Kittel!)2) verwendeten Formeln für die Zyklotronenresonanz gewinnen, 
wenn man das äußere statische Magnetfeld fortläßt. Er ergibt sich in guter 
Näherung für stark entartete Leitungselektronen unmittelbar aus (17). Den 
an anderer Stelle?) abgeleiteten Ausdruck für die Stromdichte erhält man, 
wenn man die p-Abhängigkeit von r im Nenner der Formel (17) vernach- 
lässigt. Dann gewinnt der Bruch die Gestalt 


l+ior 
l+ior' 


Zieht man ihn unter das verbleibende Integral in (17), so bekommt man bis 
auf die Wahl der Integrationsvariablen unsere Formel (10) in?). 


4. Der mittlere Energieverlust im Probekörper 


Die für die Energieverluste in einem Hohlraumresonator charakteristische 
Größe ist die „Cüte‘‘ des Resonators: 


do" 


Dabei ist w, eine Resonanzfrequenz und dw die Halbwertsbreite der Re- 

sonanzkurve, welche von einem Detektor mit quadratischer Kennlinie ange- 

zeigt wird )’). Es ist andererseits 


W 


wo W die im zeitlichen Mittel gespeicherte Energie im Seemann und W das 
zeitliche Mittel des Energieverlustes ist. Bei der Berechnung des Energie- 
verlustes hat man für die Verluste im Probekörper das Integral eee 


t= [ol => fo’ 2 (20) 
Probe e 

zu berechnen. Dies stößt auf die Schwierigkeit, daß man das Feld © im Inneren 


des Probekörpers nicht kennt. Man kann daher versuchen, mit einer geeignet 
zu definierenden Leitfähigkeit 6 = 0’ + i0’’ die Verluste durch 


1 
(22) 


G; 

Wir werden nun zeigen, daß o’ mit of [vgl. (17)] identisch ist. Wir 

haben nach (17) mit L= 0 


e? 
(24) 


7) F. Horner, T. A. Taylor, R. Dunsmuir, J. Lamb, W. Jackson, JEE Part 
III, 98, 53 (1946). 
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Andererseits ist wegen 


mit (10) 


E; 


olgt 


0, 


= 
und der Vergleich mit (25) liefert or 
(27) 


Diese Beziehung zwischen o und gy, ist in®) dadurch, daß in (26) im Nenner 
irrtümlich 07, statt o, eingesetzt wurde bzw. im Nenner von (10) fälschlich 
4, statt y,, unrichtig angegeben worden. Die Beziehung (27) ist einleuchtend: 
Die Relation zwischen Stromdichte und Feldstärke € im Körper darf nicht 
von seiner Form und Orientierung im äußeren Feld ©, abhängen. 
Gehen wir mit (27) in (24) ein, so folgt ~ 


Zerlegen wir dieses in Real- und Imaginärteil, so ergibt sich 
o 


und der Vergleich mit (28a) liefert 6) pak 


Wir fassen die Ergebnisse dieses Abschnittes folgendermaßen zusammen: 


23* 


en, 3 
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Wertet man die gemessene Güte Q des Resonators bei anwesendem Probe- 


körper unter Verwendung des Ausdruckes 
w= Par 


für die Verluste in der Probe aus und bestimmt damit oj als Funktion der 
gemessenen charakteristischen Größen des Resonators und des Probekörpers, 
so hat man dieses gemessene G7,gem, mit dem theoretischen 


(30) 


| 
zu vergleichen, worin 


5. Diskussion der Ergebnisse 


ist. 


Charakteristisch für die für den Energieverlust im Probekörper entschei- 
dende Größe oj, ist zunächst, daß sie nur die phänomenologischen Größen L, 
7, @,0’,c’' enthält und daß in ihr die Ladungsträgerkonzentration n, nicht 
explizit auftritt. Dieses Ergebnis kann insofern kaum überraschen, als wir 
zur Beschreibung der Entelektrisierung die phänomenologische Gl. (6) zu- 
grunde gelegt haben sowie die ebenfalls phänomenologische Beziehung (10) 
zwischen dem inneren und dem äußeren Feld. Im Vergleich zu den Resultaten 
von Dresselhaus, Kip und Kittel!)?) bzw. Groschwitz und Siebertz®) 
ergibt sich jedoch der unter Umständen wichtige Schluß, daß für das Auf- 
treten eines Entelektrisierungseffektes bei 
=1 und Lig (32) 
[477 
Elektronenkonzentration und Beweglichkeit wesentlich sind. Bei Dressel- 
haus, Kip und Kittel ist die für die Entelektrisierung entscheidende Größe 


die Frequenz @, =/ - n,Ln, und die Entelektrisierung wird wesentlich, 

wenn 

o, = L nn, = @” (33) 


ist. Wir werden an einem Beispiel erkennen, daß das Kriterium (33) nicht aus- 
zureichen braucht zur Entscheidung, ob bei einer gegebenen Probe mit einer 
Beeinflussung der Ergebnisse von Verlustmessungen durch die Entelektri- 
sierung gerechnet werden muß oder nicht. Da sich aus unseren Formeln (30) 
und (31) die Ergebnisse von!)?) für hohe Entartung der Leitungselektronen 
bzw. energieunabhängige Relaxationszeit 7 gewinnen lassen, wird man damit 
rechnen müssen, daß besonders im Falle nicht entarteter Leitungselektronen 
daß Kriterium (33) unzureichend sein wird. Es braucht kaum erwähnt zu 
werden, daß mit o’ und o” natürlich in erster Linie die Elektronenkonzen- 
tration n, für eine Beeinflussung der Verlustmessungen durch die Entelek- 
trisierung maßgebend sein wird. 
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In Abb. 1 haben wir die Größe 4 in Abhängigkeit von der Temperatur 


dargestellt fir Germanium mit dem spezifischen Widerstand 0° = 56 Ohm - em 
bei 20°C, dessen Gleichstromleitfähigkeit zwischen 0 und 100° C von Klin- 
der ger®) gemessen wurde. Zwischen 100°C und 200°C haben wir Klingers 
u Werte für einen konstanten Bandabstand AE = 0,77 eV extrapoliert. Nach 
(31) haben wir o unter Verwendung der empirischen Werte von o® von Klinger 
und der von Donovan und March?) angegebenen Auswertung des Inte- 
grals für rw et berechnet. Für die effektive Masse der Elektronen haben 


(30) wir den Wert m, = 0,12 m, benutzt. Alle Rechnungen gelten für eine Wellen- 

länge von 1,24 em des Hochfrequenzfeldes. 
31) 
1ei- 

cht 
wir 4 
zu- 03 
10) 

0, 
ten 
2°) 01 
- 280 300] 320] 34 360 380 400 420 440 460 Bae 

32) Abb. 1 Abb. 2 
el- In Abb. 1 bemerken wir bei etwa 380° K einen starken Abfall von oj 
Be und damit des Verlustes. Bei dieser Temperatur wird 
ch, 


33) während w, ~ w bereits bei etwa 320° eintritt. Oberhalb 450° K verschwindet 
o;/o, praktisch und die Verluste streben gegen Null. Mit diesem Verhalten 
m von oj, vergleicht man den Verlauf von Re (E/E,) und Im (E/E,) in Abb. 2 
arg und 3. Man erkennt, daß bei etwa 380° Re (E/E,) verschwindet und Im (E/E,) 


tri- ein steiles Maximum durchläuft. Oberhalb 450° geht E insgesamt gegen Null. 
30) In diesem Temperaturbereich dringt das Feld immer weniger stark in den 


Probekörper ein, er verhält sich immer mehr wie ein metallischer Körper 
nit (dieser Effekt darf nicht mit dem Skineffekt EEE werden). Wir 


er merken an, daß der Beziehung 
=0 (34) 


) Y. Klinger, Physic. Rev. 92, 509 (1953). 
) B. Donovan u. N. H. March, Proc. physic. Soc. B 69, 528 (1956). REISE 
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die Relation entspricht: 


Das bedeutet, daß der Effekt der Entelektrisierung sein Maximum erreicht 
im Gebiet negativer Dielektrizitätskonstanten e’. Wie wir gesehen haben, 
verhält sich der Probekörper aber schon in Gebieten mit noch positiver Di- 
elektrizitätskonstante anders als bei Abwesenheit der Entelektrisierung zu 
erwarten wäre. 


) 


=e te —=1-— 


w 
; 280 300 320 340 360 380 400 420 440 460 T°K 006 006 O20 030 00 060 080100 WO 200 30 40 60 80700 40 
ee Abb. 3 Abb. 4 


Für manche Fälle mag es interessant sein, die Frequenzabhängigkeit des 
Entelektrisierungseffektes zu kennen und insbesondere die Größenordnung 
des Frequenzbereiches, in welchem er auftritt. Wir haben daher in Abb. 4 
die charakteristische Größe 

1 


fiir dieselbe, den Abb. 1— "3 
zugrunde liegende Germanium- 


Mr probe für eine Temperatur 
12h von 370°K (dem entspricht 
10+ eine Elektronenkonzentration 

n, = 1,6- 1014 cm?) in Ab- 
|  hängigkeit vom Frequenz- 
067 verhältnis w/w) (w» = 1,52: 
aut 10-1 sec! ~ 1,24 cm) aufge- 


a tragen. Die Abb. 5 ist eine 
) 40 analoge Darstellung für 
-021 W Es muß betont werden, 
daß die Abb. 1—5 nur orien- 
en tierenden Charakter besitzen, 
Er da sie das Ergebnis von Rech- 
ua sind, welche nur die 
Elektronen-Gitterwechselwirkung in r berücksichtigen im Anschluß an?). 
Don ovan - March kommen beim Vergleich ihrer RER mit entspre- 
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chenden Experimenten zu dem Schluß, daß man einen komplizierteren Störungs- 
mechanismus, als dem 7 ~ e~!-Gesetz zugrunde liegt, annehmen muß, wenn man 
die Formel (31) zur Berechnung von o mit Erfolg verwenden will. Die generelle 
Temperaturabhängigkeit von o’, 0”, welche für die Abb. 1—3 wesentlich ist, 
wird jedoch durch (31) richtig wiedergegeben, die Diskrepanzen, auf welche 
Donovan und March hinweisen?) treten erst bei der Bildung von - = 
auf, deren Größenordnung jedoch ebenfalls richtig aus (31) folgt. Ferner 
haben wir die Temperaturabhängigkeit der Größe n, welche über die Dielek- 
trizitätskonstante e, der Grundsubstanz von der Temperatur abhängt, ignoriert, 
da es zur Zeit wohl keine Information über die Funktion ¢, (T) für Germanium 
gibt. 


6. Zusammenfassung _ 


Wir haben gesehen, daß es möglich ist, ohne untäinibere Annahmen über 
den Leitfähigkeitsmechanismus die Verluste in einem kleinen leitenden Probe- 
körper in einem äußeren Hochfrequenzfeld durch im Prinzip leicht zugäng- 
liche Größen darzustellen. Dabei konnten wir den Einfluß der Entelektri- 
sierung erfassen und feststellen, daß er mit abnehmender Dielektrizitäts- 
konstanten e’ wächst und schließlich bei negativem e’ dazu führt, daß das 
elektromagnetische Feld überhaupt nicht mehr in den Körper eindringt. Wir 
haben eine einfache Formel angegeben, welche es erlaubt, bei bekanntem 
Entelektrisierungsfaktor, bekannter Leitfähigkeit und Dielektrizitätskon- 
stanten für eine gegebene Frequenz für einen kleinen Probekörper den Einfluß 
der Entelektrisierung auf die im Hohlraumresonator meßbaren Verluste ab- 
zuschätzen. 

Abschließend möchte der Verfasser Herrn Dr. B. Mühlschlegel für an- 
regenden Gedankenaustausch und förderndes Interesse an der Entstehung 
dieser Arbeit herzlich danken. 


Berlin, Institut für theoretische Physik der Humboldt-Universität. 


SEE Bei der Redaktion eingegangen am 4. September 1957. IS 
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Über die elektrische Leitfähigkeit von Blei(II)-oxyd 


ER Zusätzen höherwertiger gitterfremder Kationen 
he Von HorstGrunewald und W. Neumann 


Inhaltsiibersicht 


_ Die vorliegende Arbeit behandelt die Abhängigkeit der elektrischen Leit- 
fähigkeit von gesintertem Blei(II)-oxyd (PbO) vom Gehalt an Chromoxyd, 
Neodymoxyd und Cerdioxyd. Der Gehalt an Fremdoxyd variiert dabei 
zwischen 0 und 2,0 Mol-%, das untersuchte Temperaturgebiet erstreckt sich 
von 200° C bis 500° C. Hierbei wird gefunden, daß bei einem Zusatz von bis 
zu 2 Mol-% Cr,O, oder CeO, zum PbO praktisch keine Veränderung der elek- 
trischen Leitfähigkeit gegenüber dem reinen Bleioxyd hervorgerufen wird. 
Lediglich beim Zusatz von Neodymoxyd ist eine schwache aber stetige Ab- 
nahme der elektrischen Leitfähigkeit mit steigendem Nd,O,-Gehalt feststell- 
bar, wenn man von der 200° C-Isothermen absieht. 

Auf Grund der im Vergleich mit anderen Mischphasen fast vollständigen 
Unabhängigkeit der elektrischen Leitfähigkeit des Bleioxyds vom Gehalt 
an böherwertigen gitterfremden Kationen im Temperaturbereich zwischen 
200° C und 500° C ist mit Sicherheit anzunehmen, daß hier ein anderer Fehl- 
ordnungsmechanismus vorliegen muß als die ausschließliche Ausbildung von 
Substitutionsstörstellen im Kationenteilgitter und einer äquivalenten Zahl 
von Elektronenstörstellen. 


Zur Vereinfachung der Technologie bei der Herstellung von Halbleiter- 
widerständen für die verschiedensten Anwendungen in der Elektrotechnik 
ist es zweckmäßig, von der Verwendung des anreduzierten Magniesium-Titan- 
Spinells (Urdox) abzugehen und statt dessen oxydische Mischphasen als Aus- 
gangsmaterial zu benutzen. Dieses Vorhaben hat zur Voraussetzung, daß die 
elektrischen Eigenschaften, insbesondere die Leitfähigkeit und der Tem- 
peraturkooeffizient der Leitfähigkeit für eine möglichst große Zahl oxydischer 
Werkstoffe bekannt sind. Zu den elektrischen Eigenschaften, deren Kenntnis 
bei einer technischen Anwendung des betreffenden Halbleiterwerkstoffes not- 
wendig sind, gehört ferner die Abhängigkeit der elektrischen Leitfähigkeit des 
Materials von mengenmäßig sehr geringen Zusätzen an anderswertigen Kat- 
ionen. Nach der Wagner-Schottky-Hauffeschen Fehlordnungstheorie 
kann man durch eine Veränderung der Zusätze nach Menge und Art die elek- 
trische Leitfähigkeit des Halbleiters definiert verändern!) und gleichzeitig 
produktionstechnisch etwaige Unterschiede in der Beschaffenheit des Aus- 
gangsmaterials ausgleichen. 


') H. Grunewald, Wiss. Z. Päd. Hochsch. Potsdam 1, 101 (1955). 
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OF Einleitung und Problemstellung 


Ziel der vorliegenden Arbeit war, in Analogie zu den an anderen Metall- ss 
oxyden angestellten Untersuchungen!) durch Messung der Abhangigkeit ae 
der elektrischen Leitfähigkeit vom Zusatz an anderswertigen gitterfremden ß 
Kationen den Stromleitungsmechanismus in Blei(II)-oxyd aufzuklären. Die 
Messungen wurden zunächst mit höherwertigen gitterfremden Kationen, und 7 
zwar mit Cr’-Ionen in Form von Cr,0,, Nd’"-Ionen in Form von Nd,O, +1 
und Ce’**’-Ionen in Form von CeO, durchgeführt. Unter der Annahme, daß R 
im PbO eine reine Elektronenleitung (Elektronendefektleitung, Elektronen- 
überschußleitung oder Gittereigenhalbleitung) vorliegt, ergeben sich für den 
Einbau der höherwertigen Fremdionen in das Bleioxydgitter aus der Wagner- 
Hauffeschen Fehlordnungstheorie®) die folgenden Gleichungen: 

Für den Fall, daß PbO ein Elektronenüberschußleiter ist: 


bzw | 
wa 
Nd,0, = 2Nd@ (Pb) + 2 © + 2 PbO +} Of 


bzw. 
CeO, = Ce@™ (Ph) + 26 + PbO +30. rt 


Für den Fall, daß PbO ein Elektronendefektleiter ist: ! 
0130, +2 = 2Cr@ (Pb) + 2PbO + 109) (2a) 


(1e) 


bzw. 
Nd,O, + 2@ = 2Nd@’ P(b) + 2 PhO + 40” (2b) 
bzw. 
CeO, + 2 © = Ce®™ (Pb) + PbO + 40M. (2e) 
Für den Fall, daß PbO ein Gittereigenhalbleiter ist: Kann 23 
Cr,0, = 2Cr@' (Pb) + 20 + 2 PbO + | (3a) 
Nd,0, = 2Nd@' (Pb) +2© +2 PbO + 3 OY (3b) 
bzw. 
‘ge CeO, = Ce@ (Pb) + 26 + PbO + 30%. u 
In allen diesen Fällen muß also eine Beeinflussung der elektrischen Leit- Be 


fähigkeit des PbO durch die ins Gitter eingebauten höherwertigen Fremdionen 
eintreten. Und zwar muß sich im Falle der Elektronenüberschußleitung analog 
zum Zinkoxyd®) und im Falle der Gittereigenhalbleitung analog zum Kupfer- ' 
oxyd?) wegen der Erhöhung der Konzentration an. quasi-freien Elektronen . 
mit steigendem Fremdoxydgehalt eine Zunahme der elektrischen Leitfähigkeit 

ergeben, während sie im Falle der Elektronendefektleitung analog zum 
Nickeloxyd’) wegen der Verminderung der Konzentration an Defektelek- 

tronen mit steigender Fremdionenkonzentration abnehmen muß. 


2) K. —, on H. Grunewald, Z. physik. Chem. 198, 248 (1951). 

3) H. Grunewald, Ann. Physik (6) 14, 121 (1954). 

4) H. Grunewald, Ann. Physik (6) 14, 129 (1954). 

5) K. Hauffe, ‚Reaktionen in und an festen Stoffen‘, Springer Verlag Berlin- 
Göttingen-Heidelberg 1955. 

6) K. Hauffe u. A. L. Vierk, Z. physik. Chem. 196, 160 (1950). 

7) K. Hauffe, Ann. Physik (6) 8, 201 (1950) 27, 
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Herstellung der Proben und Durchführung der Messungen 


Als Ausgangsmaterial wurde Bleioxyd DAB 6, Chromoxyd p.a., Cer- 
dioxyd p.a. und Neodymoxyd 99,8%, verwendet. Das reine Bleioxyd bzw. 
die im Mörser innig vermischten Bleioxyd/Fremdoxydmischungen wurden 
in einer Handspindelpresse zu Tabletten von 12 mm Durchmesser und etwa 
3 mm Höhe verpreßt. Der Preßdruck wurde subjektiv bei allen Pastillen in 
etwa der gleichen Größe gewählt. Anschließend wurden die Pastillen bei einer 
Temperatur von 600° C 4,5 Stunden lang in Luft von Atmosphärendruck ge- 
sintert. Sie erhielten dabei eine für die weitere Verwendung ausreichende 
Festigkeit. Die Anwendung einer Schutzgasatmosphäre beim Sintern er- 
schien nicht notwendig, da nach Hollemann- Wiberg’) eine Aufoxydation 
des PbO zu Pb,O,, die bei etwa 500° C einsetzt, bei Temperaturen über 550° C 
durch eine Zersetzung des Pb,O, unter Angabe von Sauerstoff wieder rück- 
gängig gemacht wird. 

Die Stirnflächen der Pastillen wurden dann plangeschliffen und mit 
Silber kontaktiert. Die Messung erfolgte in Luft von Atmosphärendruck in 
kleinen Kammeréfen zwischen Platinkontakten. Die L.eitfähigkeiten der 
einzelnen Proben gleicher Zusammensetzung stimmten dabei sehr gut mit- 
einander überein. 

Meßergebnisse und Auswertung 
In den Abb. 1 bis 3 ist die Abhängigkeit der elektrischen Leitfähigkeit 


von Bleioxyd vom Gehalt an Cr,O,, CeO, und Nd,O, dargestellt. Wie man 
aus den Abbildungen ersieht, bleibt in allen Fällen dienach den Gin. (1) bis (3) 


300°] 
tog 50” 004 5.09" 
0° 


4, 51 
300% 

200°C 
05 10 15 20 heil: 05 10 15 20 


ish. ry. 


—— 4% ——> Mol % Nd,0; 
Abb. 1. Abhängigkeit der elek- __ $n Abb. 2. Abhängigkeit der elektrischen 


trischen Leitfähigkeit des PbO vom tes Leitfähigkeit des PbO vom Gehalt an 
Gehalt an Cr,O, Nd,0, 


8) Hollemann-Wiberg, ,,Lehrbuch der anorganischen Chemie‘, Walter de 
& 1956, S. one. 
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zu erwartende starke Leitfähigkeitsänderung aus. Bei den Proben mit Chrom- 
oxydzusatz findet man, wie Abb. 1 zeigt, bei einer Erhöhung des Cr,O,- 
Gehaltes ein unregelmäßiges Verhalten der elektrischen Leitfähigkeit, wobei 
die Schwankungen jedoch den Faktor 2 nicht übersteigen. Insgesamt ge- 
nommen zeigen die Leitfähigkeitskurven mit wachsendem Cr,0,-Gehalt 
eine schwach ansteigende Tendenz. Dagegen beobachtet man an den Proben 
mit CeO,- und Nd,O,-Zusatz (Abb. 2 und 3) ein mit wachsendem Fremdoxyd- 
gehalt zunehmendes schwaches Absinken der elektrischen Leitfähigkeit. Es 
liegen also offenbar ähnliche Verhältnisse vor, wie sie von Grunewald 
bereits am TiO, mit niederwertigeren Zusätzen (NiO, Al,O,, Ga,0,) gefunden 
worden sind), 


| 


fog 


/ 
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Abb. 3. Abhängigkeit der elek- Abb. 4. Abhängigkeit der elektrischen 


trischen Leitfähigkeit des PbO vom 
Gehalt an CeO, 


Leitfähigkeit von reinem Blei(II)-oxyd 
eh von der Temperatur 


In Übereinstimmung mit dem Befund der Unabhängigkeit der elektrischen 
Leitfähigkeit des Bleioxyds vom Gehalt an höherwertigen gitterfremden 
Kationen steht auch die Tatsache, daß bei sämtlichen untersuchten Misch- 
phasen die Temperaturabhängigkeit der elektrischen Leitfähigkeit über- 
einstimmt. Es ergibt sich in allen Fällen der in Abb. 4 dargestellte Verlaui. 

Eine genaue Deutung der in den Abh. 1 bis 3 dargestellten experimentellen 
Ergebnisse kann heute noch nicht gegeben werden. Im Falle des Zusatzes 
von Chromoxyd könnte man annehmen, daß als Folge der sehr unterschied- 
lichen Ionenradien (Cr = 0,65 A Pb’ = 1,32 A) das Chrom nicht substi- 
tutionell in das PbO-Gitter unter Bildung von Cr@' (Pb)-Störstellen einge- 
baut wird, sondern daß sich einfach ein mechanisches Gemenge von PbO 
und Cr,O, ergibt. Dann müßte jedoch unter allen Umständen ein Zusatz von 
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Nd,O, zu der zu erwartenden Leitfähigkeitsänderung führen, da hier ausge- 
speochen gute Substitutionsmöglichkeiten als Folge der sehr gut überein- 
stimmenden Ionenradien (Pb = 1,32 A; Nd = 1,33 A) besteben. 

Aussichtsreicher erscheint daher die Annahme, daß die Erzeugung von 
Substitutionsstörstellen vom Typus Cr@' (Pb) bzw. Nd@* (Pb) nicht nur 
zu einer Veränderung der Konzentration der Elektronenstörstellen sondern 
zusätzlich noch zu einer Veränderung der Konzentration der durch die Eigen- 
fehlordnung im PbO-Gitter hervorgerufenen Ionenfehlordnung fiibrt. 

Nimmt man beispielsweise an, daß das reine PbO ein Elektronenüber- 
schußleiter ist, so kann also eine Eigenfehlordnung entweder durch Einbau 
von überschüssigen Pb’-Ionen auf Zwischengitterplätzen gemäß 


- 
PbO = PhO” +29 4 

oder aber durch Ausbildung von Sauerstoffionenleerstellen gemäß vice 
PbO = 00" +26 +30 (5) 


vorliegen. 

Nehmen wir zunächst an, daß im Bleioxyd ein Eigenfehlordnungsmechanis- 
mus nach Gl. (5) vorliegt nnd bauen wir höherwertige Fremdionen, z. B. 
Nd**’-Ionen, ay das PbO-Gitter ein, so ist denkbar, daß durch die in den Gin. (1) 
auftretenden }O¥ eine entsprechende Zahl von Snanrstetienineneteien 
im Grundgitter gemäß 


00” + 40 + 20 = Null (6) 


besetzt werden, wobei gleichzeitig pro Sauerstoffionenleerstelle zwei quasi- 
freie Elektronen vernichtet werden. Unter der Voraussetzung einer Rekom- 
bination gemäß Gl. (6) würden sich als Bruttoreaktionsgleichungen fiir den 
Einbau höherwertiger gitterfremder Kationen in PbO se se statt der 
Gin. (1a) bis (1e) die Gleichungen 


Cr,0; + 00” = 2Cr@' (Pb))+ 2 PbO (7a) 
Nd,O, + = 2Nd@ (Pb) (7b) 
bzw. 
CeO, + 00° = Ce@” (Pb) + PbO (Ze) 


ergeben. Mit diesen Gleichungen könnte die Unabhängigkeit der elektrischen 
Leitfähigkeit des Blei(II)-oxyds von Gehalt an höherwertigen gitterfremden 
Kationen verstanden werden. 

Legt man als Eigenfehlordnung des PbO die Gl. (4) zugrunde, so läßt 
sich die Konstanz der Leitfähigkeit bei wechselndem Gehalt an höherwertigen 
gitterfremden Kationen deuten durch die Annahme einer entsprechenden 
Verminderung der Konzentration der Bleiionen auf Zwischengitterplätzen 
gemäß 


PbO” + 40% +20 = PbO, Br (8) 

so daß sich jetzt als Bruttoreaktionsgleichungen errang 

Cr,0, + PbO” = 2Cr@’ (Pb) + 3 PbO (9a) 
Nd,O, + = 2Nd@° -3 PbO 
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+ GeO, + PhO" = Ce@" (Pb) + 2PHO (9e) 
ergeben. 
von Nimmt man andererseits an, daß das reine PbO ein Elektronendefekt- 
nur leiter ist, so müßte also nach den Gin. (2a) bis (2c) eine Abnahme der elek- 
ern trischen Leitfähigkeit mit wachsendem Gehalt an höherwertigen gitter- 
en- fremden Kationen eintreten, die genauso wie die Leitfähigkeitszunahme im 
Falle einer Elektronenüberschußleitung nach den an anderen Mischoxyd- 
er- systemen gesammelten Erfahrungen®) mehrere Zehnerpotenzen betragen 


müßte. Auch in diesem Falle kann eine Deutung der Konstanz der elektri- 
schen Leitfähigkeit des Bleioxyds bei wechselndem Gehalt an höherwertigen 
gitterfremden Kationen nur durch die Annahme erfolgen, daß die durch die 
(4) höherwertigen Fremdionen hervorgerufene Änderung der Defektelektronen- 
konzentration durch eine Veränderung der Ionenstörstellenkonzentration 
ausgeglichen wird. Im Falle einer Elektronendefektleitung des reinen Blei- 


au 


(5) oxyds kommen für die Eigenfehlordnung die beiden Fehlordnungsmechanismen 
409 =00”+2® (10) 
B. durch die Ausbildung von Sauerstoffionen auf Zwischengitterplätzen, bzw. 
1) 40% = PbO” + 2 + PbO (11) 
en 


durch Ausbildung von Bleiionen-Leerstellen in Betracht. 

Legt man zunächst den Eigenfehlordnungsmechanismus nach Gl. (10) 
6) zugrunde, so können durch die in den Einbaugleichungen Gln. (2a) bis (2c) 
auftretende $0 zusätzliche Sauerstoffionen auf Zwischengitterplätzen er- 
zeugt werden, so daß sich bei einer Überlagerung der beiden Effekte als 
Bruttoreaktionsgleichungen 


er Cr30, = 2Cr@* (Pb) +0C”+2PbO (12a) 


bzw. 
(12b) 


Nd,O, = 2 (Pb) + 00” + 2 PbO 


a) bzw. 
CeO, = Ce@™ (Pb) + 00” + PbO (12e) 
») ergeben. 
Nimmt man ale als Eigenfehlordnungsmechanismus im reinen Blei- 
’) oxyd eine Elektronendefektleitung durch Ausbildung von Bleiionenleerstellen 
n gemäß Gl. (11) an, so könnten durch die beim Einbau höherwertiger gitter- 
n fremder Kationen in das PbO-Gitter entsprechend den Gln. (2a) bis (2c) 


entstehenden 40% unter Erweiterung des Gitters nach außen neue Blei- 
t ionenlerstellen entstehen, so daß sich in diesem Falle als Bruttoreaktions- 
gleichungen für den Einbau höherwertiger gitterfremder Kationen in das 
Bleioxydgitter die Gleichungen 


ee Cr,0, = 2Cr@' (Pk) + PbO” + 3 PbO 
bzw. 


Nd,O, = 2 Nd@' (Pb) + PbO” + 3 PbO 


CeO, = (Pb) ++ PLO” + 2PbHO 


ergeben würden. 


®) K. Hauffe, Ergebn. exakt. Naturwiss. 25, 271 (1961). 
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In allen diskutierten Bruttoreaktionsgleichungen (Gln. (7), (9), (12) und 
(13)) tritt also durch den Einbau höherwertiger gitterfremder Kationen keine 
Änderung in der Konzentration der Elektronenstörstellen ein. Damit wäre 
demnach die Konstanz der elektrischen Leitfähigkeit des Blei(II)-oxyds bei 
Zusatz höherwertiger gitterfremder Kationen erklärbar. Auch die den Kurven 
Abb. 1 bis 3 noch zu entnehmenden geringen Schwankungen der elektrischen 
Leitfähigkeit sind erklärlich, wenn man annimmt, daß unter Umständen die 
Kompensation der Konzentrationsänderung der Elektronenstörstellen nicht 
in so idealer Weise erfolgt, wie es in den abgeleiteten Reaktionsgleichungen 
angenommen worden ist. 

Noch eine andere Möglichkeit zur Erklärung des eigentümlichen Verhaltens 
der elektrischen Leitfähigkeit des Blei(II)-oxyds bei Zusatz höherwertiger 
gitterfremder Kationen ist zu berücksichtigen. Es ist nämlich weiterhin 
denkbar, daß das Blei(II)-oxyd eine gewisse Teilleitfähigkeit durch Kationen 
oder Anionen zeigt, und daß bei angenommener Gültigkeit der Einbau- 
gleichungen (la) bis (3c) zusätzlich zu der Veränderung der Elektronenstör- 
stellenkonzentration durch die eingebauten höherwertigen Fremdionen eine 
gegensinnige Veränderung der Konzentration oder Beweglichkeit der Ionen- 
störstellen erfolgt, die dann wiederum bis zur Leitfähigkeitskonstanz führen 
könnte. 

Aufschluß über die im Blei(TI)-oxyd tatsächlich vorliegenden Fehlordnungs- 
erscheinungen können nur weitere Messungen bringen, z. B. Leitfähigkeits- 
messungen an PbO mit Zusätzen an niederwertigeren Kationen, Unter- 
suchung der Abhängigkeit der elektrischen Leitfähigkeit vom Sauerstoff- 
partialdruck. Überführungsmessungen der Kationen und Anionen sowie 
Thermokraftmessungen. 

Die Arbeiten werden fortgesetzt. 


Potsdam, Sanssouci, Pädagogische Hochschule, Institut für Experi- 
mentalphysik. 


van. Bei der Redaktion eingegangen am 31. Oktober 1957. 
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Die Bewegungsgleichungen der unitären Feldtheorie a 
in der niedrigsten Annäherung 
Inhaltsiibersicht 


Nach der Zusammenstellung der Feldgleichungen der vom Verfasser ent- 
wickelten unitären Feldtheorie werden diese annäherungsweise integriert. 
Die Bewegungsgleichungen werden in der niedrigsten Annäherung nach der 
bekannten Methode von Einstein und Infeld abgeleitet. Sie haben die 
richtige Form, denn sie enthalten nicht nur die Newtonsche, sondern auch 
die Coulombsche Kraft. 


Einleitung 


In den vorangehenden Arbeiten!) hat der Verfasser eine neue Form der — 
unitären Feldtheorie entwickelt, ausgehend aus einem heterogenen Hamilton- 
ian. Die physikalische Interpretation der einzelnen geometrischen Felder 
ist in dieser Theorie besonders einfach. Die Theorie enthält die Maxwellschen 
Gleichungen in ungeänderter Form. Die strenge statische kugelsymmetrische 
Lösung der Feldgleichungen führt zu den von der allgemeinen Relativitäts- 
theorie erwarteten Beziehungen. Die vorliegende Arbeit befaßt sich mit der 
Ableitung der Bewegungsgleichungen direkt aus den Feldgleichungen. 

Seit der Arbeit von Einstein und Grommer?) ist es bekannt, daß man 
zur Ableitung der Bewegungsgleichungen in der allgemeinen Relativitäts- 
theorie keine weitere Hypothese braucht, sondern daß diese aus den Feld- 
gleichungen folgen. Das Problem der Bewegung eines Probeteilchens im 
reinen Gravitationsfeld haben Infeld und Schild) gelöst. Diese Verfasser 
haben dabei nur von den Feldgleichungen für die leere Raum-Zeit Gebrauch 
gemacht und sind zu denselben Resultaten gekommen, die man mit Hilfe 
des Variationsprinzips der stationären Teilchenwirkung erhält®). Die Be- 
wegungsgleichungen eines Probeteilchens im unitären Feld wurden, soweit 
dem Verfasser bekannt, noch nicht direkt aus den Feldgleichungen abgeleitet. 
Deshalb hat der Verfasser die Bewegungsgleichungen eines Probeteilchens vor- 
läufig auch aus dem genannten Variationsprinzip in der Arbeit II berechnet 
und in III integriert, um festzustellen, inwieweit sie im Einklang mit der 
Relativitätstheorie stehen. 


1) J. Pachner, Ann. Physik 19, 353 (1957) 20, 368 (1957); 1, 112 (1958). Diese 
Arbeiten werden hier als bzw. I, II, III bezeichnet. Ebenso werden die Gleichungen aus 
diesen Arbeiten als bzw. (I...), (II...), (III...) zitiert. 

2) A. Einstein, J. Grommer, Sitz.-Ber. Berlin. Akad. 2 (1927). 

3) L.Infeld, A. Schild, Rev. mod. Physics 21, 408 (1949). 

4) Siehe z.B.: L. Landau, E.Lifschitz, Teorija polja, 2. Aufl., Moskva-Lenin- 
grad 1948, S. 280-283, 286. 
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Für die Bewegung von vergleichbaren Teilchen haben Einstein und 
Infeld5) eine Methode entwickelt, die später Infeld®) auch auf die ältere 
Form der Einsteinschen unitären Theorie erweitert hat. Dabei hat er fest- 
gestellt, daß die aus den Feldgleichungen berechnete Bewegung in der nied- 
rigsten Annäherung durch das elektromagnetische Feld unbeeinflußt wird. 
Zu demselben Ergebnis ist auch Callaway’) bei der Behandlung der neueren 
Einsteinschen unitären Theorie gekommen. 

Diese Resultate haben Bonnor*) veranlaßt, den Hamiltonian der Ein- 
steinschen unitären Theorie auf Grund einer ausführlichen Diskussion so 
abzuändern, daß in den Bewegungsgleichungen nicht nur die Newtonsche, 

sondern auch die Coulombsche Kraft erscheint. Dabei identifiziert er — 
ebenso wie Infeld®) und Callaway’) — die Komponenten g,, mit dem elek- 
trischen und die g,, mit dem magnetischen Feld (i, k = 1, 2, 3). 

Der Hamiltonian der vom Verfasser in den Arbeiten I und II entwickelten 
unitären Theorie unterscheidet sich von dem Bonnorschen nur durch ein 
Glied, das die kosmologische Konstante enthalt. Trotzdem kann man nicht 
die Bonnorschen Bewegungsgleichungen übernehmen, da die strenge statisch 
kugelsymmetrische Lösung der Feldgleichungen in der Arbeit I durch Ver- 
gleich mit dem entsprechenden Fall der allgemeinen Relativitätstheorie 
gezeigt hat, daß die Komponenten g,, umgekehrt mit dem magnetischen und 
die g,, mit dem elektrischen Feld zu identifizieren sind. Es ist daher unsere 
Aufgabe, die Bewegungsgleichungen unter dieser Voraussetzung zu berechnen. 


Die Feldgleichungen 
Zuerst stellen wir die Feldgleichungen der unitären Theorie zusammen. 
In der Arbeit II haben wir gefunden: 


Quv,e — — ev Qua 0, os (4) 
wobei 
Ruy = t+ + hes (5) 
A Gu v + x Gu B 9x Jap Ju Iu a (6) 
und Run, oe) bzw. Aus, oe) die zyklische Divergenz von Ruy bzw. H ur bedeuten. 
® 
Strich bezeichnet die partielle Ableitung. 
Die der Gl. (3) vollständig äquivalente Beziehung lautet 


Wir nehmen an, daß der Tensor g,, dem elektromagnetischen Feld (H, D), 
der Vektor J’, dem Maxwellschen Vektorpotential und der durch Rotation 


. 5) A. Einstein, L. Infeld, Canad. J. Math. 1, 209 (1949). 
8%) L.Infeld, Acta Phys. Polon. 10, 284 (1950). 


f 

9) J.Callaway, Physic. Rev. 92, 1567 (1953). an ee 
_ 8) W. B. Bonnor, Proc. Roy. Soc. (London) A 226, 366 (1954). t EK 
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nd | aus dem Vektorfeld J", entstehende Tensor Fy, 8” 
ere 


st- F,=T,u-Tur (2) 

ad - dem elektromagnetischen Feld (B, E) proportional sind. Das Feld F,,, wird 

rd. mit der Gleichung 

en = Ruy + Pak: 

in- 

so Angenäherte Lösung der Feldgleichungen 

he, Infeld®) löst das System von 64 algebraischen Gln. (4) durch sukzessive 

em, Approximation. Nach ihm ist die affine Konnexion J‘,, durch die Beziehung 

ad Mi» (9) 
pen 


ein | gegeben, das bekannte Christoffelsche Symbol ist und mit 


eh der Formel 


er- hee — — Gus My. — M,.) (10) 
bestimmt wird, wobei 


Luve = 3 (11) 


en. Der Doppelpunkt in Gl. (10) bezeichnet die Riemannsche kovariante Ab- 
leitung. Der symmetrische kontravariante Tensor h*¢ wird aus dem kova- 
rianten symmetrischen Tensor he, 


en. = Jeo (12a) 
nach der Formel ‘gas 
(1) log det hag + gt? noel (12b) 
(2) *e bras 
gebildet. 
(3) Auf Grund der strengen statischen kugelsymmetrischen Lösung der Feld- 


(4) gleichungen aus der Arbeit I setzen wir in der niedrigsten Annäherung voraus, 
GAD 


(sin. (1 I. 93 ~ — Hy, Jo3 m Aj, 
~ 


(5) (13) 
(6 wiows, web mi 9a ~D,, 924 ~D,, ~ D3. | 
- Den fundamentalen Tensor g,, entwickeln wir in die folgenden Potenz- 
; reihen nach einem sehr kleinen Parameter e: Pre 
= — Oman + & Imn + A + (14a) 
IP 
Ima = Ima + & Ima . (14e) 


> 
a) 

ere 

Ann. Physik. 7. Folge, Bd. 1 


The = Mis = 


354 Annalen der Physik. 7. Folge. Band 1. 1958 


Die lateinischen Indizes nehmen nur die Werte 1, 2, 3 ein. Die Ballen unter 
den einzelnen Zeichen geben die Ordnung der Annäherung an. Im Ausdruck e* 
bezeichnet n den Exponent, nicht einen kontravarianten Index. 

Die Entwicklungen (14d, e) unterscheiden sich von den in den Arbeiten 
von Infeld®), Callaway’) und Bonnor®). Der Grund für diese Änderung, 
die Infeld®) zuläßt, liegt in einer anderen Annahme über die physikalische 
Bedeutung der Komponenten g,,. Aus der Annahme (13) folgen unmittelbar 
die Entwicklungen (14d, e). ee 

Die Bewegungsgleichungen werden in der quasistatischen Annäherung 
bestimmt, d.h. unter der Voraussetzung, daß die Teilchen sich so langsam 
bewegen, daß durch die Differentiation nach der Zeit (nach x,) der Exponent 
des Parameters im betreffenden Glied um Eins erhöht wird. 


Nach den vorangehenden Gleichungen haben wir die M,„, bestimmt, 


die wir bei der weiteren Berechnung brauchen: _ ten I une 
4 ig’ 
Mi, = — Tins, (15a) 
ia = Ma = — Gia + ths (15b) 
dng Mix =—Mis=— inn + Tien, (15e) 
thd Miu = = Mir = = 9k4,4; (15d) 
4 or 
Mis =—2 (es. ar Ian). reins! (15f) 


Dadurch haben wir alles zur Lésung der Feldgleichungen vorbereitet. 
In der zweiten und dritten Ordnung der Annäherung werden diese Gleichungen 
in die des Gravitationsfeldes und in die des elektromagnetischen Feldes ge- 
spaltet. Mit der Berechnung der Gleichungen des Gravitationsfeldes werden 
wir uns nicht befassen, da sie mit den Gleichungen der allgemeinen Relativitäts- 
theorie identisch sind. Wir nehmen deshalb die Ergebnisse der Lösung aus 
der Arbeit von Einstein und Infeld°®) über. 

Die Gin. (3a) des elektromagnetischen Feldes nehmen in der zweiten und 
dritten Ordnung der Annäherung folgende Form ein: 
(16a) 


Mi, +t nt = — — = = 0. 
Wir suchen Lösung in welcher die Singularitäten 
im Feld elektrisch geladene Massenpunkte darstellen. Mit Rücksicht auf die 


Annahme (13) finden wir 


=®, (17a) 


®) Siehe Fußnote 5, S. 227-229. 
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wobei die skalare Funktion ® der Laplaceschen Differentialgleich ung 

A®=0 (18) 

genügt. Im Falle des Zweikörperproblems, auf das wir uns der Einfachheit 


wegen beschränken (obwohl die Lösung des Mehrkérperproblems keine 
Schwierigkeiten macht) lautet ihr partikulares Integral 


wo 
=) 


die räumliche Entfernung zwischen dem Punkt mit den Koordinaten %, und 


dem k-ten Teilchen mit den Koordinaten gm bezeichnet. In Gin. (19) und (20) 
geben die Indizes in der Klammer an, um welches Teilchen es sich handelt. 
Die Größen 2 sind die Integrationskonstanten. 


2 
Der mit Gl. (6) definierte Tensor H,, ist in der zweiten und dritten Ord- 
me mit den Beziehungen 


bestimmt. Dabei haben wir das die kosmologischen Konstante enthaltende 
Glied vernachlässigt. 

Die weiteren Gin. (2) des elektromagnetischen Feldes gehen nun in die 
folgende Form über 


Durch Einsetzen von Gin. (17a, b) stellen wir fest, daß die gefundene Lösung 
der Gin. (16a, b) auch den Gin. (22a, b) genügt. Paes 


Die Bewegungsgleichungen in der niedrigsten Annäherung 


Callaway’) hat gezeigt, daß nur die die symmetrischen Tensoren ent- 
haltenden Feldgleichungen die Bewegungsgleichungen der unitären Feld- 
theorie von Einstein bestimmen. Dasselbe gilt auch in der vorliegenden 
Form der unitären Theorie. Es folgt weiter aus den Arbeiten von Einstein 
und Infeld5)*), daß die Integrabilitätsbedingung der Gin. (1) ist, daß ge- 
wisse Integrale verschwinden, die über eine beliebige geschlossene zwei- 
dimensionale, das k-te Teilchen umhüllende Fläche 0 genommen werden. 
Und gerade diese Bedingungen sind die gesuchten Bewegungsgleichungen, 
die in der niedrigsten Annäherung lauten: 
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wobei 
— An, = Pas + Ome P 44> (24a) 


Pm Ar 
Pmt) — Ome (R (Ru Pa) (24b) 
Hier bezeichnet P,, den symmetrischen verjüngten Riemann-Christoffel- 


schen Krümmungstensor, der mit a nach Gl. (5) gebildet wird. 


Wenn das Integral in Gl. (23) von der Form der Fläche o) nicht ab- 
hangen soll, so miissen die Bedingungen 


ag Amr,r =0 dick (25c) 
Mit der Berechnung des Integrals 
1 
2 Ay My 
4a f 


werden wir uns nicht befassen, da wir ihren Wert in der Arbeit von Einstein 
und Infeld®) finden, woher wir ihn in die resultierende Bewegungsgleichung 
einsetzen. 

Setzen wir die mit Gln. (17a, b) gegebene Lösung der Feldgleichungen 
in die Formeln (15a—f) ein, so finden wir, daß 


— 4 — 


m (le din 8 n 8 
ug Ru Pu ‘ge Miss —M,, Mn = 3 @, D 
4 4 4 2 2 2 2 a 
und weiter nach Gl. (24b) <A 
den Amr — Ome a ait 


Der Ausdruck in der eckigen Klammer, den wir Y,,,„ bezeichnet haben, ist 
in den Indizes (r, n) schiefsymmetrisch : 


mnr* 
Durch Differentiation stellen wir fest, daß A;,,, die Bedingung (25b) erfüllt. 


Wie in der genannten Arbeit von Einstein und Infeld!®) ausführlich be- 


gründet, dann das 


ax iz {24 2 Amr n, do en? do = 0. ab (86) 

iehe hie 

10) Siehe Fußnote 5, S. 213—214. at 
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Der mit Gl. (6) definierte Tensor H,,, ist in der vierten Ordnung der An- 


näherung mit dem Ausdruck 
| ore 
-b) Har =2x (Ims 3 Ges (27) 
ce) gegeben. Setzen wir ihn in Gl. (24c) ein, so erhalten wir mit Rücksicht auf 
Gl. (17a) 
>] - 2 Amr = — 4% (OP, ,). (24e’) 
= 2’ 2 2’ 2 
b- Durch Differentiation stellen wir fest, daß auch der Tensor Ay, der Be- 


dingung (25c) genügt, da die Funktion ® das Integral der Laplaceschen 


2 
a) Differentialgleichung (18) ist. Bei der Berechnung des letzten Flächen- 
integrals benützen wir die Beziehungen aus der wiederholt genannten Arbeit 


b) von Einstein und Infeld!!). Wir finden ; 
2 Aims do = — 87) 
Hier bezeichnet o die Entfernung eines Teilchens vom anderen: 
in Nun erhalten wir die gesuchten Bewegungsgleichungen, wenn wir in 


ng Gl. (23) den Wert aus der Arbeit®) und die Gln. (26) und (27) einsetzen. Nach 
einer kleinen Umformung, bei der wir in der niedrigsten Annäherung (im Gauß- 
en schen Maßsystem) schreiben 


finden wir folgende Bewegungsgleichungen für das erste von én zwei Teilchen, 


b’) die wir uns beschränkt haben, 
= — mn” (2) (1) ge ga) 93 


In der dreidimensionalen nobilitas geht diese Gleichung in die be- 
kannte Form 


2 


(23b) 


über. 
Die Feldgleichungen bestimmen daher in der vorliegenden unitären Theorie 
6) die Bewegung des Teilchens ohne irgendeine weitere Hypothese. Die Gin. 
(23a, b) die die niedrigste Annäherung darstellen, enthalten auf der rechten 
Seite nicht nur die Newtonsche (das erste Glied) sondern auch die Cou- 


4 11) Siehe Fußnote 5, S. 238. 


| 
| 

2: 
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lombsche Kraft (das zweite Glied). Die letztere erscheint hier als Folge der 
Einführung des Tensors H,, in die Feldgleichungen. In der Einsteinschen 
unitären Theorie (die wir aus unseren Gleichungen erhalten, wenn wir H,, 
gleich Null setzen) erscheint die Coulombsche Kraft nicht®)?). 

Die Bewegungsgleichungen (23a, b) hat schon Bonnor®) erhalten, obwohl 
er aus einer anderen Annahme über die physikalische Bedeutung des Feldes 
Gus und aus anderen Entwicklungen fiir Imn und Ing ausgegangen ist. 

Vergleichen wir nun Gin. (30a, b) mit (LE 69),so finden wir, daß die dimensions- 
lose Konstante b den Wert 

b=+1 «S00, F20 (31) 
2 


haben sollte. Auf der anderen Seite hat sich die Formel (III.24) aus der Lé- 
sung der aus dem Variationsprinzip der stationären Teilchenwirkung abge- 
leiteten Bewegungsgleichungen ergeben, nach der N 


b=-}4, x>0. (III.24) 

Die Ergebnisse der Untersuchung in der vorliegenden Arbeit und in der Achs 
Arbeit III haben daher zwar die erwartete Form, aber die Konstante 6 soll War 
nach jeder von diesen Arbeiten einen anderen numerischen Wert annehmen. wird 


Eine vollständige Klärung dieses Widerspruchs kann erst die Ableitung der Glim 
Bewegungsgleichungen eines Probeteilchens aus den Feldgleichungen bringen. 


rr Prag, Physikalisches Institut der Technischen Hochschule. 
7 der Redaktion i angen am 6. Dezember 1957. a 
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er Der Potential. und Feldstärkeverlauf Pr 
Im Kathodenfallgebiet von Glimmentladungen’) 
Von Ginter Brederlow 
Mit 16 Abbildungen 


pa 


i Inhaltsiibersicht 


Mit einer Gliihkathode wird im Kathodenfallgebiet normaler und anomaler 
He-, Ne-, A-, H, und N,-Entladungen die Potentialverteilung bestimmt und 
durch Differentiation der Potentialkurven der Feldstärkeverlauf in der 
Achse der Entladung ermittelt. Die MeBergebnisse werden mit der von 
Warren angegebenen Beweglichkeitstheorie verglichen, wobei festgestellt 
wird, daß diese Theorie in der ersten Hälfte des Kathodenfallgebietes zwischen 
Glimmlicht und Kathode die Verhältnisse richtig beschreibt. In unmittel- 
barer Kathodennähe treten jedoch Abweichungen auf, für die versucht wird, 
eine Erklärung zu finden. 


Zur Ermittlung der Feld- bzw. Potentialverteilung im Kathodendunkel- 
raum von Glimmentladungen wurden bisher verschiedene Meßmethoden heran- 
gezogen. Ausgehend von der statischen Sonde benutzte man neben der von 
Brown und Thomson?) angegebenen Schattenmethode den Elektronen- 
strahl®), den Starkeffekt*) und die Glühsonde®). Leider weichen die bisher 
veröffentlichten Meßergebnisse teilweise stark voneinander ab. Ziel der vor- 
liegenden Arbeit ist es daher, mittels einer wenig benutzten aber zuverlässig 
erscheinenden Meßmethode den Potential- und Feldstärkeverlauf im Fall- 
gebiet normaler und anomaler Glimmentladungen zu bestimmen. Nach 
Angaben von W.H. Ernst) wurde mit einer Glühsonde die Potentialver- 
teilung in der Achse von He-, Ne-, A-, H, und N,-Entladungen ermittelt und 
durch Differentiation der Feldstärkeverlauf gewonnen. 


Eine Glühsonde stellt einen Elektronen emittierenden Körper dar. Bringt 
man diese in den Fallraum einer Glimmentladung und verändert ihr Potential 
gegenüber der Kathode, dann setzt Elektronenemission ein, wenn die Sonde 
negativ gegenüber dem Raumpotential geladen ist. Vernachlässigt man 
die thermische Energie der ausgesandten Elektronen, dann darf das Potential, 

1) Auszug aus der Dissertation Greifswald 1957. 

2) W.L. Brown u. E. E. Thomson, Philos. Mag. 8, 919 (1929). 

8) F. W. Aston, Proc. Roy. Soc. (London) A 84, 526 (1911); T. Harris, Philos. 
Mag. 80, 182 (1915); A. E. M. Geddes, Proc. Roy. Soc. (Edinburgh) 46, 136 (1926); 
R. P. Stein, Physic. Rev. 89, 134 (1953); P. E. Little u. A. v. Engel, Proc. Roy. Soc. 
(London) A 224, 209 (1954); R. Warren, Physic. Rev. 98, 1650 (1955). 

4) E. Brose, Ann. Physik 58, 731 (1919); W. Steubing, Ann. Physik 10, 297 (1931). 

5) W.H. Ernst, Helv. phys. Acta 9, 381 (1935). 
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bei dem Elektronenemission einsetzt, dem Raumpotential ne w wien 
Messungen mit dieser Sonde kénnen ohne genaue Kenntnis des Entladungs- 
mechanismus durchgeführt werden, wobei die Möglichkeit besteht, die Zuver- 
lässigkeit der Meßmethode nachzuprüfen. 


Apparatur und Meßmethode 


Abb. 1 zeigt das Entladungsrohr. Um in verschiedenen genauestens meß- 
baren Abstäuden an der Kathode Sondenmessungen durchzuführen, wurde 
die Kathode verschiebbar angebracht, während die Glühsonde feststehend war. 
Die Sonde befindet sich hier in der Achse der Entladung. Es wurden auch 
Messungen mit einer in radialer Richtung verschiebbaren Sonde durchgeführt. 


“er 


achires 


Abb.1. Skizze des Entladungsrohres 
Wurde diese Sondenform verwandt, dann fielen am Entladungsrohr die beiden 
senkrecht verlaufenden Ansätze fort und statt dessen wurde das Rohr mit 
einem 10 cm langen waagerechten Ansatz versehen, der die Sondenführung 
aufnahm. Der Kathodendurchmesser wurde hier auf 52,8 mm erhöht. Die 
Kathode war dann volumenfüllend. Als Kathodenmaterial konnte nur Wolf- 
ram verwendet werden und zwar deshalb, weil von dem stark geheizten Sonden- 
draht in geringen Mengen Wolfram verdampfte und sich auf der Kathode 
niederschlug. An diesen bedampften Stellen setzte bei Verwendung anderer 
Materialien die Entladung nicht an. 

Die Kathode wurde mittels eines Elektromagneten durch Drehen einer 
Schraubspindel von 0,5 mm Ganghöhe, die zugleich auch als Meßspindel 
diente, bewegt. Es bereitete zunächst Schwierigkeiten, eine einwandfrei funk- 
tionierende Anordnung zu schaffen. Störungsfrei arbeitete die Verstellvor- 
richtung für die Kathode erst dann, als eine im Hochvakuum ausgeglühte 
Molybdänspindel, die mit Kolloidgraphit bestrichen worden war, und eine 
Nickelmutter verwendet wurden. 

Der genaue Aufbau des Bewegungsmechanismusses und der Meßanordnung 
für die Abstandsbestimmung ist aus der Skizze Abb. 1 ersichtlich. Mit der 
Kathode konnte ein Bereich von 45 mm Länge überstrichen werden, wobei die 
jeweilige Lage mit einer Genauigkeit von + 0,01 mm bestimmt werden konnte. 
Die Zuverlässigkeit der Meßvorrichtung wurde jeweils vor dem Einbau in das 
Entladungsrohr mittels einer Zeiß-Meßuhr geprüft. Wie man aus Abb. 1 
ersehen kann, erfolgt die Stromzuführung durch einen Nickel-Eisen-Stab, 
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en. der zugleich auch die Aufgabe hatte, das ganze System zu halten, über die 
gs- Halterungsscheibe aus Kupfer und zwei an ihr befestigten Schleifkontakte 
er- auf die Spindel, die die Kathode tragt. An diesen Kupferplatten waren Glas- 
scheiben befestigt, deren Größe nahezu dem Rohrquerschnitt entsprachen. 
Diese Scheiben hatten die Aufgabe, das ganze System zu zentrieren und an der 
Kathodenseite den Entladungsraum möglichst gut von dem den Bewegungs- 
mechanismus der Kathode enthaltenden Teil des Entladungsrohres abzu- 


x trennen, um bei hoch anomalen Entladungen einen Ansatz der Entladung an 
den unter Spannung stehenden Teilen der Verstellvorrichtung zu verhindern. 

2 Es war notwendig, zur Abschirmung zwei Glasplatten in der in der Abb. 1 

’ gezeigten Anordnung anzubringen. Die der Anode abgewandte Seite der 

” Kathode und ein Teil der Spindel wurden durch eine Abschirmung aus Raso- 
thermglas verdeckt. Diese Anordnung erwies sich auch in hochanomalen 
Entladungen als sehr zweckmäßig. 

Die Messungen wurden mit zwei Sondenformen durchgeführt. Zunächst 
sei die verschiebbare Sonde beschrieben. Ihr Aufbau und ihre Abmessungen 
sind aus der folgenden Skizze — Abb. 2 — ersichtlich. 

Glaskapillaren 
A 

7 
Eisenkerne 
0 1 in 
Abb. 2. Skizze der beweglichen Sonde 

Die Sonde konnte durch einen Elektromagneten in jede gewünschte Lage 
en gebracht werden, wobei diese Lageveränderung durch Eisenkerne. die an der 
it Sonde befestigt waren, ermöglicht wurde. Eingeschmolzene Führungs- 
1g stützen verhinderten ein Drehen des Sondendrahtes aus der durch die Ka- 
ie thodenoberfläche bestimmten Ebene. Der Glühdraht bestand aus einem 
if - Wolframdraht von 0,027 mm Durchmesser. Geheizt wurde mit einem Strom 
n- von 0,4 A, wobei längs desselben ein Spannungsabfall von 2,8 V auftrat. 
le Die Abmessungen der Sonde wurden möglichst klein gehalten. Daher waren 


er die Störungen der Entladung, wie später gezeigt wird, sehr gering. Es war 
allerdings nicht möglich, den Durchmesser der Glaskapillaren, der an der 


er Spitze 0,25 mm betrug, noch kleiner zu wählen, da die Stabilität der Sonde 
el dann nicht ausreichend war. 

k- Da der Sondendraht nicht unmittelbar bis auf die Kathodenoberfläche 
r- gebracht werden konnte, war es erforderlich, den Nullpunkt der Abstands- 
te skala durch Messen von außen festzulegen. Die Abstandsbestimmung wurde 
1€ bei geheizter Sonde vorgenommen, und zwar wurde mit einem mittels einer 


Mikrometerschraube verschiebbaren Mikroskop durch eine geschwärzte 
gz Scheibe der Abstand des gliihenden Drahtes von seinem durch die polierte 
er Kathodenoberfläche erzeugten Spiegelbild gemesssen. Aus diesem Wert wurde 
ie unter Beriicksichtigung des Betrachtungswinkels der Abstand Sonde—Ka- 
e. thode bestimmt. Der mögliche Fehler dürfte hier + 0,03 mm nicht über- 
1S schreiten. Ein Einfluß elektrostatischer Kräfte machte sich nicht störend 
1 durch etwaige Abstandsänderungen bemerkbar, wie Untersuchungen mit er- 
b, heblich verschiedenen Sondenpotentialen ergaben. Auch ein Durchbiegen 
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der Sonde unter dem Aufprall positiver Ionen in der Entladung war nicht feder 


feststellbar. griffe 

Mit dieser Sonde war es möglich, den gesamten Entladungsraum vor der den ¢ 
Kathode von einer Rohrwand bis zur anderen zu erfassen. Die Messungen wider 
wurden in der Hauptsache jedoch nur in der dem Ansatz abgewandten Seite zwise 
der Entladung durchgeführt, da in Wandnähe am Ansatz des Führungs- Präzi 
rohres Abweichungen der Aquipotentialflachen von der sonst auftretenden wird 
Form möglich waren. Trotz 


Als Vorbild für die feststehenden Sonden diente die von W.H. Ernst notw 
angegebene Sonde. Jedoch wurde diese nicht in allen Einzelheiten getreu die ir 
nachgebaut. Die Abmessungen und der Aufbau sind aus Abb. 1 ersichtlich. nung 
Die elektrischen Daten und der Durchmesser der Glaskapillaren unterschieden nach: 
sich nicht von denen der beweglichen Sonde. Um die beim Heizen der Sonde wider 
durch Ausdehnen des Wolframdrahtes möglichen Lageveränderungen zu Lage. 
verhindern, wurde diese Sonde durch eine Stahlfeder von 0,1 mm Stärke ge- Span 
spannt. Trotz dieser mechanischen Beanspruchung konnte der Wolfram- der n 
draht auf helle Gelbglut erhitzt werden, ohne zu zerreiBen. Der Nullpunkt ermit 
der Abstandsskala wurde nur bei geheizter Sonde bestimmt, und zwar wurde genaı 


mit der Stromstärke geheizt, bei der in der Entladung gemessen wurde. differ 

Der Fülldruck im Entladungsrohr wurde vor jeder Messung mit einem erfaß 
McLeod bestimmt. Die Meßgenauigkeit lag bei + 0,01 Torr. Jedes Ent- D 
ladungsrohr wurde vor Inbetriebnahme und nach jedem Sondenwechsel in Meßl 


einem elektrischen Ofen ausgeheizt. Anschließend wurden beide Elektroden abfal 
ausgewirbelt. Alle Metallteile des Bewegungsmechanismusses für die Kathode diese: 


wurden einzeln vor dem Zusammenstellen im Hochvakuum geglüht. Scha! 
Der elektrische Aufbau der Meßanordnung, Meßmethode nd 

Die Schaltung des Entladungsstromkreises mit der damit verbundenen Stell 
Sonde ist der Abb. 3 zu entnehmen. inder 
Da mit der Sonde zur Feldstärkemessung Spannungsdifferenzen bis auf Entls 
weniger als 0,1 V genau erfaßt werden mußten, konnte mit einer einfachen treffe 


Spannungsteiler - Schaltung paral 


nicht gearbeitet werden, da dem 

höheren Spannungen die hang 

Ke erforderliche Meßgenauigkeit besti: 
& wore 98 nicht zu /erreichen war. Es den 

Mon wurde deshalb auf dem der war. 

 Spannungsquelle parallel ge- | Teile 

schalteten Widerstand zwei Bei « 

in festem Abstand verschieb- bis E 

Ian bare Schleifkontakte so an- abfal 

ZIKR gebracht, daß bei den wäh- emiss 

10k2 rend des Betriebes herr- poter 

2? schenden Bedingungen am trat, 

@ Amperemeter Voltmeter © Galvonometer Widerstand zwischen beiden } 

von Abb. 3. Schaltskizze Kontakten ein Spannungs- Char 

ie abfall von ungefähr 70 V auf- teris 

trat. Über einen zweiten windungsreichen Widerstand von 3k2 wurden beide posit 
Schleifkontakte verbunden. Auf diesem Widerstand konnte über eine Kontakt- bei 
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feder durch eine Drehvorrichtung an beliebigen Stellen die Spannung abge- 
griffen werden. Die durch diesen Nebenwiderstand gegebene Spannung gegen 
den dem Kathodenpotential am nächsten liegenden Kontakt auf dem Haupt- 
widerstand wurde mit einem Präzisionsvoltmeter gemessen. Die Spannung 
zwischen diesem Kontakt und der Kathode wurde ebenfalls mit einem 
Präzisionsvoltmeter bestimmt. Die Sondenspannung gegen die Kathode 
wird durch die Summe der beiden gemessenen Spannungswerte angegeben. 
Trotz des veränderlichen zusätzlich durchfließenden Sondenstromes war die 
notwendige Spannungskonstanz gewährleistet, da die Sondenstromänderung, 
die im allgemeinen nur einige uA betrug, in ungünstigsten Fällen eine Span- 
nungsänderung kleiner als 0,1 V ergab. Während einer Meßreihe wurde zu- 
nächst die Sondenspannung mittels des verschiebbaren Abgriffes am Neben- 
widerstand verändert. Der Schieber mit den Doppelkontakten behielt seine 
Lage. Konnte der neue Meßpunkt so nicht mehr erfaßt werden, dann wurde die 
Spannung durch den Doppelschieber verändert. Jedoch wurde nun nicht 
der neue Meßwert bestimmt, sondern es wurde der letzte Wert noch einmal 
ermittelt, da es mit dem Doppelschieber nicht möglich war, eine genügend 
genaue Spannungsänderung zu erreichen. Nur so konnten die Spannungs- 
differenzen zwischen allen Meßpunkten mit der erforderlichen Genauigkeit 
erfaßt werden. 

Der Sondenstrom wurde mit einem Mikroamperemeter mit verschiedenen 
Meßbereichen gemessen. Durch das Heizen der Sonde trat ein Spannungs- 
abfall längs derselben auf, der die Potentialmessungen verfilschte. Um 
diesen Fehler möglichst herabzusetzen, wurde die von W. H. Ernst angegebene 
Schaltung, die der einer Wheatstonschen Brücke gleicht, angewandt. In 
jeder Stromverzweigung gibt es in den einzelnen Zweigen Punkte, zwischen 
denen keine Potentialdifferenz auftritt. Diese sich im Potential entsprechenden 
Stellen galt es zu ermitteln. Dies war nur auf indirektem Wege möglich, 
indem auf einem möglichst steilen Stück der Charakteristik bei brennender 
Entladung und einem Sondenpotential, das größer war, als das an dem be- 
treffenden Ort herrschende Raumpotential, der Schleifkontakt des der Sonde 
parallel geschalteten Widerstandes so lange verschoben wurde, bis der mit 
dem Galvanometer G gemessene Sondenstrom von der Heizleistung unab- 
hängig war. Das Potential der Sondenmitte wurde so mit dieser Schaltung 
bestimmt. Der Heizstrom der Sonde war stets so zu bemessen, daß nur in 
den allernächsten Bereichen der Sondenmitte Elektronenemission möglich 
war. Es darf angenommen werden, daß längs des Elektronen emittierenden 
Teiles des Wolframdrahtes ein Spannungsabfall von höchstens 2 V auftrat. 
Bei den Potentialmessungen wurde das Sondenpotential so lange verringert, 
bis Elektronen emittiert wurden. Wegen des vorher beschriebenen Spannungs- 
abfalles am Sondendraht war es möglich, daß das Potential, beidem Elektronen- 
emission einsetzte, im ungünstigsten Falle um 1 V höher war als das Raum- 
potential. Da dieser Fehler bei allen Messungen als konstanter Faktor hinzu- 
trat, verfälschte er die Feldstärkewerte nicht. 

Ein Vergleich der mit geheizter und ungeheizter Sonde aufgenommenen 
Charakteristiken liefert das Raumpotential. Abb. 4 zeigt solche Charak- 
teristiken, die im Kathodendunkelraum aufgenommen wurden. Ist die Sonde 
positiv gegenüber dem Raumpotential geladen, dann unterscheiden sich die 
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Wird dagegen das Raumpotential erreicht, dann werden Elektronen aus- 
gesandt; der Sondenstrom steigt sprunghaft an. Die Charakteristik läßt er- 
kennen, daß das Heizen der Sonde die Entladung nicht störte, da bei Sonden- 
spannungen, die höher sind als das Raumpotential, die Charakteristik der ge- 
heizten und der kalten Sonde zusammenfallen. In manchen Fällen, wenn der 
Heizdraht zu kurz war, traten Störungen 
auf, die auf eine positive Aufladung der 


SPÄTER Glasisolation der Sonde zurückzuführen 
MEE at sind. Diese Stérungen wurden jedoch bald 
erkannt und durch Verlängern des Sonden- 

m 30 wi“ drahtes beseitigt. Mit dieser Methode, die 


von W.H.Ernst als „Methode des Ver- 
zweigungspunktes‘‘ bezeichnet wurde, waren 


0: jedoch nur Potentialbestimmungen im 
Kathodendunkelraum und in manchen 
' Fallen in angrenzenden kathodennahen 
Abb. 4. Sondencharakteristiken im Teilen des negativen Glimmlichtes möglich. 


Kathodendunkelraum Im negativen Glimmlicht — zumindest in 
kathodenferneren Bereichen — und im 
Faradayschen Dunkelraum war diese Methode nicht anwendbar. Hier 
wurde wegen der vorherrschenden hohen Konzentration an Elektronen der 
Sondenstrom bei Sondenspannung in unmittelbarer Nähe des Raumpoten- 
tials hauptsächlich durch Elektronen bestimmt. Die Sondencharakteristik 


verlief sehr steil. Sie zeigt den in Abb. 5 an- 


J | gegebenen Verlauf. Die ungeheizte statische 

u A Sonde lädt sich schwach negativ gegenüber 
7 dem Plasmapotential auf. Wird jetzt die Sonde 

50 geheizt, dann wird sich das Aufladepotential 
zu höheren Spannungswerten verschieben, da 

ie ein Teil des einfließenden Elektronenstromes 
a — halte Sonde durch die ausgesandten Glühelektronen kompen- 
» ~~-gpneineSone| siert wird. Wird der gesamte Elektronen- 
anlaufstrom durch genügend starke Heizung 

5 ausgeglichen, dann gibt das Potential, bei dem 

: der Sondenstrom bei geheizter Sonde Null ist, 

das Raumpotential an. Bei Sondenpotentialen, 

2 die höher sind als das Raumpotential, erfolgt 

’ keine Elektronenemission. Jedoch fallen hier 

die bei geheizter und kalter Sonde aufgenom- 


u menen Charakteristiken nicht zusammen. Die 
Bi ı 320 340 \W Ursache dieser Abweichungen konnte nicht 


mM 


Abb. 5. Sondencharakteristik 
im negativen Glimmlicht 


festgestellt werden. Diese Meßmethode wird 
von W.H. Ernst als ‚„Nullstrommethode‘‘ be- 
zeichnet. Die thermische Belastbarkeit des 
Sondendrahtes wird nicht immer erlauben, 
die Temperatur desselben so zu erhöhren, daß die erforderliche hohe Elek- 
tronenemission gewährleistet ist. In der Praxis wird man daher mit Null- 


strommethode oft einen Potentialwert ermitteln, der kleiner als das Raum-: 
re ist. Diese ng wird aber um so Br sein, je steiler der 
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Verlauf der Charakteristik und je héher die Temperatur der Sonde ist. Durch 
eine mögliche Ortsabhängigkeit der Steilheit der Sondencharakteristik kann 
ein Fehler in die MeBergebnisse hineingetragen werden. Messungen in ver- 
schiedenen Bereichen des Glimmlichtes und des Faradayschen Dunkel- 
raumes in Stickstoffentladungen ergaben, daß dieser Fehler wegen seiner 
Geringfügigkeit nicht berücksichtigt zu werden braucht. 

In einigen Fällen konnte das Potential am kathodenseitigen Rand des 
negativen Glimmlichtes sowohl mit der ‚Methode des Verzweigungspunktes“ 
als auch mit der „Nullstrommethode‘‘ gemessen werden. Es wurde die ,,Null- 
strommethode‘“ angewandt, obgleich die mit ihr erhaltenen Meßergebnisse 
kleiner als die wirklichen Werte des Raumpotentials waren. In ungünstigen 
Fällen waren diese Abweichungen jedoch nicht größer als 2 V. 


Meßgenauigkeit, Reproduzierbarkeit 


Die Meßgenauigkeit dieses Verfahrens wird bestimmt durch die Störungen 
der Entladung durch die Sonde, durch die Zuverlässigkeit mit der das an 
dem betreffenden Ort herrschende Potential ermittelt wurde und durch die 
Genauigkeit, mit der der Abstand Sonde—Kathode gemessen wurde. 

Die Störungen der Entladung durch die Sonde waren hier leicht feststellbar, 
da die verschiebbare Sonde vollkommen aus dem Entladungsraum entfernt 
werden konnte. Es zeigte sich, daß sich beim Einführen der Sonde die 
Brennspannung erhöhte. Durch Messen des Kathodenfalles an zwei zur Rohr- 
achse symmetrischen Punkten bei gleichzeitiger Kontrolle der Brennspannung 
wurde festgestellt, daß die Erhöhung der Brennspannung durch Änderungen 
der Größe des Kathodenfalles bewirkt wurde. In He-Entladungen, wo die 
Störungen am größten waren und die Brennspannung sich in ungünstigsten 
Fällen beim Einführen der Sonde um 6%, vergrößerte, wurde an zwei zur Rohr- 
achse symmetrischen Punkten der Potentialverlauf mit der verschiebbaren 
Sonde aufgenommen. Die Meßergebnisse zeigten, daß die Störungen durch die 
Sonde im Fallraum praktisch vernachlässigt werden können. Unmittelbar 
vor dem Glimmlicht und in diesem selbst bewirkte die Sonde eine Potential- 
erhöhung. Den Potentialkurven konnte entnommen werden, daß die Feld- 
stärkewerte im Glimmlicht nicht merklich verfälscht wurden, da die Potential- 
kurven hier nahezu parallel verliefen. Im wesentlichen bewirkte das Ein- 
bringen der Sonde eine Verlängerung des Kathodendunkelraumes und hatte 
daher keinen die Verhältnisse verzerrenden Einfluß auf den Verlauf der 
Feldstärkekurven. Bei der feststehenden Sonde dürften die Störungen der 
Entladung denen der beweglichen entsprechen. Als Beweis für die Zuver- 
lässigkeit dieser Meßmethode darf man die Tatsache werten, daß die bis zur 
Kathode verlängerten Potentialkurven stets durch den Nullpunkt gehen. 
In normalen bzw. schwach anormalen He-Entladungen bei geringen Drucken, 
wo die Astonsche Dunkelraumlänge größer als 0,25mm war, konnte die 
Zuverlässigkeit der Meßmethode auch daran geprüft werden, ob am Rande 
des Astonschen Dunkelraumes die Anregungsspannungen der im folgenden 
Entladungsteil ausgesandten Linien gemessen werden. Solche Messungen 
wurden durchgeführt. Die Abweichungen der Meßergebnisse von den An- 
regungsspannungen waren nicht größer als 2 V. Diese Abweichungen können 
dadurch bedingt sein, daß es nicht möglich war, die Sonde wegen der in diesem 
Bereich hohen Feldstärke mit der erforderlichen Genauigkeit an den Rand 
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des Astonschen Dunkelraumes zu bringen. Ob die mit der Nullstrommethode 
ermittelten Werte das tatsächliche Raumpotential wiedergaben, konnte 
nur mit der von Langmuir und Mott-Smith®) angegebenen Sonden- 
methode geprüft werden. Es zeigte sich, daß in manchen Fällen bei nicht ge- 
nügend starker Elektronenemission der Sonde das gemessene Potential um 
höchstens 2 V kleiner war als das tatsächliche. 
Um möglichst gut reproduzierbare Verhältnisse zu erzielen, waren saubere 
Entladungsbedingungen erforderlich. Die Durchführung der vorher be- 
schriebenen Reinigungsprozesse des Entladungsrohres und der Elektroden 
gewährleistete bei Verwendung sehr reiner Gase saubere Verhältnisse. Das 
zeigt die Tatsache, daß in He- und Ne-Entladungen stets die Aureole auftrat. 
Beim Brennen der Entladung in Edelgasen machte sich eine Erscheinung 
bemerkbar, die bereits von R. Seeliger und W. Bartholomeyezyk’) und 
F. M. Penning und J. H. A. Moubis?) als ‚„Inseleffekt‘‘ beschrieben wurde. 
Der Entladungsansatz zog sich hier nach einer von den jeweiligen Verhältnissen 
abhängigen Brenndauer auf einen Teil der Kathode zusammen, wobei die 
Entladungsstromstärke sich stark erhöhte und die Brennspannung sich ver- 
ringerte. Bei Steigerung der Entladungsstromstärke vergrößerte sich diese 
„Insel‘‘ und bedeckte nach genügender Erhöhung des Stromes nach einiger 
Zeit die gesamte Kathodenfläche. Nach Untersuchungen von F. M. Penning 
und J. H. Moubis wird diese Erscheinung dadurch hervorgerufen, daß eine 
monomolekulare Oxydschicht auf der Kathode die Austrittsarbeit für Elek- 
tronen erhöht und daher zunächst den Entladungsansatz auf die Teile der 
Kathode beschränkt, die schon frei von der Oxydschicht sind. 
In Wasserstoff und Stickstoff konnte dieser Effekt nicht beobachtet werden. 
Da bei Entladungen in diesen Gasen die Brennspannung sich allmählich er- 
höhte und nach einer gewissen Zeit einen konstanten, aber im Vergleich zu 
den entsprechenden Edelgasentladungen einen sehr hohen Wert erreichte, 
wurde angenommen, daß beide Gase mit der Wolframkathode reagierten und 
wie Sauerstoff eine monomolekulare Oberflächenschicht bildeten. Für diese 
Annahme spricht die Tatsache, daß die Kathodenzerstäubung in Wasser- 
stoff- und Stickstoffentladungen erheblich geringer war als in Edelgasent- 
ladungen, die Kathode also durch eine entsprechende Oberflächenschicht ge- 
schützt wurde. 
Bei Beachtung dieser Tatsachen war es möglich, gut reproduzierbare 
Entladungsbedingungen zu schaffen. Kontrollen des Kathodenfalles mit neuen 
Gasfüllungen und manchmal auch mit verschiedenen Entladungsrohren er- 
gaben in keinem Falle Abweichungen, die größer als 5% des vorher gemessenen 
Wertes waren. Kontrollen des Potentialverlaufes im Fallraum ergaben ähn- 
liche Ergebnisse. 


Potentialverteilung im Fallgebiet von He-, Ne-, A-, Ho- und No-Entladungen 


Mit der beweglichen Sonde wurde im Fallgebiet anomaler und normaler 
bzw. schwach anomaler Entladungen in He, Ne, A, H, und N, die Potential- 
verteilung ermittelt. Die Ergebnisse dieser Messungen sind im Ausschnitt 


s) J. Langmuir u. H. Mott-Smith, Gen. Elektr. Rev. 27, 449, 538, 616, 762, 810 
(1924). 

7) R. Seeliger, W. Bartholomeyczyk, Ann. Physik 1, 241 (1947). 

8) F.M. Penning, J. H. A. Moubis, Phil. Res. Rep. 1, 119 (1945/46). 
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in Abb. 6 wiedergegeben. Es ist hier nur der Potentialverlauf in Ne-Ent- 
ladungen gezeigt, da die Potentialverteilung in diesem Gas charakteristisch 
für alle untersuchten Gase in entsprechenden Stromstärke- und Druck- 
bereichen ist. In den Abb. 6 sind die Linien gleichen Potentials angegeben. 
Die Zahlen an der Rohrwand bezeichnen den Abstand von der Kathode, 
gemessen in mm, während der Rohrdurchmesser in cm unterteilt ist. Das @ 
an der Rohrwand gibt die Grenze des Glimmlichtes an. Die Zahlenangaben 
zwischen den einzelnen Äquipotentiallinien stellen die Spannungsdifferenzen 
zwischen denselben dar. Die Ka- 


thode war hier volumenfüllend. 
Die Kurven zeigen, daß das 2 
Potential in einer Ebene parallel ur —‘—is 34 
Kathode bei den höheren und mitt- = = | 
leren Drucken in Wandrichtung zu- 
nächst abnimmt und nach Durch- 6 
laufen eines Minimums in unmittel- é— 10 
barer Wandnähe dann wieder steil 57 
ansteigt. Die Tiefe dieses Minimums 
hing bei allen untersuchten Gasen \e/\ 
vom Fülldruck ab, und zwar nahm 
die Größe dieser Potentialeinsenkung J 30 
mit steigendem Druck zu. 
Die jeweiligen Lagen dieser 24 
Potentialminima geben die Grenzen 30 50 
des Bereichs an, in dem der Strom- | i 1 
transport erfolgte. Das zeigt die 
Tatsache, daß der Ansatz der Ent- re rn. 
ladung auf der Kathode sich nur bis -2 Ea) 
zu den durch sie angegebenen Ab- Ne anomal Ne schwach anomal 
ständen von der Rohrwand er- p=106 Tor ng p =101 Torr 
streckte, was dadurch bewiesen i= 60mA ve i =025 mA 
wird, daß nach Abschalten der Ent- Abb.6. Potentialverlauf in einer Ne-Ent- 
ladung ein Ring am Kathodenrand ladung. p = 1 Torr 


erkennbar war, dessen Breite dem 

nicht von der Entladung bedeckten Teil der Kathode entsprach. Bemerkens- 
wert ist die Tatsache, daß die Entfernung dieser Potentialeinsenkung von der 
Wand nur vom Fülldruck abhängig ist und nicht von der Entladungsstrom- 
stärke. Sie nahm mit abnehmendem Druck zu. In Wasserstoff und Helium, 
also bei den leichteren Gasen, war der Wandabstand bei entsprechenden 
Druckverhältnissen größer als bei den übrigen Gasen. Bei niedrigen Drucken 
verschwand besonders bei normalen Entladungen dieses Potentialminimum 
vollständig. 

Aus dem Abfall des Potentials in Richtung auf die Rohrwand folgt, daß 
die Stromdichte mit wachsendem Abstand vom Kathodenmittelpunkt ab- 
nahm. In Stickstoffentladungen durchlief die Stromdichteverteilung bei 
höheren Drucken ein Maximum, das in der Mitte zwischen der Wand und dem 
Kathodenzentrum lag. 

Der Abb. 6 kann man entnehmen, daß im Kathodendunkelraum das 
Potential in unmittelbarer Wandnähe stets steil anstieg, und die Wand Mr 
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dessen stark positiv geladen sein mußte. Es war jedoch nicht möglich, die 
Äquipotentiallinien in der nächsten Umgebung der Kathode vollständig ein- 
zuzeichnen, da sie sich hier in nächster Nähe der Wand auf einen schmalen 
Bereich zusammendrängten. In manchen Edelgasentladungen, besonders bei 
niedrigen Drucken und dann, wenn die Rohrwand stark bestäubt war, wurde 
in Kathodennähe, 0,1 bis 0,3 mm von der Wand entfernt, ein Potentialwert 
gemessen, der oft nur um einige Volt niedriger als der des Kathodenfalles war. 
Brannten die Entladungen in Gefäßen mit schwach bestäubten Rohrwänden, 
dann standen die Meßergebnisse in zufriedenstellender qualitativer Über- 
einstimmung mit denen von Beck und Emeleus?), die mit einem in die Rohr- 
wand eingeschmolzenen Draht im Fallgebiet auf der Rohrwand ermittelt 
wurden. Hatte aber eine starke Kathodenzerstäubung stattgefunden, dann 
waren die Wandladungen besonders in Kathodennähe größer und der Anstieg 
des Wandpotentials von der Kathode zum Glimmlicht war weniger steil. 
Beck und Emeleus stellten zwischen dem Verlauf des Wandpotentials und 
dem des statischen Aufladepotentials einer in den Fallraum gebrachten Sonde 
eine gewisse Ähnlichkeit fest. Diese Beziehung konnte hier nicht bestätigt 
werden. Es zeigte sich, daß keine Parallelität zwischen beiden Meßwerten 
besteht. In den meisten Fällen war das Wandpotential erheblich höher als 
das statische Aufladepotential einer Sonde. Es darf angenommen werden, 
daß diese hohe Wandladung durch einen Ladungstransport längs einer lei- 
tenden Oberflächenschicht aufrechterhalten wird, wobei die in hoher Kon- 
zentration im Glimmlicht vorhandenen wenig energiereichen geladenen Teil- 
chen einen wesentlich stärkeren Einfluß auf die Höhe der Wandladung haben, 
als die im Fallraum sich hauptsächlich in einer Vorzugsrichtung bewegenden 
energiereichen Ionen und Elektronen. 


2 Der Feldstärkeverlauf im Kathodenfallgebiet 


In den kathodennahen Entladungsteilen wurde durch Differentiation der 
Potentialkurven der Feldstärkeverlauf in der Achse von Edelgas-, Wasser- 
stoff- und Stickstoffentladungen ermittelt. Die Messungen erstreckten sich 
über einen Stromdichtebereich von 3,3 mA/cm?—0,004 mA/em? und einen 
Druckbereich von 0,1 Torr—6 Torr. Die Feldstärkeverteilung wurde bei 
konstantem Druck in der normalen Entladungsform und in anomalen Ent- 
ladungen bestimmt. Bei hohen Stromdichten waren nur dann Messungen 
möglich, wenn der Kathodenfall kleiner als 750 V war, da mit der MeBanord- 
nung höhere Werte nicht erfaßt werden konnten. 

Der zu untersuchende Druckbereich war durch die Druckabhängigkeit der 
Kathodendunkelraumlänge bestimmt. Messungen wurden mit volumen- 
füllender Kathode oder mit einer solchen, deren Durchmesser um 1 cm kleiner 
als der Rohrdurchmesser war, durchgeführt. 

In Abb. 7 wird der Feldstärkeverlauf in einer He-Entladung in Abhängig- 
keit vom Kathodenabstand angegeben. Ein @ bezeichnet die. Grenze des 
Glimmlichtes und ein Stern gibt den Punkt an, an dem die Meßmethode ge- 
wechselt wurde. Der hier gezeigte Feldstärkeverlauf ist charakteristisch für 
alle aufgenommenen Feldstärkekurven. Am anodenseitigen Rand des nega- 


tiven Glimmlichtes wurde in allen Gasen, mit Ausnahme von Stickstoff, der 


®) J. W. Beck, K. G.Emeleus, Philos. Mag. 11,55(1931), 
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Nullwert durchlaufen, in Stickstoff jedoch schon häufig im negativen Glimm- 
licht selbst. Bemerkenswert ist noch die Tatsache, daß stets am kathoden- 
seitigen Rand des Glimmlichtes die Meßmethode gewechselt werden mußte. 
Da in diesem Bereich die Meßgenauigkeit sich verringerte, ist es naheliegend 
anzunehmen, daß der hier auftretende plötzliche starke Feldstärkeabfall nur 
vorgetäuscht wird. Da aber 

kurz vor dem kathoden- fr a} 
abgewandten Ende des Ka- ~ 
thodendunkelraumes das 
Potential noch sicher mit 

der Methode des Verzwei- & 
gungspunktes erfaßt wer- 
den kann und in der Mitte 
des Glimmlichtes mit der 
Nullstrommethode, so folgt 
mit Sicherheit aus dieser 
mit großer Genauigkeit 1 
Iststellkuren Spannungs- Abb. 7. Feldstärkeverlauf in einer He-Entladung. 
differenz, die der oben be- p = 1 Torr, i=1,0mA 

zeichneten - „Fläche . ent- 

spricht, daß die Feldstärkekurve am kathodenseitigen Rand des Glimmlichtes 
mit stärkerer Neigung als im Kathodendunkelraum verläuft. 


400 


23 ¢ 5 6? 8 9 0 %lm] 
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Wegen der Vielzahl der Messungen ist es nur möglich, die Feldstärkekurven 
in bestimmten Gruppen wiederzugeben. Es werden jeweils die bei gleichem 
Druck aufgenommenen Kurven zusammengefaßt. Mit Rücksicht auf die 
später folgenden Vergleiche mit der Beweglichkeitstheorie sind die Feld- 
stärkewerte — gemessen in V/em — durch I”: dividiert, wobei J die Ent- 
ladungsstromstärke in mA bei nichtvolumenfüllender Kathode ist. Der 
Nullpunkt der Abstandsskala ist an den Punkt gelegt, an dem die Feldstärke 
Null oder nur wenig von diesem Wert verschieden war. Bei den Messungen 
mit volumenfüllender Kathode, deren Ergebnisse im gleichen ‚Diagramm 
angegeben werden, sind die entsprechenden Kurven in ihren Unterbrechungen 
durch verschieden gekennzeichnete Kreise kenntlich gemacht. Die dazu- 
gehörigen Stromstärkewerte sind unter Berücksichtigung der verschiedenen 
Kathodenflächen auf die Stromstärken umgerechnet, die bei gleicher Strom- 
dichte durch die kleinere Kathode fließen würde. Bei diesen Umrechnungen 
wurde die durch die Äquipotentiallinien bezeichnete Stromdichteverteilung 
auf der Kathode beachtet. 


Die in allen Abbildungen eingezeichnete Parabel gibt den durch die Be- 
weglichkeitstheorie geforderten Feldstärkeverlauf im Kathodenfallgebiet an. 
Der Vergleich zwischen Theorie und Experiment und die sich darauf ergebenden 
Folgerungen werden an späterer Stelle ausführlich behandelt. aa 


MeBergebnisse in He-, Ne-, A-, Hs- und No-Entladungen 


Im folgenden sind die Meßergebnisse in H,- und N,-Entladungen im Aus- 
schnitt angegeben. Da die in verschiedenen Druckbereichen aufgenommenen 
Kurven sich kaum unterscheiden, seien hier nur die in Edelgasentladungen 
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24 66 0 Rh 6 B 2 4 [na] 
Abb. 10. Feldstärkeverlauf in Ne. p = 0,25 Torr 
Die Raumladungsdichteverteilung 
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bei niedrigstem und héchstem Druck aufgenommenen Kurven wiedergegeben, 
in Wasserstoff- und Stickstoffentladungen nur die im mittleren Druckbereich 


1000 Rot? / 
© 03 mA x10 mA 
012 mA 2 o 
© 1 2 kai 
Abb.8. Feldstärkeverlauf in He. Abb.9. Feldstärkeverlauf in He. 
p=1 p = 6 Torr 
FERN 


Aus den Feldstärkekurven werden mit der Poissonschen Gleichung die 
Raumladungsdichtekurven berechnet. Die Ergebnisse dieser Rechnungen 
sind im Auszug in der Abb. 16 mitgeteilt. Da der Verlauf der Raumladungs- 
dichte in den kathodischen Entladungsteilen der verschiedenen Gase keine 
prinzipiellen Unterschiede zeigt, sind nur die in Ne-Entladungen geltenden 
Werte angegeben. Es muß jedoch beachtet werden, daß die Kurven die 
Ergebnisse zweimaliger Differentiation wiedergeben und daher die Abweichun- 
gen bis zu 10% des angegebenen Wertes betragen können. In unmittelbarer 
Kathodennähe, wo die Fehlergrenze noch größer ist, ist die mit Sicherheit 
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angebbare Mindestgröße der positiven Raumladung eingezeichnet. Da der 
genaue Kurvenverlauf hier nicht bestimmt werden konnte, ist die N 


Kurvenform gestrichelt einge- A. 

In einergewissen Entfernung 

von der Kathode, der in allen 

gen unabhängig vom Druck- 

und Stromstärkebereich ein die 

Potential von 80—120 V gegen 

die Kathode entsprach, änderte 

sich die Neigung der Raum- / / 


ladungsdichtekurven stark. Sie 
fielen weniger steil ab, um 
später am kathodenseitigen 
Rand des negativen Glimm- 
lichtes ein Maximum zu durch- 
laufen. Der Punkt, in dessen 
näherer Umgebung die Neigung 
der Raumladungsdichtekurve 500} 
sich stark änderte, ist mit 


001 mA 
dem dazugehörigen Raumpoten- mA 
tial in die obigen Kurven ein- 
gezeichnet. Mit Sicherheit ist 
den berechneten Kurven zu ent- ee Eu 
nehmen, daß das Potential Abb. 11. Feldstärkeverlauf in Ne. p= 1 Tor 


dieses Punktes sich mit ab- 
nehmender Stromstärke bei konstantem Druck verringerte. Diese Ab- 
hängigkeit ist bei allen untersuchten Gasen bemerkbar — auch bei Wasser- 
stoff und Stickstoff — ; hier nahm jedoch das Potential des oben bezeichneten 


. 
j 
5 0 15 20 25 
Abb. 12. Feldstärkeverlauf in A. p = 0,1 Torr nigh 


Punktes in Stickstoff Werte zwischen 350 V und 160 V an, in Wasserstoff 
zwischen 250 V und 70 V. Ein gesetzmäßiger Zusammenhang zwischen dem 
Druck und diesem Potential konnte bei keinem Gas’festgestellt werden. 


4 


ben, 
eich 
> 
; 
i Va i 
+ 50 mA 
x 10 mA jan 
nm) 
“yg 
le. 
] 
r 
| ; 
F 
en 
ngs- 
eine 
den 
die 
lun- 
arer 
heit 4 


32 Ammalen der Physik. 7. Folge. Band. 058 


Diskussion 


Die Meßergebnisse wurden mit der von Warren”) angegebenen Beweg- 
lichkeitstheorie verglichen. Warren geht bei der Berechnung des Feld- 
stärkeverlaufs im Kathodenfallgebiet von der Annahme aus, daß keine Ioni- 


£ £4 sation durch ElektronenstoB erfolgt, 
iF PR wohl aber die Ionen in diesem Ent- 
ladungsteil Zusammenstöße erleiden. 
Ihre Geschwindigkeit in Abhängigkeit 
\ von der Feldstärke wird durch 
k’ (E |p)" 
| gegeben. 
\ \y 
\ 
2000} 
# 
| 
mA 
a } 
o7 
a Abb. 13. Feldstärkeverlauf in A. Abb. 14. Feldstärkeverlauf in H,. 
p= 4 Torr Te p = 2,5 Torr 


Mit diesen Annahmen ergibt sich als Lösung der Poissonschen Gleichung 
unter Vernachlässigung des Beitrages der Elektronen zur Raumladung 


"ls 

wobei F die Größe der von der Entladung bedeckten Kathodenfläche in cm*, 
k’ den Beweglichkeitswert der positiven Ionen und J die Entladungsstrom- 
stärke angibt. 

Aus dieser Gleichung ist ersichtlich, daß, wenn die Beweglichkeitstheorie 
die Verhältnisse im Fallraum richtig beschreibt, alle bei konstantem Druck 
aufgenommenen Feldstärkekurven mit der von der Beweglichkeitstheorie 
geforderten Kurve zusammenfallen müssen, trägt man E/I*/s über den Ab- 
stand vom Glimmlicht auf. 

Beim Vergleich der experimentell ermittelten Feldstärkekurven mit dieser 
Theorie wurde ein k’-Wert angenommen, der eine möglichst gute Anpassung 
der durch die Theorie vorgeschriebenen Kurvenform an die experimentell 
ermittelte ergab. Die Berechnungen werden unter der Annahme durchgeführt, 
daß der in Gl. (1) auftretende Faktor 1 + y = 1 ist. Diese Näherung ist, wie 


10) R. Warren, Physic. Rev. 98, 1658 (1955). 
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Messungen von verschiedenen Autoren zeigen, durchaus gerechtfertigt. Die 
veg. | 80 bestimmten k’-Werte unterschieden sich in Übereinstimmung mit den 
ae Meßergebnissen von Warren?) in Edelgasentladungen von den durch Horn- 
ae beck?) für Atomionen gefundenen Werten ungefähr um den Faktor zwei. 
gt, Sie seien in nachfolgender Tabelle zusammen mit den Hornbeckschen 


Ent- In 


E =10Torr 


005} 
mA o 15 mA 20 1 2 3 4% 5 6 
u Abb. 15. Feldstärkeverlauf in N, Abb. 16. Raumladungsdichte- ie: 
p=1Torr  verlauf in Ne. p=10Torr 
Bi. 


Beweglichkeitswerten, die für E/p-Werte größer als 100 V/em mm Hg gültig 
ung | sind, mitgeteilt. In der eben genannten Arbeit sind keine Beweglichkeitswerte 
für Molekülionen angegeben, doch wurde darauf hingewiesen, daß ihr Wert 
ungefähr doppelt so groß ist wie der der Atomionen. 


(1) 
Tabelle 1 
‚m?, k’-Werte in em’ sec!» (mm Hg)'* - 10% (T = Torr) 
om- 
k ren Nt a nbeck Aus dem Feldverlauf ermittelte Werte 
orie 
uck | He 4 p=6T |p=4T p =1T 
orie k = 6,5 k=6 k’ = 8,7 
Ab-] Ne 1,6 p=1T |p =05T | p =0,.26T 

en A 0,8 p=4T |p=1T | p=08T | p =09T 
= |k=16 = 1,6 b= 16 
itell 
rt, Die aus den Meßwerten berechneten Beweglichkeiten kommen den der 
wie 


Molekülionen näher als den der Atomionen. Die sich hieraus ergebende 


1) J. A: Rev. 84, 616 
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Folgerung, daß der Stromtransport der positiven Teilchen hauptsächlich en 
von Molekülionen getragen wird, erscheint gewagt, doch soll ihre Berechtigung | Bei k 
näher untersucht werden. nähe 
Massenspektroskopische Untersuchungen von Morris™) in Edelgasent- | gag t 
ladungen zeigen, daß im Bereich der positiven Säule bei Drucken von 1Torr | 99er 
bis zu 8% des Stromes der positiven Träger zur Wand von Molekülionen ge- | gen } 
tragen werden. Aus den Ergebnissen dieser Untersuchungen kann geschlossen | ‚treel 
werden, daß die Molekülionen hauptsächlich bei dem Zusammenstoß hoch Entla 
angeregter Atome mit neutralen Atomen gebildet werden, ebenfalls auch, | gieseı 
aber nicht mit so großer Wahrscheinlichkeit, durch den Dreierstoß von Atom- | gie D 
ionen und neutralen Atomen. Da im negativen Glimmlicht die Voraussetzun- | trone 
gen für den Ablauf dieser Bildungsprozesse außerordentlich günstig sind, | yie ¢ 
darf man erwarten, daß hier Molekülionen in großer Zahl gebildet werden. Es E 
ist daher durchaus möglich, daß diese wegen ihrer größeren Beweglichkeit geste 
hauptsächlich den positiven Ionenstrom im Fallraum von Edelgasentladungen | dicht 
bestimmen. stärk 
Tabelle 2 tialb 
k’-Werte in sec”! (mm Hg)/*-10¢ (T= Torr) Kath 
Aus dem Feldverlauf ermittelte Werte m 
Bew: 
H, p =6T p =2,5T » =1T | p =0,5T dure 
k’ = 35 = 2 k’ = 22,3 k’ = 20 
oder 
N, 0,97 p =2,5T eth p =05T | p =0,25T kanr 
2 
Die für das Stickstoffmolekül angegebene Beweglichkeit wurde von Hor 
Varney®) für E/p-Werte, die größer als 100 V/em mm Hg waren, gemessen. mm] 
Beweglichkeitswerte von Wasserstoffmolekiilen in dem hier untersuchten zutr 
E/p-Bereich waren leider nicht greifbar. Auffallend ist bei beiden Gasen die E|p- 
Tatsache, daB die Beweglichkeitswerte mit abnehmendem Druck sich ver- Bere 
ringern und daß eine erhebliche Differenz zwischen den von Varney für N, Ent 
gemessenen und den hier bestimmten Werten besteht. Diese Besonderheiten war 
können dadurch hervorgerufen worden sein, daß in Stickstoff neben N¢- Teil 
auch Nj-Ionen, die von Varney bei seinen Beweglichkeitsmessungen beob- ] 
achtet wurden, auftraten und die Beweglichkeitswerte in der oben genannten Rau 
Weise beeinfluBten. In Wasserstoff könnten entsprechende Ionenarten dan 
ebenfalls von Bedeutung sein. mer 
Zusammenfassend darf gesagt werden, daß aus den Meßergebnissen ge- um 
schlossen werden kann, daß die Beweglichkeitstheorie in allen untersuchten Ver 
Druckbereichen in der ersten Hälfte des Kathodenfallgebietes zwischen Glimm- son: 
licht und Kathode gültig ist. | 
Diese Meßergebnisse stehen teilweise in Widerspruch zu den von Warren?) unv 
ermittelten Gesetzmäßigkeiten. Warren untersuchte den Feldstärkeverlauf sch: 
in anomalen Glimmentladungen in He, H,, N, und A in einem Druckbereich tho 
von 0,03—1 Torr mit einem magnetisch kompensierten Elektronenstrahl. auf: 
#2) D. i . Soc. A 
RN Varney, Payne. Rev: 0,10 0908). 
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Er fand, daß bei hohen Drucken bei den leichteren Gasen die Beweglichkeits- 
theorie die Verhältnisse im gesamten Kathodenfallgebiet richtig wiedergibt. 
Bei kleineren Drucken traten Abweichungen auf. Hier wurde in Kathoden- 
nähe ein zu steiler Anstieg der Feldstärke verzeichnet. Hierzu ist zu bemerken, 
daß bei den höheren Drucken die Kathodendunkelraumlänge stets kleiner als 
2,2 cm war und der Bereich, in dem der steile Feldstärkeverlauf entsprechend 
den hier gefundenen Werten auftreten müßte, sich nur über einige mm er- 
strecken dürfte. Es ist zweifelhaft, ob mit einem Elektronenstrahl, der im 
Entladungsfeld stets von der Kathode weggelenkt wird, in Kathodennähe 
dieser Verlauf festgestellt werden kann. Mit abnehmendem Druck wächst 
die Dunkelraumlänge, und es verbesserte sich die Möglichkeit, mit dem Elek- 
tronenstrahl den steilen Anstieg der Feldstärke in Kathodennähe festzustellen, 
wie die Meßergebnisse von Warren zeigen. 

Es bleibt noch zu untersuchen, wodurch der bei allen Entladungen fest- 
gestellte steile Anstieg der Feldstärke und damit auch der Raumladungs- 
dichte in Kathodennähe hervorgerufen wird. In allen Druck- und Strom- 
stärkebereichen wurde in Edelgasentladungen festgestellt, daß in einem Poten- 
tialbereich von 80—120 V die Raumladungsdichtewerte in steilem Anstieg in 
Kathodenrichtung sich stark erhöhten und in Kathodennähe Werte annahmen, 
die die des Dunkelraumes um ein Beträchtliches übertrafen. Diese Raum- 
ladungsdichtezunahme kann folgende Ursachen haben: Zunächst kann der 
Beweglichkeitswert in diesem E/p-Bereich abnehmen, es kann aber auch 
durch eine kurzwellige Strahlung, die aus kathodennahen Entladungsteilen 
oder von der Kathode ausgesandt wird, eine Trägerbildung erfolgen. Auch 
kann Ionisation durch Elektronenstoß oder durch den Stoß positiver Ionen 
auftreten. 

Zur ersteren Annahme, die von Warren für E/p-Werte, die den von 
Hornbeck untersuchten E/p-Bereich überschreiten, also größer als 1000 V/em 
mmHg sind, ausgesprochen wurde, darf gesagt werden, daß sie sicher nicht 
zutrifft, da der Raumladungsdichteanstieg auch dann auftritt, wenn die 
E/p-Werte in Kathodennähe noch in dem von Hornbeck aufgenommenen 
Bereich enthalten sind. Ionisation durch kurzwellige Strahlung, die aus der 
Entladung kommt, ist nicht sehr wahrscheinlich, da nicht einzusehen ist, 
warum sich das Ionisierungsgebiet ausschließlich auf einen scharf begrenzten 
Teil des Kathodendunkelraumes beschränken soll. 

Ionisation durch den Stoß positiver Ionen kann ebenfalls nicht für den 
Raumladungsdichteanstieg verantwortlich gemacht werden, da dieser auch 
dann auftritt, wenn die Energie der positiven Ionen nicht ausreicht, um eine 
merkliche Ionisation hervorzurufen. War der Kathodenfall hoch genug, 
um eine Ionisation durch positive Träger zu ermöglichen, dann zeigte der 
Verlauf der Raumladungsdichtekurven keine Veränderungen gegenüber der 
sonst auftretenden Form. 

Ionenbildung durch Elektronenstoß nur vor der Kathode ist ebenfalls 
unwahrscheinlich, denn aus dem Verlauf des atomaren Ionisierungsquer- 
schnittes in Abhängigkeit von der Elektronenenergie folgt, daß auch in ka- 
thodenfernen Bereichen des Kathodendunkelraumes starke Ionisation 
auftreten müßte, da in Edelgasen für Elektronenenergien größer als 100 eV 
der Ionisierungsquerschnitt in dem hier untersuchten Bereich praktisch 
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soviel Energie durch Stöße verbunden mit Anregung und Ionisation ver- 
lieren, daß sie später nicht ionisieren können, kann widerlegt werden, da in 
den Gebieten geringerer Raumladung in den meisten Fällen noch Spannungs- 
differenzen zwischen dem Glimmlicht und den entsprechenden Punkten auf- 
treten, die 100 V weit überschreiten. Auch ergaben Messungen von Brewer 
und Westhaver!#4), daß ein großer Teil der von der Kathode ausgesandten 
Elektronen das Glimmlicht nahezu mit Kathodenfallenergie erreicht. 

Nimmt man an, daß in Kathodennähe durch den vorher beschriebenen 
Bildungsprozeß Molekülionen entstehen, dann kann die starke positive Raum- 
ladung vor der Kathode erklärt werden; denn bei den Molekülionen zeigt der 
Verlauf des Bildungsquerschnittes in Abhängigkeit von der Elektronen- 
energie eine Ähnlichkeit mit dem der Anregungsquerschnitte verschiedener 
Linien des entsprechenden Atoms. Die Bildungsquerschnitte durchlaufen, 
wie Messungen von Hornbeck und Molnar”) zeigen, in einem Bereich von 
20—40 eV in He, Ne und A ein Maximum. Bei 60 eV ist ihr Wert in He und A 
schon auf !/, und in Ne bei 80 eV auf ?/, ihres Maximalwertes gesunken. Die 
hier gefundene Raumladungsdichteverteilung ist daher mit diesem Verhalten 
gut zu vereinbaren. Auch spricht für die obige Annahme die Tatsache, daß 
in den untersuchten Entladungen die Ausdehnung des Bereiches hoher Raum- 
ladung stets mit dem der ersten Kathodenschicht zusammenfällt. 

Auch in Wasserstoff- und Stickstoffentladungen war ein derartiger Zu- 
wachs der Raumladung vor der Kathode zu verzeichnen. Welche Prozesse 
für diesen Raumladungsdichteverlauf verantwortlich gemacht werden müssen, 
kann hier nicht entschieden werden. Mit Sicherheit ist nur festzustellen, 
daß der Verlauf des Ionisierungsquerschnittes für Molekülionen wegen seiner 
Ähnlichkeit mit dem des Ionisierungsquerschnittes für Atomionen in Edelgasen 
einen derartigen nur auf die Kathodennähe beschränkten Raumladungs- 
dichtezuwachs nicht zuläßt. Man darf vielleicht vermuten, daß hier kompli- 
ziertere Ionenarten eine Rolle spielen. 


Herrn Professor Dr. R. Seeliger möchte ich für die Anregung zu dieser 
Arbeit und sein stetes Interesse für ihren Fortgang herzlich danken. Herrn 
Fricke bin ich für die Unterstützung und die technischen Ratschläge beim 
Aufbau der Versuchsanordnung zu Dank verpflichtet. ir rk 
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wh Zur Theorie der Relaxation 


I. Die thermodynamischen Grundlagen und ihre mathematische Lésung 


Von Siegfried Kästner 


Inhaltsübersicht 


In der Arbeit werden zunächst ganz allgemein die thermodynamischen 
Grundgleichungen der Relaxation aufgestellt, wobei neben der reinen Re- 
laxation gleichzeitig auftretende Effekte der inneren Reibung Berücksichtigung 
finden. Bei der mathematischen Formulierung der Theorie wird eine Tensor- 
darstellung im Raum der Zustandsvariablen verwendet, die sich auch bei der 
anschließend behandelten Lösung der Grundgleichungen bewährt. Die Er- 
gebnisse werden tabellarisch zusammengestellt. Schließlich wird ihre An- 
wendung auf die mechanische und dielektrische Relaxation diskutiert. 


nies: 
4. Einleitung 

Bei der gesetzmäßigen Erfassung des makroskopischen Verhaltens der 
Materie spielen neben den Gleichgewichtszuständen die Zeiteffekte 
eine wichtige Rolle!). Von diesen Zeiteffekten läßt sich eine große Gruppe 
unter dem Begriff Relaxationserscheinungen. zusammenfassen. Dabei 
handelt es sich um irreversible Vorgänge, bei denen in dem betrachteten System 
beim Anlegen äußerer Kräfte (allgemeiner Art) auf Grund der inter- und 
intramolekularen thermischen Wechselwirkung Rückstell- und Reibungs- 
kräfte auftreten, während die mikroskopischen Trägheitseffekte vernach- 
lässigbar klein sind. 

Bekannte Erscheinungen dieser Art sind z. B. die mechanische und die 
eng mit ihr verkoppelte akustische Relaxation, die dielektrischen Eigen- 
schaften polarer Substanzen, das para- und ferromagnetische Relaxations- 
verhalten und die Kernresonanz-Relaxation. Außerdem fällt hierunter auch 
die Kinetik chemischer Reaktionen und verwandter Vorgänge wie Löslich- 
keit, Ionisation, Dissoziation, Luminiszenz usw. 

Die folgenden Ausführungen dienen speziell der Aufklärung des mecha- 
nischen und dielektrischen Relaxationsverhaltens der Materie. Sie 
werden aber thermodynamisch so allgemein gefaßt, daß sie auch für andere 
Relaxationserscheinungen, bei denen keine magnetischen Felder wirksam 
sind?), Gültigkeit besitzen. Dabei wird neben der reinen Relaxation das 


dr 


1) W. Walcher, Kolloid-Z. 184, 1 (1953). 

2) Beim Auftreten magnetischer Felder sind gewisse Symmetriebeziehungen, die 
bei den folgenden Ableitungen eine wichtige Rolle spielen, nicht erfüllt. 
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gleichzeitige Auftreten von innerer Reibung behandelt. ai findet in reiner 
Form z. B. im mechanischen Fall als Viskosität und im elektrischen Fall als 
Ohmsche Leitfähigkeit ihren Ausdruck und kann im Gegensatz zu weiteren 
Transporterscheinungen wie Diffusion und Wärmeleitung nicht durch nach- 
träglich hinzugenommene Gleichungen berücksichtigt werden. 

Vom mathematischen Standpunkt fügen sich die Beziehungen der inneren 
Reibung zwanglos in die Theorie der reinen Relaxation ein und ergeben eine 
Bereicherung der auftretenden Spektren hinsichtlich gewisser Grenzterme. 
Im folgenden soll deshalb auch bei Berücksichtigung der inneren Reibung 
nur von Relaxation gesprochen werden, wobei gegebenenfalls für die Relaxation 
ohne innere Reibung der Begriff reine Relaxation verwendet wird. 

Die Ausführungen beschränken sich auf homogene Systeme. Außerdem 
wird vorausgesetzt, daß die Abweichung vom thermodynamischen Gleich- 
gewicht stets hinreichend klein ist. Daraus folgt nach Meixner’) sowohl 
die Anwendbarkeit der Theorie der irreversiblen Prozesse*) als auch die 
Gültigkeit des Boltzmannschen Superpositionsprinzips. Damit kommt man 
zur linearen Theorie der Relaxation, deren thermodynamische Grund- 
lagen im folgenden als Tensorgleichungen im Raum der Zustandsvariablen 
dargestellt werden. 

Für den Spezialfall der mechanischen Relaxation ohne innere Reibung 
gehen die hier abgeleiteten Beziehungen in die bekannte Theorie der reinen 
Relaxation von Meixner®) über. Betrachtet man andererseits mechanische 
Vorgänge, die isotherm verlaufen, so erhält man Übereinstimmung mit den 


von Biot®) angegebenen Formeln. rt 
(ewe 3 
2. Die Grundgleichungen der reinen Relaxation Fe 


Die thermodynamische Theorie der reinen Relaxation nimmt ihren 
Ausgang von der inneren Energiedichte U des betrachteten Systems, die eine 
Funktion von n Zustandsvariablen 2* ist. Dabei wird zunächst jeder Zustands- 
variablen 2* eine konjugierte Variable Z, zugeordnet, die nach Gibbs 
durch das totale Differential der Energiedichte U definiert ist?) : A re 


aU 
und demzufolge aus der Beziehung 


berechnet werden kann. Die n Zustandsvariablen gliedern sich im einzelnen 
ch a äußere Zustandsvariable a* und i innere Zustandsvariable i*. 


2* = a* 0=* 


2,3 
27 el 0 = a* »=a+l,a+2...n (2,5) 

3) J. Meixner, Kolloid-Z. 184, 3 (1953). 

4) I. Prigogine, Etude Thermodynamique des Phénoménes irreversibles. Paris— 
Liege 1947; S.R.de Groot, Thermodynamics of Irreversible Processes. Amsterdam 
1951; K. G. Denbigh, The Thermodynamics of the steady state. London—New York 
1951; R. Haase, Ergebn. exakt. Naturwiss, 26, 56 ET R. O. Davies, Rep. 


Progr. 19, 326 (1956). 
5) J. Meixner, Z. Naturforschg. 9a, 654 (1964), 
6) M.A. Biot, ‘I. appl. Phys. 26, 1385 ash 
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Für die konjugierten Variablen Z, besteht — abgesehen von einem auftre- 


tenden negativen Vorzeichen — eine analoge Aufspaltung in äußere kon- 
jugierte Variable A, und innere konjugierte Variable I,. 
= A, . 0 = I, 1,3,...6 
(2,4) 
Z,=-L <A, x=at+la+2,..." 
Die Gl. (2,1) läßt sich damit auch in folgender Form schreiben: 
a n 
dU= 3 A,da*— 2 I, di*. (2,5) 
| x=a+1 


Spezialisiert man weiterhin hinsichtlich der Reihenfolge der Zustandsvariablen 
a! als die Entropiedichte S mit der Temperatur T als konjugierte Variable A, 


| 
so geht Gl. (2, 5) in die Theres 


über. Dabei stellt die erste Summe der Gl. (2,7) die dem System zugeführte 
differentielle Arbeitsdichte dA dar. 

Die zweite Summe G. (2,7) beschreibt die im 
System. Unter Verwendung des ersten Hauptsatzes der Thermodynamik 

(2,9) 


kommt diese Tatsache in der Beziehung u jae 


n 
TdS=dQ+ 2 (2,10) 


zum Ausdruck, wobei Q die dem System zugeführte Wärmedichte bezeichnet. 
Die inneren konjugierten Variablen J, verschwinden in einem als Bezugs- 
zustand gewählten ungehemmten Gleichgewichtszustand Ip,,%). 

Tox = 0. (2,11) 
In der Umgebung des Bezugszustandes kann man die innere Energiedichte U 
in eine Reihe entwickeln. Zu diesem Zweck führen wir die Abweichungen der 
Zustandsvariablen z* vom Bezugszustand 2 als die Verschiebungen des 
Systems q* 

ein, die sich entsprechend Gl. (2,3) in die äußeren Verschiebungen .x* und die 
inneren Verschiebungen y* gliedern. 

= x O= 1, 3,...6 

= y" 0 = 2 x=at+l1,a+2,. 
Andererseits nennen wir die entsprechend gebildeten Differenzen der kon- 
jugierten Variablen Z, die Kräfte des Systems iso 


(2,13) 
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mit den äußeren Kräften K, und den inneren Kräften L,. twang 
P, = K,, 0 = L,, 
P,=-L 0=K, »=a+la+2...n. 


Außerdem wird angenommen, daß sämtliche Bezugswerte der äußeren 
Verschiebungen x* und äußeren Kräfte K, außer x! und K, verschwinden 
sollen. 


(2,15) 


Kox = 0 (2,16) 
Unter dieser Voraussetzung lautet die Entwicklung der inneren Energie- 
dichte U | 
U=-U,+U, +0, (2,17) 
mit dem konstanten Glied U,, dem linearen Glied U, ” re Py 
U, = (2,18) 


das die Temperatur 7, des Bezugszustandes enthält, und der quadratischen 
Form U, 

“ey 


x=liA= 


(2,19) 


aus der man nach Formel (2,2) unter Verwendung der Gin. (2,6) und (2,16) 


die Kräfte des Systems P, berechnen kann. _ ’ 
eq” ite 


Weitere Glieder der Entwicklung (2,17) werden im Rahmen der linearen 
Theorie vernachlässigt. 
Die Matrix A,, ist symmetrisch 


(2,20) 


Axa = Aix, (2,21) 
da ganz allgemein gilt 


Die zugehörige quadratische Form U, ist nach Meixner) positiv definit. 
Dabei ergibt sich zunächst aus der Tatsache, daß die innere Energiedichte 
für konstante Entropiedichte ein absolutes Minimum besitzt, daß U, für 
q=0 positiv definit ist. Aus der gleichen Eigenschaft der freien Energie- 
dichte F für konstante Temperatur findet man anschließend, daß F,(T,) 

nn A, 4A 

un) 
ie > Au 
ebenfalls positiv definit sein muß. Schließlich kann man U, ganz allgemein 
durch die Beziehung 


F, (T,) =+ (2,23) 


P 
darstellen, woraus folgt, daß U, in sämtlichen Variablen g* positiv definit 
ist, da A,, immer verschieden von null und positiv sein muß. Die letztere Aus- 
sage ist eine Folge davon, daß die Dichte 0, und die Temperatur T, des Bezugs- 
zustandes sowie die spezifische Wärme c bei konstanten äußeren Verschie- 


bungen 2* (x > 2) ebenfalls immer positiv und endlich sind, wobei der Zu- 
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sammenhang zwischen A,, und den angegebenen Größen allerdings erst im 
Laufe späterer Rechnungen deutlich wird. 


Eine genaue Untersuchung zeigt jedoch, daß U, auch positiv semidefinit 
sein kann. Dies bedeutet, daß die quadratische Form, die zunächst eine 
Funktion von n Verschiebungen ist, durch eine geeignete lineare Transfor- 
mation nur noch von rn — k Verschiebungen abhängig wird. Im Falle der 
reinen Relaxation hat dieses Ergebnis wenig Bedeutung, da man die über- 
flüssigen Verschiebungen eliminieren kann und schließlich wieder zu einer 
(n — k)-reihigen Matrix A,, mit positiv definiter quadratischer Form U; 
gelangt. Im Falle der später durchgeführten Verallgemeinerung der Theorie 
ist dies jedoch nicht möglich und muß besonders beachtet werden. 


Nach dieser für die mathematische Behandlung wichtigen Betrachtung 
über die Eigenschaften der Matrix A,, werden nunmehr die irreversiblen Vor- 
gänge des Systems näher untersucht. Nach der Theorie der irreversiblen 
Prozesse) bestehen hierbei lineare Beziehungen zwischen den ersten Ab- 
leitungen nach der Zeit der inneren Zustandsvariablen i* und den konjugierten 
inneren Variablen J,, die im folgenden unter Verwendung der Gln. (2,11) 
bis (2,15) unmittelbar für die entsprechenden inneren Verschiebungen y* 


und inneren Kräfte L, angegeben werden sollen. 


Dabei wurde die Ableitung nach der Zeit ¢ durch einen Punkt gekennzeichnet. 
Die Matrix E** ist nach der Onsagerschen Reziprozitätsbeziehung symme- 
trisch 4) 


(2,25) 


‘T, ET — = (2, 26) 
worin die Reversibilität der mikroskopischen Vorgänge zum Ausdruck kommt’). 
Außerdem ist nach Meixner ) die als Dissipationsfunktion D bezeichnete 
zugehörige quadratische Form 


Lu In (2,27) 
positiv definit. Es folgt dabei zunächst, daß D nicht negativ sein kann, da 
die Entropiedichte eines adiabatisch isolierten Systems nicht abnehmen kann, 
denn es gilt nach Gl. (2,10) 
TS=Q+2D. (2,28) 
Wäre D jedoch nur positiv semidefinit, so würde die Determinante |E*?| 
verschwinden und die j* wären linear abhängig, so daß eine lineare Kom- 
bination der inneren Zustandsvariablen i* zeitlich konstant und verschieden 
von null sein könnte. D muß also positiv definit sein. Schließlich folgt daraus, 
daß auch die Matrix B,,, der inversen Beziehung 


ya 
L, = m Bya y ER (2,29) 


7) N.G. van Kampen, Fortschr. Physik 4, 405 (1956). 
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symmetrische Gestalt besitzt wh desse 
Bar = Bax (2,30) zum 
und die zugehörige quadratische Form itil Ayı 

ebenfalls positiv definit ist. wi Agi 
Kine genaue Untersuchung zeigt jedoch, daB die unter Verwendung der « 
Matrix B,; gebildete quadratische Form der Gl. (2,31) ganz analog dem ent- u 

sprechenden Verhalten der quadratischen Form der Energiedichte (2,19) 

im Spezialfall auch positiv semidefinit sein kann. I 

Unter Verwendung der quadratischen Form D kann man schließlich die 
Beziehung (2,29) in folgender Weise schreiben: A, 
aD 
L, =—. (2,32) 

‘ 
Damit sind die Grundgleichungen der reinen Relaxation gewonnen. Die ach 
formelmäßigen Ergebnisse sollen im folgenden kurz zusammengestellt werden. beh 
Dabei geht man von den quadratischen Formen U, und D aus. gan 
= erze 


Zwischen den Kräften und den Verschiebungen bestehen lineare Differential- 
gleichungen, die durch die Beziehungen (2,20) und (2,32) gegeben sind. 


au. éD We 
P,= L, = (2,34) 
Eliminiert man die inneren Kräfte L,, so kann man dafür auch folgendes 
Gleichungspaar angeben: 
aU, aU, , 
=— —=0 2,35 i 


Berücksichtigt man schließlich, daß D von den &* unabhängig ist und die 
äußeren Kräfte K, nach Gl. (2,15) für <> a verschwinden, so lassen sich de 
_ die zwei Formeln der Gl. (2,35) durch eine in Lagrangescher Form auf- Fc 

_ tretende Beziehung zusammenfassen. 


— 
aU, 


Führt man die Differentiation aus, so folgt das Gleichungssystem b merits 


tli 
= Axa K, 1 s Ss ot au 
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dessen spezielle Gestalt am deutlichsten in der folgenden Matrizendarstellung 


zum Ausdruck kommt. 


tri 
| 
Ar --- Aen | |x, 
Aa+11-=- Ir! Ba+1,a+1--- Ba+ı,n] +! 0 
A,ı Ang B, a+1 0 
(2,38) 


3. Die Grundgleichungen der reinen inneren Reibung 


Nachdem im letzten Kapitel die Grundgleichungen der reinen Relaxation 
aufgestellt wurden, soll nunmehr der Fall der reinen inneren Reibung 
behandelt werden. Dabei gilt nach der Theorie der irreversiblen Prozesse*) 
ganz entsprechend Gl. (2,10) (für x > 2) eine Beziehung für die Entropie- 
erzeugung, die wegen Gültigkeit von Gl. (2,16) sofort für die im letzten Ka- 
pitel eingeführten äußeren Kräfte K, und die äußeren Verschiebungen 2* 

TS=G+ 
Weiterhin besteht für «> 2 ein zur Gl. (2,29) analoger Zusammenhang 
zwischen den äußeren Kräften K, und den Verschiebungsgeschwindigkeiten <* 


a 
K, == By, (3,2) 
= 
mit der symmetrischen Matrix B,, — ont 
= Bix (3,3) 


der unter Verwendung der positiv definiten bzw. semidefiniten quadratischen 
Form D 


D= + = &* (3,4) 
auch in der Gestalt 


ikea werden kann. Dagegen wird K, wie im Fall der reinen Relaxation 
aus der Form der U, 


= 
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wobei hier U, jedoch nur von 2! abhängt. sh 


1 
Ay (=). (3,7) 


4. Die Grundgleichungen der Relaxation (unter Berücksichtigung der inneren 
Reibung) 


Im Falle der reinen Relaxation — der im folgenden durch einen Stern 
gekennzeichnet wird — werden die äußeren Kräfte X, nach Gl. (2,35a) be- 
stimmt 


Kr = (4,1) 
dx” 
während im Falle der reinen inneren Reibung — der entsprechend zwei 


Sterne trägt — die äußeren Kräfte X, nach den Gin. (3,5) und (3,6) berechnet 
werden. 


(4,3) 


0x” 


Treten Relaxation und innere Reibung gleichzeitig auf, so addieren sich beide 
Anteile hinsichtlich der äußeren Kräfte 


K,, = K* + K**, 4) 
was auf die Beziehung 
0x” 08” 


führt, wobei die quadratische Form der Energiedichte U, durch die Summe 
U, = US + UM (4,6) 
gebildet wird. Nimmt man die Beziehungen (2,35b) hinzu, welche die inneren 


Vorgänge beschreiben, so lauten die allgemeinen Grundgleichungen der 
Relaxation schließlich formal genau wie im Falle der reinen Relaxation 


+ — = K, 
mit der entsprechend Gl. (4,6) gebildeten Dissipationsfunktion D. 
D = D* + D**., (4,8) 


In der Matrizendarstellung der Gl. (2,38) bedeutet dies, daß man an Stelle 


von der bei der reinen Relaxation verwendeten Matrix 


0 
(4,9) 


Bo +1.041--- Ba+ı,n 
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TEM) Ba+1,a+1 --- Ba+ı,n 
+1 Bun 


ausgehen müßte. Berücksichtigt man jedoch. daß zwischen irreversiblen 

Vorgängen verschiedener Art nach der Theorie der irreversiblen Prozesse®) 

ganz allgemein eine Kopplung auftritt, so wird man schließlich die symme- 


| 0 

| = 

benutzen, bei der nur noch in der ersten Zeile und Spalte Nullen then. 
Die Dissipationsfunktion D wird also nunmehr durch die Beziehung Ye 
D=+ 22 Bad (4,12) 


gebildet. Fiihrt man die Differentiation u Gl. (4,7) aus, so folgt im einzelnen 


n 


n n 

Ang + Bag =0 atli<ucn 


analog zur Matrizendarstellung. 

Hinsichtlich der Eigenschaften der quadratischen Formen U, und D wird 
dabei entsprechend der im Kapitel 2 durchgeführten Untersuchungen voraus- 
gesetzt, daß sie positiv definit oder positiv semidefinit sein können. Allerdings 
bestehen für das semidefinite Verhalten gewisse naheliegende Ausschließungs- 
bedingungen zwischen U, und D, die zum Teil noch im einzelnen erwähnt 
werden. 

Schließlich wird noch der Spezialfall zugelassen, daß die Dissipations- 
funktion D verschwindet, um auch un itinmaaatnenetts Medien zu erfassen. 
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5. Tensordarstellung der Grundgleichungen der Relaxation 


Die im letzten Kapitel abgeleiteten Grundgleichungen der Relaxation 
lassen sich als Tensorgleichungen in einem n-dimensionalen Raum deuten, 
falls man die n Verschiebungen g* als Komponenten eines Vektors ansieht. 
Diese Tensordarstellung bringt neben einer übersichtlichen Formulierung den 
Vorteil, daß man bei den weiteren Rechnungen von wichtigen Transforma- 
tionsgesetzen des allgemeinen Tensorkalküls®) Gebrauch machen kann. 
Über die Metrik des zugrunde liegenden Raumes bzw. sogar über seine 
affinen Eigenschaften sind dabei keine speziellen Annahmen notwendig. 
Um jedoch gleichzeitig gewisse anschauliche Vorstellungen mit den mathe- 
matischen Untersuchungen verknüpfen zu können, wird speziell angenommen, 
daß die Verschiebungen q* die (kontravarianten) kartesischen Komponenten 
eines Vektors in einem n-dimensionalen euklidischen Raum sind, der zur Ab- 
kürzung der Verschiebungsraum (V-Raum) genannt werden soll. 
In Tensordarstellung lauten die Grundgleichungen der Relaxation 

Ay g’ + Ky, (5,1) 
wobei entsprechend den allgemeinen Rechenregeln des Tensorkalküls die 
griechischen Indizes die Werte 1 bis n annehmen können und über einen be- 
stimmten Index, der gleichzeitig oben und unten auftritt, summiert werden 
muß. Die bisherigen Matrizen A,, und B,, treten als symmetrische ko- 
variante Tensoren auf, während die äußeren Kräfte X, die Komponenten 


eines kovarianten Vektors sind. Dabei ist zu beachten, daß in dem Ausgangs- 
koordinatensystem die Beziehungen 


| 
und 


(5,3) 


gelten. Die quadratischen Formen 


ck 


D =+ By; q* g* (5,5) 


sind skalare Größen im V-Raum. Dabei werde die Dissipationsfunktion D 
nunmehr im Gegensatz zur Definition (4,12) als Funktion sämtlicher g* 
betrachtet, wobei sie — da B,, immer verschwindet — im allgemeinsten Fall 
positiv semidefinit oder null ist. 


6. Periodische Vorgänge 
Zur Gewinnung der allgemeinen Lösung der Gl. (5,1) werde von periodi- 
schen Vorgängen ausgegangen. Dabei setzt man 
= Re {q* e') K, = Re 2%, (6,1) 
= ow obei Re den Realteil kennzeichnet. z ist eine komplexe Konstante, die im 


Spezialfall periodischer Vorgänge mit zeitlich konstanter Amplitude rein ima- 
ginär wird:2=iw. Dann ist = 22 f die Kreisfrequenz und f die Frequenz. 


8) A. Benmeorichd, Mechanik der deformierbaren Medien. Leipzig 1949. S. 358; 
E. Schrödinger, Space- Time Structure. Cambridge 1950. 
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Für die komplexen Amplituden der Verschiebungen q* und der Kräfte 
8, folgt mit dem Ansatz (6,1) aus Gl. (5,1) die ee 


mit dem komplexen Tensor W,; Bow 

= Ana + 2 Ban. Shab (6,3) 


7. Die Verschiebungen als Funktion der äußeren Kräfte 

Zur Berechnung der Verschiebungen als Funktion der äußeren Kräfte 
wird im V-Raum ein neues Koordinatensystem eingeführt, das die beiden 
symmetrischen Tensoren A,, und B,, mit positiv definiter bzw. semidefiniter 
quadratischer Form auf Hauptachsen transformiert. Ein derartiges Koor- 
dinatensystem läßt sich durch eine lineare Transformation immer finden, 
denn es ist anschaulich verständlich, daß man zunächst durch eine reine Drehung 
einen der beiden Tensoren auf Hauptachsenform bringen kann. Durch eine 
anschließende Dehnung erhält dieser Tensor Kugelgestalt. Eine weitere 
Drehung bringt schließlich auch den zweiten Tensor auf Hauptachsen. Da 
Drehung und Dehnung aber lineare Transformationen sind, kann man diese 
drei nacheinander durchgeführten Transformationen durch eine lineare 
Transformation ersetzen. 

Wird das neue Koordinatensystem durch einen Wu gekennzeichnet, 
so bestehen die 5: 


4x” 


q = cy q *— og’, 
mit den reellen Transformationskonstanten cy und cy. Nach der Transforma- 
tion bekommt die Gl. (6,2) die Gestalt 
(Aw) + 3 Buy x’) == Re (7,2) 
wobei die Klammern an den Indizes andeuten, daß auf Grund der Haupt- 
achsengestalt der Tensoren keine Summation mehr notwendig ist. Sämtliche 
A, und Bi, sind positiv oder null. Löst man nach den q’* auf und 


transformiert auf die ursprünglichen ungestrichenen Koordinaten zurück, so 
gewinnt man schließlich den gewünschten Zusammenhang 


q* — pri 
mit dem kontravarianten Tensor $*? 


x’ x’ u 


wobei zur Abkürzung die Frequenzen der Verschiebungsrelaxation s,°) 


9) Genauer gesagt handelt es sich hierbei nicht um Frequenzen, sondern um Kreis- 
frequenzen. 
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eingeführt wurden. Insgesamt treten n Relaxationsiboquensen auf, die jedoch 
nicht alle voneinander verschiedene Werte haben müssen. Im übrigen können 
auch die Werte Null bzw. Unendlich auftreten, falls Aj,,, bzw. Bi,» ent- 
sprechend dem positiv semidefiniten Verhalten der quadratischen Form (5,4) 
bzw. (5,5) null wird. Dabei wird jedoch der Fall, daß Aj,,, und By, für 
ein bestimmtes x’ gleichzeitig verschwinden, ausgeschlossen, da er lediglich 
bedeutet, daß überflüssige Zustandsvariable vorhanden sind, die in jedem 
Fall eliminiert werden können. 

Wie man der Produktform der Gl. (7,4) entnehmen kann, ist der kontra- 
variante Tensor ®**, der mit Ausnahme der Punkte z = — 3, eine analytische 
Funktion von z ist, wieder symmetrisch. 


(7,8) 


Er steht mit dem Tensor W,, der Ausgangsbeziehung (6.2) in einem einfachen 
Zusammenhang. 
Dabei ist 6; der Kronecker-Tensor im V-Raum. a re 


Betrachtet man die Terme mit s, = 0 und s, =oo gesondert, so kann 
man den Tensor ®** in drei Glieder zerlegen. 
Bx* (z 2) = byt + + (2). 9) 
Hierbei sind b*? wo bx? Konstanten, und der verkürzte Tensor gr ent- 
halt nur noch Terme mit endlichen Relaxationsfrequenzen s,,, deren charak- 
teristischer Verlauf über @ fiir 2=iw mit entsprechender Aufspaltung in 
Real- und Imaginärteil in Abb. 1 dar- 


Im ist. Dagegen ergibt das erste 

| PR Glied von Gl. (7,9) einen Anteil —1/w 
We. zum Imaginärteil und das zweite Glied 

ba -m einen konstanten Beitrag zum Realteil. 
Interessiert man sich nur für die 

: äußeren Verschiebungen r*, so kann man 

> x die Ergebnisse als Tensorgleichung im 


schiebungen (dem AV-Raum) deuten, der 
Abb.1. Frequenzabhängigkeit von zum Unterschied vom V-Raum durch 
Real- und Imaginärteil eines dis- k 
kreten Terms der Komponenten des lateinische Indizes gekennzeichnet werden 
verkürzten Modultensors der Ver- soll. Dabei ergibt sich die Beziehung 


schiebungsrelaxation — (7,10) 


wobei die lateinischen Indizes die Werte 1 bis a annehmen können. Durch 
Vergleich der Formel (7,10) mit (7,3) erkennt man, daß der Modultensor 
der Verschiebungsrelaxation R*' mit den entsprechenden Komponenten von 
B** für x,4 <a identisch ist. Die oben angeführten Symmetrie- und Zer- 
legungseigenschaften von ®** sind deshalb auch für R*! zutreffend. Schließ- 
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‚ch lich erkennt man an der Produktgestalt der Gl. (7,4), daß die quadratische 
1en Form R 

nt- R (z) = K, (7,11) 
4) für reelle X, immer durch eine Bezieh 

| mit nicht- -negativen P, gebildet wird, wobei zur Charakterisierung dieser 
ch Eigenschaft ein neuer Begriff eingeführt werden soll: Die quadratische Form 
she 


(7,11) ist, genau wie jedes Diagonalelement des Modultensors der Verschie- 
bungsrelaxation R**, eine spektralpositive Funktion. Dabei werde noch er- 
‚6) wähnt, daß für z = i w der Realteil einer spektralpositiven Funktion niemals 
negative Werte bzw. der Imaginärteil niemals positive Werte annehmen kann, 
wobei auf eine präzisere Formulierung dieses Tatbestandes an dieser Stelle 
jedoch verzichtet werden soll. 


8. Die äußeren Kräfte als Funktion der äußeren Verschiebungen Te 


Während in der Ausgangsbeziehung (6,2) die äußeren Kräfte als Funktion 
'‚8) sämtlicher Verschiebungen auftreten, sollen nunmehr die äußeren Kräfte 
nur als Funktion der äußeren Verschiebungen dargestellt werden. Dabei 
wird von einem — durch einen Querstrich gekennzeichneten — neuen Koor- 
dinatensystem ausgegangen, daß die inneren Verschiebungen 4* bereits in 
Hauptachsenform enthält. Die ersten a Gleichungen der Beziehung (6,2) 
9) lauten dann wie bisher 


nn 


+ Wea = Ke, (8,1) 
nt- wobei lediglich an Stelle der Verschiebungen q* unter Beriicksichtigung der 
ık- Gl. (2,13) in ausführlicher Schreibweise die äußeren Verschiebungen r* und 
an die inneren Verschiebungen 9* eingeführt wurden. Die weiteren i Gleichungen 
i haben auf Grund der oben gemachten Voraussetzung die einfache Form 
ste 
Wat? + = 0. (8,2) 


Setzt man die Beziehung (8,2) in die Gl. (8,1) ein, so kann man die inneren 
eil. Verschiebungen * eliminieren und erhält nach Übergang auf lateinische 


die Indizes die gewünschte Tensorgleichung im AV-Raum 

an 

im = Kirt! 
er- mit dem kovarianten Tensor Re: 


Wie man Gl. (8,4) entnehmen kann, ist der Modultensor der Kraftrelaxation 
10) R,; wieder symmetrisch. 

ch Rır- (8,5) 
oud Führt man die reellen Tensorkomponenten A,, und B,, ein, so ergibt sich 
schließlich eine zur Gl. 9) analoge Aufspaltung 


B- Rz: = Tor: + Tot 6) 
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mit den Konstanten 1,7; Tooxı 

Tort = Bei — = (8,8) 

»=4+1 Bun if 

und dem verkürzten Modultensor der Kraftrelaxation t,, baa = 
kt ’ ( ’ ) 

Bun 8, + 2 

wobei zur Abkürzung die Frequenzen der Kraftrelaxation 3,°) Be 

3 


- 
8,10 
Bun 
sowie analog zur Gl. (7,4) die Transformationskonstanten Ch a shied 
4 
eingeführt wurden. Der verkürzte Tensor ‘tn der nur noch endliche von null 
verschiedene Relaxationsfrequenzen enthält, besteht maximal aus i Termen, 
es wobei jedoch wiederum verschiedene der 
Ra. Im cs “8, den gleichen Wert annehmen können. 
t ey Dazu treten die Glieder mit den Kon- 
stanten (8,7) und (8,8),die den Relaxations- 
frequenzen null und unendlich entsprechen. 
Rs; Der charakteristische Verlauf des Real- 
und Imaginärteils eines einzelnen Terms 
des verkürzten Modultensors der Kraft- 
| 0 0 relaxation t,, über » für z=iq zeigt 
—~ %, Abb. 2. Dagegen liefert das erste Glied 
Abb.2. Frequenzabhängigkeit von von Gl. (8,6) einen konstanten Beitrag 
Real- und Imaginärteil eines dis- zum Realteil und das zweite Glied einen 


kreten Terms der Komponenten des |: - 
verkürzten Modultensors der Kraft- linear anwachsenden Anteil zum Ima- 


ginärteil. Bei der quadratischen Form R 
R (z) = (2) a! (8,12) 


‘ist für reelle x* die entsprechende durch z dividierte Größe R/z wie jedes analog 


gebildete Diagonalelement des Modultensors der Kraftrelaxation R,;/z 
wiederum eine spektralpositive Funktion. Für den verkürzten Modultensor 
T,; ergibt sich diese Tatsache aus der Produktform der Gl. (8,9), für die abge- 
spaltenen Grenzglieder folgt sie unmittelbar aus den speziellen Voraussetzungen 
für die quadratischen Formen U, und D. 

Zwischen dem Modultensor der Kraftrelaxation R;; — der wieder mit Aus- 
nahme der Punkte z= —s, eine analytische Funktion von z darstellt — 
und dem im letzten Kapitel definierten Modultensor der Verschiebungs- 
relaxation R*! besteht folgender einfache Zusammenhang: 

Dabei ist or der Kronecker-Tensor im AV-Raum. 
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9. Kontinuierliche Spektren 


Bei den bisherigen Ableitungen galt als Voraussetzung, daß nur eine end- 
liche Anzahl ö innerer Verschiebungen y* auftreten. Führt man in den Be- 
ziehungen (7,4) bzw. (8,7) bis (8,9) einen entsprechenden Grenzübergang 
durch, so lassen sich die Ergebnisse jedoch auch auf den Fall erweitern, daß 
die Zahl der inneren Verschiebungen gegen unendlich geht. Dabei können 
kontinuierliche Relaxationsspektren, bzw. Spektren mit kontinuierlichen 
Anteilen auftreten. 


Eine übersichtliche Darstellung für sämtliche möglichen Fälle gewinnt 
man, wenn man ganz allgemein an Stelle wu Summen ee Inte- 


sich dam 
= 


(9,1) 


ehe der Modultensor der Kraftrelaxation Rez durch eine entsprechende 
Formel 


gegeben wird, wobei der Faktor s in Gl. (9,1) aus historischen Gründen einge- 
fügt wurde. Die dabei definierten Spekraltensoren P*', P,, sind symmetrisch 


pki — P!k P,, = Py, (9, 3) 


und lassen sich formal analog den Gin. (7,4) und (8,9) in Produktform darstellen 


Pr 8 =: = 9,5 


wobei eine auftretende Entartung allerdings noch besonders berücksichtigt 
werden muß. Ganz allgemein kann man jedoch feststellen, daß die quadra- 


tischen Formen P, P 


(9,7) 


für reelle K,, x* und reelle s> 0 niemals negativ werden. Speziell folgt 
daraus, daß die Diagonalelemente der Spektraltensoren P**, P,, ebenfalls 
niemals negative Werte annehmen. 

Zu diskreten Spektren gelangt man bei der Integralform (9,1), (9,2) unter 


Verwendung der Diracschen Deltafunktion. Dabei sind für die Darstellung 
der Relaxationsfrequenzen Null und Unendlich Grenzübergänge notwendig, 
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wobei man von einer endlichen sehr kleinen bzw. sehr groBen Relaxations- 
frequenz ausgeht. 


Eine Zusammenstellung der wichtigsten Ergebnisse der bisherigen Rech- 
nungen, die das Relaxationsverhalten bei periodischer Erregung betrafen, 
zeigt Tab. 1. Sie enthält außerdem in der letzten Zeile das unter Anwendung 


der (zweiseitigen) Laplace-Transformation™) gewonnene Superpositions- 
gesetz, das vom Bezugszustand ausgehend den Verlauf der Verschiebungen 
bzw. der Kräfte beschreibt, falls der zeitliche Verlauf der Kräfte bzw. der Ver- 
schiebungen vorgegeben ist. 


10. Das Relaxationsverhalten bei sprungförmiger Erregung 2 


Eine ebenso große theoretische und praktische Bedeutung wie das Relaxa- 
tionsverhalten bei periodischer Erregung besitzen die entsprechenden Eigen- 
schaften der Materie bei sprungförmiger Erregung. Die zugehörige Lösung 
läßt sich analog dem periodischen Fall in elementarer Weise entwickeln. 
Man kann aber andererseits an die Rechnungen für periodische Erregung an- 
knüpfen und im Superpositionsgesetz der Tab. 1 eine sprungförmige Erregung 
einsetzen. Die entsprechenden Ergebnisse dieser Transformation sind in 


1) G.Doetsch, Theorie und Anwendung der Laplace-Transformation. Berlin 


1937. 


EN 


Tabelle 1 
Die Grundgleichungen der Relaxation für periodische Erregung 
Verschiebungs-Relaxation Kraft-Relaxation 
Erregung K,, (t) = Re a* (t) = Re {r* e*} 
kn — k zt K zt 
Wirkung x ( ) Re {r e } k (t) ~ {®, e } 
Rg, R, = Ker 
R [ Ber) == [ 
0 0 
Modul = 
, = 
1 i+oo+.a 1 +to+a 
a* (t) = fe z) | R z) (z) dj 
Superpositions- (t) ( ( x (t) ( 
gesetz 
+00 +00 
R,.(z) = [ re= [ 
-- 00 
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Tab. 2 zusammengestellt. Dabei wurde zur Charakterisierung der Sprung- 
erregung der Einheitssprung e, (t) 


N t>0 
verwendet. Außerdem ist zu beachten, daß die in Tab. 2 und den weiteren 
Tabellen auftretende Größe ¢ als Grenzwert e gegen null aufzufassen ist. = 
Tabelle 2 
Die Grundgleichungen der Relaxation fiir sprungförmige Erregung ine 
| Verschiebungs-Relaxation Kraft-Relaxation 
Erregung = Rye (t) (t) = 
Wirkung a* (t) = R** (t) K, K, (t) = () 


— 
=f Pe [1 (t) =f Py, (8) ds t>0 
0 0 


tensor 


= — 1) Ry» (t) dt = OF 
Superpositions- t | tre _ 
Ky) = f Ry dr 
-00 -00 


Führen sämtliche äußeren Kräfte X, im Zeitpunkt t= 0 einen Sprung 
der Höhe K, durch, so wird der zeitliche Verlauf der äußeren Verschiebungen 


x* durch das Produkt von K „mit einem zeitabhängigen symmetrischen Tensor 
zweiter Stufe R*! beschrieben, der als Fundamentaltensor der Verschiebungs- 
relaxation bezeichnet werden soll. Werden andererseits die äußeren Ver- 
schiebungen vorgegeben, so tritt ganz entsprechend der zeitabhängige Funda- 


mentaltensor der Kraftrelaxation R,, in Erscheinung. Analog den Modul- 


Abb. 3. Zeitabhängigkeit Abb. 4. Zeitabhängigkeit 
eines diskreten Terms der eines diskreten Terms der 
Komponenten des ver- Komponenten des ver- 
kiirzten Fundamental- kürzten Fundamentalten- 
tensors der Verschie- sors der Kraftrelaxation 
bungsrelaxation 


tensoren spalten die Fundamentaltensoren ebenfalls in drei Glieder auf, 
wobei der charakteristische Verlauf eines diskreten Terms der verkürzten 
Fundamentaltensoren in der Abb. 3 bzw. 4 dargestellt ist. Neben den Re- 
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laxationsfrequenzen s führt man bei der sprungförmigen je zweck- 
mäßig die dazu reziproken Relaxationszeiten r 

(10,2) 
ein, auf die sich z. B. die Indizierung der ersten beiden Glieder der Funda- 
mentaltensoren bezieht. 

Außerdem enthält Tab. 2 die Integro-Differentialgleichung, welche die 
beiden Fundamentaltensoren verkniipft. SchlieBlich wird in der letzten Zeile 
das Superpositionsgesetz für sprungförmige Erregung angegeben, das die 
fundamentale Bedeutung der eingeführten zeitabhängigen Tensoren zum Aus- 


druck bringt. Die quadratischen Formen R, R ef 
R (t) = Rk(t) K, K, 
R(t) = (10,4) 


werden für reelle K,, 2* und reelle t > 0 ebenso wie die Diagonalelemente der 


Fundamentaltensoren R**, R,, niemals negativ. Genauer gesagt verlaufen 
hierbei die entsprechend gebildeten quadratischen Formen der verkiirzten 
Fundamentaltensoren 7, 7 sowie die Diagonalelemente der verkiirzten Fun- 
damentaltensoren 7**, 7,,, totalmonoton?). 

Der mathematische Zusammenhang zwischen den Fundamentaltensoren 
und den zugeordneten Modultensoren zeigt Tab. 3. Zwischen den Konstanten 
der Grenzglieder bestehen dabei spezicll folgende einfache Beziehungen: 


Pin Tabelle 3 even} 4 


Zusammenhang zwischen Modultensor und Fundamentaltensor 


Verschiebungs-Relaxation Kraft-Relaxation 
© co 

(z) =z | R*' (t) dt Rer (2) = 2 f Ry, (the dt 

0 -e 
1 +ioo+a 1 — @) 
BY) = | de Ru) = f de 

-ioo+.a -too+e 


11. Die Berechnung der Spektraltensoren 


Während der Modultensor und der Fundamentaltensor über spezielle 
Erregungsformen definiert sind, ermöglicht der Spektraltensor eine von der 
Erregungsform unabhängige Charakterisierung des zugrunde liegenden rela- 
xierenden Mediums. Es besteht deshalb Interesse, den Spektraltensor bei 
vorgegebenem Modultensor bzw. Fundamentaltensor berechnen zu können. 

In Tab. 4 sind die entsprechenden Ergebnisse zusammengestellt, wobei 
in der ersten Zeile zunächst die Gliederung der Spektraltensoren angegeben 
wird. Danach folgen die Formeln für vorgegebenen Modultensor, die nach den 
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eck- Tabelle 4 
Berechnung des Spektraltensors bei vorgegebenem Modultensor bzw. Fundamentaltensor 
0,2) Verschiebungs-Relaxation | Kraft-Relaxation 
6 (8s —8 
(8) = lim + Py, (8) = Ooxz lim 6 (8 — %) + 
die & 8 —>0 
eile +05 lim 6(s — %) + 0*' (s) + lim (s — 8) + (8) 
die 00 >00 
us- 1 1 
P¥* (s) = + — Im ie} (8) = + — Im {Re (— 8 +ie)} 
| (a) = — (e F is} Ger (8) = +2, Im (e 
), 4) 2 3 2 
= — Re {rij (e tie} (8) = — Re {Eyer + 
en i 
en P* (s) (ty dt P, (8) TR (t) dt 
| = ı(8) = 
INn- 2218 > 228 
2 2 pki D (a) do k 
en (0) Pin (0) | f 
0 
Regeln der Stieltjes-Transformation™) gewonnen werden. Dabei kann 
man den verkiirzten Spektraltensor unter Verwendung der Kramers- 
.- Kronig-Beziehung!*) bereits berechnen, wenn nur Real- oder Imaginärteil 
des verkürzten Modultensors r*! (i w) bzw. f,, (i @) bekannt sind, die durch 
ws die eingeklammerten Indizes r und i gekennzeichnet werden!2). Anschließend 
sind die Beziehungen für vorgegebenen Fundamentaltensor dargestellt, die 
sich nach den Regeln der Laplace-Transformation™) ergeben. Dabei ist 
zu beachten, daß die mit einem Stern versehenen Fundamentaltensoren zum 
Unterschied von den entsprechenden bisher betrachteten Größen ohne Stern 
für £< 0 nicht verschwinden, sondern in diesem Bereich die analytische 
Fortsetzung dieser Größen für t> 0 darstellen. Für die Berechnung der Kon- 
stanten der Grenzglieder gelten folgende einfache Zusammenhänge: 
of = = rf! = = Took 
e = t= 7h Door ı = Toort = 
r SchlieBlich wird noch die Querbeziehung zwischen den beiden Spektral- 
- tensoren angegeben. 
1) B. Gross, J. appl. Phys. 18, 212 (1947); T. J. Stieltjes, Recherches sur les 
: fractions continues, Ann. de la faculté de Toulouse 8, 1 (1894). 
1 12) H.A. Kramers, Atti Congr. dei Fisici, Como 1927, S. 545; R. de L. Kronig, 
J. opt. Soc. Amer. 12, 547 (1926). 
| 12a) ist der negative Imaginärteil von 
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ws 12. Anwendung auf die mechanische und dielektrische Relaxation 


Die Grundgleichungen der mechanischen und dielektrischen Relaxation für periodisd 
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Bei einer speziellen Anwendung der in den Tab. 1 bis 4 dargestellten 
Grundgleichungen der Relaxation muß man die geometrische Bedeutung der 
Modul-, Fundamental- und Spektraltensoren im gewöhnlichen dreidimen- 
sionalen Raum berücksichtigen. Für die mechanische und dielektrische 
Relaxation sind die entsprechenden Ergebnisse in den Tab. 5 und 6 zusammen- 
gestellt. Dabei wird nur auf das Verhalten bei periodischer Erregung einge- 
gangen, bei dem die Modultensoren in Erscheinung treten. Sämtliche Aus- 
sagen lassen sich jedoch unmittelbar auf die Fundamental- und Spektral- 
tensoren übertragen. 


Tabelle 5 


Erregung 


Verschiebungs-Relaxation Kraft-Relaxation 
xa _ We In 
Mechanische M+M t=ME HM, 
j Dielektrische Ss = Dt + » e, t =. Ds + D,d 
| Relaxation = D*t + De, e, = 9,3 + D,,b 
3=Rt+ R%e, + R78 t= Rs +R, d* + R,, 
Mechanisch- 
dielektrische = + + | & = R, R,; 
| 2 t +3 eu # 7 = Run Ryan? 
1 Se 
Tabelle 6 
ae Die Grundgleichungen der mechanischen und dielektrischen Relaxation isotroper Medi 
für periodische Erregung 
Verschiebungs-Relaxation Kraft-Relaxation 
| Mt + Ny" t= MS + 
+ (Br + (B 9,9, + Jur)” 
Dielektrische $= Dt t = Ds 
Relaxation A a 
A. Mechanische Relaxation 


_ dichteabweichung vom Bezugszustand © — Sp = 3 und die Komponenten des 
_ Deformationstensors d** auf, während als äußere Kräfte $, die Temperatur- 


Im mechanischen Fall treten als äußere Verschiebungen r* die Entropie- 


abweichung T— Q,=t und die Komponenten des Spannungstensors 3,; 
wirksam sind. Dabei wurden sofort die bei periodischer Erregung maß- 
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gebenden — durch Frakturbuchstaben gekennzeichneten — komplexen 
Amplituden angegeben. 

Bei der Relaxation eines beliebigen anisotropen Mediums wird der im 
AV-Raum definierte Modultensor der Verschiebungsrelaxation R*! im drei- 
dimensionalen Raum durch einen Skalar M, einen Tensor zweiter Stufe M*? 
und einen Tensor vierter Stufe M***” charakterisiert. Ganz analog gliedert 
sich der Modultensor der Kraftrelaxation R,., in die entsprechenden Größen 
M, Maras 

Davon ist M bzw. M mit der (frequenzabhängigen) spezifischen Wärme 
bei konstanten Spannungen c, bzw. bei konstanten Deformationen c, ver- 
knüpft 

— 90% m — Ton 

M = T, M nee ); (12,1) 
wobei neben der Dichte o, die Temperatur 7, des Bezugszustandes in die 
Rechnung eingeht. Bei der Ableitung dieser Formeln muß man berück- 
sichtigen, daß im Rahmen einer linearen Näherung von Beziehungen der 
Art (2,10), (2,28) die Zunahme der Entropiedichte dS nur durch die zuge- 
führte Warmedichte dQ bestimmt wird 


T, dS = dQ, (12,2) 


während die spezifische Wärme c andererseits ganz allgemein durch die 
differentielle Verkopplung von Wärmedichte und Temperatur definiert ist. 


dQ = 0, ¢ aT. (12,3) 


Die sechs unabhängigen Komponenten des symmetrischen Tensors M*’ 
bestimmen die durch die Temperatur erzeugten Deformationen bei kon- 
stanten Spannungen, während der Tensor M,; die durch die Entropiedichte 
erzeugten Spannungen bei konstanten Deformationen beschreibt. Anderer- 
seits werden durch den in dreifacher Hinsicht symmetrischen Tensor vierter 
Stufe Mir” bzw. Myru,, der 21 unabhängige Komponenten besitzt, die 
isothermen bzw. adiabatischen Spannungs-Deformations-Beziehungen erfaßt. 
Dabei ist im letzteren Fall zu beachten, daß im Sinne der linearen Näherung 
der Gl. (12,2) kein Unterschied zwischen adiabatischen und isentropen Vor- 
gängen gemacht wird. 

Ist das relaxierende Medium isotrop, so müssen die Tensoren M*?, M,; 
und M**“", M2,» bis auf einen skalaren Faktor durch den Metrik-Tensor 
g** dargestellt werden. Da es hierbei bei einem Tensor zweiter Stufe nur 
eine Möglichkeit der Zuordnung gibt, bei einem dreifach symmetrischen 
Tensor vierter Stufe aber zwei, bedeutet dies, daß das Relaxationsverhalten 
eines isotropen Mediums durch vier skalare Modulfunktionen der Verschie- 
bungsrelaxation M, N, PB, Q bzw. durch vier skalare Modulfunktionen der 
Kraftrelaxation M, N, $8, OQ charakterisiert wird. Daneben kann man aber 
auch von diesen acht Funktionen vier beliebige auswählen und schließlich 


13) Diskutiert man die Beziehung (12,1b) speziell für & = 0, so erkennt man, daß 
die im Anschluß an Gl. (2,24) aufgestellte Behauptung, daß A,, immer verschieden von 
null und positiv sein muß, zu Recht besteht. 
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noch dav on abgeleitete Funktionen verwenden, wie z. B. die Ilse Kom-] Ve 


bination, die nur spektralpositive Funktionen enthält). Vo 
Dabei sind c,, c, die spezifischen Wärmen bei konstantem Druck bzw. bei “te 
konstantem Volumen, $,, bezeichnet den adiabatischen Kompressionsmodul oe 


und & den Schermodul, der fiir adiabatische und isotherme Vorgänge denselben 
Wert besitzt. 

Schließlich werde darauf hingewiesen, daß man im allgemeinen bei der 
mechanischen Relaxation für die Begriffe Verschiebungsrelaxation und Kraft- 
relaxation die Bezeichnungen Retardation und Relaxation verwendet. Man 
muß dann zwischen Relaxation im weiteren und im engeren Sinne unter- 
scheiden. La 


B. Dielektrische Relaxation 


Im dielektrischen Fall treten als äußere Verschiebungen r* neben der 
bereits im mechanischen Fall eingeführten Entropiedichteabweichung 8 
die Komponenten des dielektrischen Verschiebungsvektors d* auf, während 
als äußere Kräfte §,, neben der Temperaturabweichung t die Komponenten 
des elektrischen Feldstärkevektors e, maßgebend sind. 

Bei anisotropen Medien wird der Modultensor der Verschiebungsrelaxation 
im dreidimensionalen Raum durch einen Skalar D, einen Vektor D* mit drei 
Komponenten und einen symmetrischen Tensor ®** mit sechs unabhängigen 


‘dl Komponenten vertreten. Ganz entsprechend ordnet man dem Modultensor 
7 der Kraftrelaxation die Größen D, D,, D,, zu. 
ao Davon ist D bzw. D durch die spezifische Wärme bei konstanter Feld- 
stärke c, bzw. bei konstanter dielektrischer Verschiebung c,, bestimmt. 
T 
D = 12,5 
T, 00 Cav ( 


Der Vektor ®* beschreibt die durch die Temperatur erzeugte dielektrische 


Verschiebung bei konstanter elektrischer Feldstärke und der Vektor ®, die 
entsprechend durch die Entropiedichte erzeugte Feldstärke bei konstanter 
_ dielektrischer Verschiebung. Andererseits wird durch den symmetrischen 
Tensor D** bzw. D,, die isotherme bzw. adiabatische Feldstärke-Verschie- 
bungs-Beziehung erfaBt. 

Bei isotropen Körpern verschwinden ®*. und ®,, und das Gleichungs- 
system zerfällt deshalb in je eine Beziehung zwischen den zwei thermischen 
und elektrischen Variablen, wobei die Verkopplung jeweils durch-einen Skalar 

erfolgt. Das dielektrische Verhalten isotroper Medien ist also unabhängig 
von gleichzeitig stattfindenden thermischen Zustandsänderungen. 
: Schließlich sei darauf hingewiesen, daß man auf dielektrischem Gebiet 
in allgemeinen oft nur die Verschiebungsrelaxation betrachtet und dafür 
abgekürzt die Bezeichnung Relaxation verwendet. Dies ist besonders beim 


14) Genauer gesagt sind hierbei neben c,, c, nur 8,a/2 und @/z spektralpositiv. 
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Vergleich der dielektrischen Relaxation mit entsprechenden mechanischen q 
Vorgängen zu beachten, da man im letzteren Falle — wie oben angegeben ee 
wurde — darunter etwas anderes versteht. oy 
| 
C. Mechanisch-dielekrische Relaxation gust 
Eine gekoppelte mechanisch-dielektrische Relaxation kann nur in aniso- 
tropen Medien auftreten. Dabei werden die Modultensoren im dreidimensiona- 
len Raum durch je einen Skalar, einen Vektor, zwei Tensoren zweiter Stufe, 
einen Tensor dritter Stufe und einen Tensor vierter Stufe charakterisiert, 
wobei insgesamt 55 unabhängige Komponenten auftreten. 
Dagegen zerfällt beim isotropen Körper das Gleichungssystem in einen 
thermisch-mechanischen und einen elektrischen Teil. 


Berlin-Adlershof, Deutsche Akademie der Wissenschaften zu Berlin. 
Laboratorium für Kunststoffe. 


Bei der Redaktion eingegangen am 23. November 1957. 
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s Milnesche Problem bei anisotroper Streuung 
(Anwendung auf thermische fans 1) 


Es wird das ber ühmte te Milnesche Problem von neuem aufgegriffen und 
gezeigt, daß sich das ,,Standard-Problem“ (s. Text) ohne besondere analyti- 
sche Hilfsmittel (wie Funktionentheorie) in sehr guter Näherung lösen läßt, 
aber nur, wenn man von vornherein Anisotropie des elementaren Streuaktes 
annimmt. Die Dichteverteilung am Rande des Mediums wird nach Exponential- 
integralen entwickelt. Der ,,isotrope‘‘ Fall läßt sich durch Extrapolation dar- 
aus gewinnen. Nach dieser Methode ergibt sich auch ein trivialer Beweis 
der bekannten ,,Hopf-Bronstein- Beziehung“. Es wird begründet, weswegen 
die Behandlung bei isotroper Streuung so außerordentlich kompliziert ist. — 
Die Abweichungen von der Isotropie, soweit sie nicht schon in der Trans- 
portweglänge einbegriffen sind, müssen sich experimentell merklich äußern. 


I. 
a). 

a Problemstellung ist dieselbe wie bei Placzek?) (,,Standardproblem‘), 
nur mit dem Unterschied, daß das dort vorausgesetzte isotrope Streugesetz 
hier durch ein allgemeineres ersetzt wird. Die Streufunktion enthält bei uns 
noch ein Glied, das quadratisch vom Kosinus des Streuwinkels 6, um den das 
thermische Neutron aus seiner ursprünglichen Bewegungsrichtung abgelenkt 
wird, abhängt. Ein Parameter kennzeichnet den Grad dieser Art von Aniso- 
tropie. In je einem Spezialfall ergibt sich entweder das isotrope Gesetz wie 
bei Placzek bzw. das der Rayleighschen Phasenfunktion entsprechende, 
das von Chandrasekhar) untersucht wurde. Ein Glied in der Streu- 
funktion, das linear vom Kosinus des Winkels 6 abhängt, stellt sich für unsere 
Untersuchung als von untergeordneter Bedeutung heraus, und wird nach- 
träglich leicht durch eine Einführung der Transportweglänge berücksichtigt. 
Wir wollen diesen Punkt nur flüchtig berühren und lassen dann besagtes 
Glied in der Streufunktion fort. Die Feldfunktionen werden dadurch nicht 
beeinflußt. Z. B. ist die Winkelverteilung am Rande vollständig unabhängig 
von diesem linearen Glied. Voraussetzung ist natürlich, daß die Absorption 
im Medium vernachlässigt wird. 


1) Mitteilung aus dem Kernphysikalischen Institut Zeuthen-Miersdorf der D. A.d.W. 
zu Berlin. 

2) G. Placzek u. W. Seidel, Physic. Rev. 72, 550 (1947). 
S. Chandrasekhar, ,,Radiative Transfer“. Oxford 1950. ln 
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Wir wiederholen kurz die Problemstellung. Untersucht werden die Trans- 
porterscheinungen für Neutronenstrahlen, die in einem absorptionsfreien 
Medium nach dem bestimmten Gesetz gestreut werden. Das Medium ist ein- 
seitig, etwa nach rechts hin, unbegrenzt, während sich links als Begrenzungs- 
fläche eine seitlich unendlich ausgedehnte ebene Oberfläche befindet, auf die 
von der dem Medium abgewandten Seite, also von links, keine Neutronen auf- 
treffen sollen. Dagegen bedindet sich auf der anderen Seite in einiger Ent- 
fernung von der Begrenzungsebene (also einige freie Weglängen weit) eine ther- 
mische Neutronenquelle. Man denke dabei an die Verhältnisse in einem Pile. 
Die Neutronenquelle dient dazu, einen konstanten Nettostrom thermischer 
Neutronen in Richtung der Oberfläche aufrechtzuerhalten. Die gemachte 
Voraussetzung, daß nämlich die Absorption der Neutronen in dem Medium 
zu vernachlässigen ist, läßt sich rechtfertigen, indem man das Neutronenfeld 
nur in Abständen vom Rande betrachten will, die klein gegen die Diffusions- 
länge sind. In diesem interessierenden Feldraum sind dann auch räumliche 
thermische Neutronenquellen zu vernachlässigen. 

Die wesentlichsten Ergebnisse unserer Untersuchung sind folgende: Die 
Neutronendichte an der Oberfläche ist nicht mehr genau y3 pro Einheits- 
Nettostrom, wie die bekannte Hopf-Bronstein-Beziehung das für isotrope 
Streuung fordert. Genauer ergibt sich, daß die Dichte größer als dieser Wert 
ist, wenn cos? 0 > !/,, dagegen kleiner, wenn cos* 0 < !/,. Die Beträge dieser 
Abweichungen lassen sich leicht angeben, wenn eine verallgemeinerte Hopf- 
Bronsteinsche Beziehung benutzt wird. Aus dieser neuen Beziehung kann 
man dann weiter durch einfache Abschätzung entnehmen, daß sich die soge- 
nannte Extrapolationslänge um weniger als 1%, von der für isotrope Streuung 
abgeleiteten unterscheidet. (Voraussetzung sind natürlich gleiche Transport- 
weglängen). Der Dichteverlauf in unmittelbarer Nähe der Oberfläche hängt 
dagegen wesentlich von dem Unterschied zwischen cos? und dem ,,iso- 
tropen‘“ Wert !/, ab. 


Die Aufstellung der Fundamentalgleichung, die den Diffusionsvorgang 
beschreibt, ist etwas verwickelt). Wir haben zwei Grenzfälle: 
lw=0, w=!]|; 


2w,=1, w,=0. Vokuum Medium 


Uns wird später nur 2. interessieren. w, 

und w, sind Streu- und Absorptinswahr- 

scheinlichkeit. 
0 

Grenzfall 1 


Ein Beispiel. — Auf die Begrenzungs- 

ebene z = 0 (Abb. 1) fallt ein vorgegebenes . 1. 
; 4 A strahl tritt in das Medium ein. # ist 

Spektrum ® (0, u) von links ein. m ist der Winkel, den der Strahl mit der 
die Geschwindigkeitskomponente, also positiven z-Achse einschließt. Die Ge- 
4 = cos #, wenn der Betrag von vgleich 1. schwindigkeit v ist Eins, so daß u = 
Es werden keine Strahlen zurückgestreut, °0s® = Komponente der Geschwin- 
digkeit in Richtung z Fi 
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‘) W. Verde u. G. C. Wick, Physic. Rev. 71, 852 (1944). a re 
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die Wahrscheinlichkeit, einen Strahl innerhalb du hinter der Ebene z anzu- 
treffen. z ist in Einheiten der freien Weglänge zu messen. 


Die Fundamentalgleichung (Absorptionsgleichung) lautet 


Wenn w, + 0 ist, so sind die Verhältnisse natürlich wesentlich komplizierter. 
Wir haben (Boltzmann -Gleichung): 


op , 
wo a als Operator gedacht ist, der auf den Geschwindigkeitsraum wirkt und 
vorhin einfach Einheitsoperator war. Jetzt trägt a allgemeinen tensoriellen 
Charakter. Es ergibt sich, daß gilt: 


Hier sind 
R=zW#-;, A= 
ww 1 
(mi P, Py du = by) 


die Legendreschen Kugelfunktionen, und die Eigenwerte bestimmen sich zu 


b,=w, + (1 _ P, (cos 6), 
wo cos 6 der Kosinus des elementaren Ablenkwinkels ist und der Mittelungs- 
strich sich auf das elementare Streugesetz bezieht. Die Formeln ergeben sich 
aus einfachen Wahrscheinlichkeitsüberlegungen (s. etwa Wick®)). 


Grenzfall 2 


| © Hier hat man nach leichter Umformung 

Dens un. t+ Lo, Py, 

I 


ae = [een Pıdu 


und 


IH 

= P, (cos 6). 
Die offensichtliche Eigenschaft des „Nettostromes“, nämlich durch das Me- 
dium hindurch konstant zu sein, erlaubt eine wesentlich allgemeinere Formu- 
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lierung der Randbedingung bei z = 0 als in 1. Wir setzen diges et ; a 
4 


und können verlangen Aa 
(0, fir 
Dadurch ist das Problem eindeutig bestimmt. Die Randbedingung besagt 
(Abb. 1), daß von links keine Strahlen in das Medium eintreten. 

Das Problem in dieser Fassung können wir als Verallgemeinerung des 
Placzekschen ,,Standardproblems“ auffassen. Der Lösung wenden wir uns 
jetzt zu. 

Die Aufgabe der Lösung des Problems besteht im wesentlichen in dem Auf- 
finden der ,,Dichte“ y,, wo nach dem 


ZU - 


ind Alle anderen Fragen lassen sich dann leicht nach Lésung dieser Aufgabe beant- Ta 
len worten, etwa die nach der Auffindung der Winkelverteilung am Rande (z = 0). 
Was den asymptotischen Verlauf der Dichte fiir groBe z betrifft (in Zukunft 
denken wir uns immer z in Einheiten der Streuweglänge gemessen), so wird 
schon der gewöhnlichen Diffusionstheorie zufolge die Dichte linear mit z 
anwachsen. Dies geschieht aber natürlich nur so lange, bis die hier vernach- 
lässigte Absorption und eventuelle Quellen das verhindern. Man bekommt 


Y~ (1 — cos 6) z. 
Um auch den Verlauf der Dichte für kleinere z (einschließlich 2 = 0) zu be- 
schreiben, fügen wir (nach dem Vorgang von Milne) eine Funktion q von z 
hinzu. Da sich aber unsere Theorie auf ein allgemeines elementares Streu- 
zu gesetz bezieht, wollen wir das berücksichtigen, indem wir qg noch abhängen 


lassen von den Koeffizienten w,, wobei sich aber später herausstellt, daß der 
zweite praktisch ausreichend ist. Wir führen an Stelle von w, eine „charak- 


IS- teristische Zahl“ x ein, die einfach mit w, zusammenhängt und gegeben sein 
ch soll durch: 
(I, 1) 
Wir haben dann Sara. 


(1 — cos 6) 2+q (z, 2). 


Der Fall x= 5 ist von verschiedenen Autoren näherungsweise behandelt 

worden (,,isotroper Fall“). Eine kritische Zusammenstellung der Ergebnisse 

findet man in dem Buch von Kourganoff5). Wir greifen heraus. 
Eddington setzt 


+ E, (2) = | (Exponentialintegrale). 
1 
e- 
*) s. V.Kourganoff, ,,Basic Methods in Transfer Oxford 1952. 
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Das ergibt gute Resultate innerhalb des Intervalls von z = 0 bis un ist aber 
an den Grenzen ungenau. 


Placzek nimmt als Ausgangslösung für ein Iterationsverfahren : 


1 —az ~ 
a= %)—( (x 3). 


Diese Näherung ist so gewählt, daß sich an den Grenzen die genauen Werte 
ergeben. Placzek findet 
2, = 0,7104,4509. (1, 2) 


C. Mark®) hat auf Grund einer exakten Theorie von Placzek die 
genauen Zahlenwerte für die Funktion q für einige Werte von 2 berechnet, 
von denen wir später gelegentlich Gebrauch machen werden. 

Der Fall x = 6 (Rayleighsches Streugesetz) ist von Chandrasekhar 
behandelt worden. Er findet 

% = 0,71139. 


In der gegenwärtigen Arbeit wird nun für alle möglichen ganzen Zahlen 
für x das gq berechnet (x = 3 bis 9), und zwar, indem gq (ähnlich wie bei 
Eddington) durch eine Summe von Exponentialintegralen dargestellt wird. 
Die Entwicklungskoeffizienten in dieser Darstellung werden exakt genug 
nach einem einfachen Verfahren bestimmt, das allerdings in dem erwähnten 


Fall x = 5 versagt (s. Abschnitt V Schluß). wb 


Aus der Boltzmann _-Gleichung 
D+19= II, 1 

und 
21 +1 d 

P, (cos 6); D=— (Operator), (II, 3) 
kénnen wir eine Reihe von Gleichungen ableiten, die wir als ,,Transport- 
gleichungen“ oder auch „Fokker-Plancksche Gleichungen‘ bezeichnen 
werden. 


Wir setzen 
1 
„()= Seem duu"; Yo A = 


Dann lauten die ,,Planckschen Gleichungen‘, wenn py = 4, w, = 0 ist, 
+H = 0 | 


| 2 (II, 4) 
3 Dys + (5 — 2 we) = 0 
Im by a 3 D2 (yı ya + 3 ys) — (5 — 2 w,) (7 — 2 Ws) Y2. = 0 
usw 
Die dritte Plancksche Gleichung z. B. kann man etwa so herleiten: 


6) c. Mark, Physic. Rev. 72, 558 (1947). 
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D. Das Milnesche Problem bei anisotroper 


Wir haben 


und bilden 
f du Po Yo + We (2 P, + Py. 


Die weiteren Glieder auf der rechten Seite sind fortgelassen, da sie wegen der 
Orthogonalitätsrelationen für die Kugelfunktionen nichts — Wir 
bekommen 


Dt i P§ du +3 3 Ya f P3 du. 


Da aber 


Pidu == 


2 1 
3 Y2 Yo 


ergibt sich 
2 3 2 1 2 2 iis FR 
5 Dy + Dy, + 3 + 


Wegen w, = (fehlende Absorption) und Dy, = 0 (erste 
chung an 4) wird die letzte reduziert auf 


Da + zn = 


was offenbar mit der dritten nn Gleichung (II,4) übereinstimmt. 

Die beiden ersten Planckschen Gleichungen (II, 4) lassen sich unmittelbar 
integrieren. Die erste besagt, daß die Divergenz des Stromes verschwindet, 
d.h. in unserem Falle, daß — wie wir bereits vorausgesetzt — der Strom kon- 
stant ist. Wir hatten gesetzt (Stromnormierung): 


=-1. (II, 5) 
Die zweite Gleichung liefert darauf unmittelbar 
Po = (2 + %). (II, 6) 


Die Integrationskonstante z, ist dieselbe Zahl (0,71 . .), die uns bereits begegnet 
ist. Sie ist durch die Randbedingung bei z = 0 bestimmt. Wir kénnen das 
aber auch anders ausdriicken. Wir sagen, bei richtiger Wah! der Konstanten 
2, müssen sich die übrigen Funktionen y, in den Planckschen Gleichungen 
so bestimmen lassen, daß sie dann ein richtiges asymptotisches Verhalten für 
z > cozeigen; sie müssen nämlich für z > coin genügend starkem Maße gegen 
Null streben, — derart, daß die Umkehrung D-! der Operation D einen ein- 
deutigen Sinn bekommt. Es ist 
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(II, 7’) 


Nach diesen Behauptungen (die sich streng begründen ließen), können wir 
_ neue Operatoren einführen durch 


3 1 (5 
usw. 


und haben dadurch die Möglichkeit, alle Funktionen y, mit / > 2 mit Hilfe 
der Planckschen Gleichungen durch y, allein auszudrücken; y, = A, y». 
Da die Dichte mit y, über die Gleichung 


Yo + 2 pe = 3(z + %) (II, 9) 
zusammenhängt, bleibt dann zur Lösung unseres Problems nur noch eine 
Gleichung aufzufinden für y, allein, die y, eindeutig zu berechnen gestattet. 

Zu diesem Zwecke stellen wir eine (verallgemeinerte) Milnesche Integral- 
gleichung für y, auf, die wir folgendermaßen gewinnen. 

Wir gehen von der Boltzmann-Gleichung (II, 1) aus, die sich bezüglich z 
als gewöhnliche Differentialgleichung erster Ordnung auffassen läßt. Sie ist 
linear und hängt von dem: „Parameter yu‘, der Geschwindigkeitskomponente, 
ab. Der Inhomogenteil (die rechte Seite von (II, 1) ist offenbar so gewonnen 
worden, daß die unabhängige Funktion 9 (z, u) über den Parameter inte- 
griert wurde, was auf die Funktionen y, führte. Die formale Lösung dieser 
Differentialgleichung, bei der wir den Inhomogenteil uns als bereits bekannte 
Funktion von 2 vorstellen wollen, können wir so schreiben: 


(II, 8) 


(II, 10) 
Der rechterhand stehende Umkehroperator bekommt einen Sinn, wenn die 
Randbedingung bei z = 0, d.h.: eae 


(0,u)=0, fir 0<p< (II, 11) 
in die Betrachtung einbezogen wird. 

Wir integrieren beide Seiten von (II, 10) über # von —1 bis +1 und nehmen 
die Integration über « in zwei Schritten vor, wobei wir den Punkt u = 0 
des ursprünglichen Integrationsintervalls jeweils in die eine Integrations- 
grenze legen. Dadurch wird erreicht, daß der „kleine Parameterwert“ u = 0 
nicht in Erscheinung tritt. Wir erhalten dann einen Ausdruck für wo, den 
wir so schreiben wollen: 


(IT, 12) 
worin die Operatoren M, folgende Bedeutung haben: at i 
12-2] 


“ [sign (2—¢)¥. (11,19) 


Analoge Formeln ergeben sich dann für die Umkehrungen der höheren Ab- 


D 


Ww 


06 Annalen der Physik. 7. Folge. Bndl. 158° 
0 
ak 
Q, ül 
R 
m 
ei 
J 
ze 
+ 
d 
ate G 
st 
r 
a 
N fe 
s 
RE 
t 
g 
YE 
= 
ul d du P 
My=z | | „PıwWe 
0 0 


Ab- 


2") 


13) 


D. Lyons: Das Milnesche Problem bei anisotroper Streuung 407 


Die Vorzeichenfunktion sign (x) ist definiert durch EEE 
hoist on Altar 


Die Operatoren M, sind hermitisch. 
Man könnte auf ähnliche Weise die Funktionen 9, (2) konstruieren fir 
n> 1. Eine nähere Untersuchung der so entstehenden Gleichungen würde 


aber nichts Neues bringen; sie würde lediglich zeigen, daß unsere Annahme _ 3 ‘ 


über das asymptotische Verhalten der Funktionen y, richtig ist (die einzelnen 
Rechnungen sind zu bekannt, als daß sie hier vorgeführt werden müßten). 


Wir führen jetzt den (im allgemeinen Falle nicht hermitischen) Operator Q > 


ein durch die Gleichung 2 

Q2=2 +23 M,2,, (II, 14) 

2Q,=1-M,- (II, 15) 


Indem wir y, in der Gl. (II, 12) mittels der Beziehung (II, 9) eliminieren und 
%, aus (II, 6) einsetzen, erhalten wir nach leichter Umformung in 


Q 3 (Q—Q)) (2 + %) (11,16) 


eine Milnesche Integralgleichung allgemeinen Typs für die Dichte y,. Im 
„isotropen‘‘ Fall fallen die Operatoren 2, und 2 zusammen, und man über- 
zeugt sich leicht, daß in diesem Spezialfall unsere Integralgleichung äquivalent 
der gewöhnlichen Milneschen Integralgleichung ist. 

In dem Vorhandensein eines Inhomogenteils unterscheidet sich unsere 
Gleichung von der Milneschen. Die in (I) eingeführte Funktion q (unter der 
unwesentlichen Voraussetzung, daß cos 6 = 0 ist) läßt sich jetzt formal so 
schreiben: 


q (2) = - (z + %). (II, 17) 


Die Einklammerung der —1 an dem Symbol für den reziproken Operator Q 


soll andeuten, daß es sich um eine „bedingte‘‘ Umkehrung handelt, was be- 


sagt, daß diese Umkehrung nur möglich ist, wenn der Zahlenwert von 2, ar 


richtig gewählt ist (s. später). In den obigen Formeln soll natürlich 2, (2 + 2)) 
als Funktion von z aufgefaßt werden, auf die die Operatoren wirken. 

Wir haben das Ergebnis: Durch unsere Umformungen ist es gelungen, eine 
formale Lösung unseres Problems hinzuschreiben. Es bleibt noch zu unter- 
suchen, wie die Umkehrung des Operators 2 zu erfolgen hat. Um diese Frage 
beantworten zu können, gehen wir von der allgemeinen Theorie der Integral- 
gleichungen vom Milneschen Typ aus, die im nächsten Abschnitt behandelt 
werden soll. 


II. 


In diesem Abschnitt besprechen wir den mathematischen Formalismus 
der Theorie der Integralgleichungen vom Milneschen Typ. Es sollen die 
Gesichtspunkte hervorgehoben werden, nach denen man zu verfahren hat, 
um eine exakte Theorie zur Lösung unseres speziellen Problems zu gewinnen, 
ohne — wie sonst meist üblich — funktionentheoretische Hilfsmittel heran- 
ziehen zu müssen. Wenn dadurch zwar eine gewisse Schwerfälligkeit der 
Beweisführung (und vielleicht — soweit sie u. durchgeführt wird — 
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auch eine mathematisch nicht ganz einwandfreie) mit in Kauf genommen Es | 


werden muß, so sind wir dafür in der Lage, die physikalische Seite um so mehr weis 
hervorkehren zu können. Es gelingt uns, Beziehungen aufzudecken, die bei 
abstrakterer Formulierung leicht unbemerkt geblieben wären. Diese neu- 


artigen Zusammenhänge zeigen eine gewisse (wie der Verfasser bemerkt zu Es ı 
haben glaubt) Analogie mit einer anderen physikalischen Disziplin, nämlich Rec 
der Elektrostatik. Es erscheint (schon im Hinblick auf die leichte Verständi- zusé 


gung) als zweckmäßig, diese Analogie mit in die Bezeichnungsweisen auf- 
zunehmen. 

Zunächst spalten wir unseren Operator 2 (II, 14) additiv in zwei Teile 
auf, von denen der erste (und wichtigere) hermitisch sein soll. Er werde, 
ohne Verwechslungen befürchten zu müssen, ebenfalls mit 2 bezeichnet. 


Der restliche Teil 2, könnte später mit Hilfe eines Störungsv erfahrens mit- | vee 

beriicksichtigt werden. ri 

25-24 2,. Pur 

2 hat damit die Gestalt Vash af 

WO = 1 ist. M,,, sind die in (II, 13) eingeführten Operatoren. In unserem vall 

„Musterbeispiel‘ haben wir dementsprechend zu setzen: For 

Hier bedeutet x die „charakteristische‘“ Zahl von früher. Wir führen eine Zw 
„charakteristische‘‘ Funktion durch = 

Gon Pos (t). (III, 1) unc 


(Die Bezeichnungsweise ist dem Buch von Chandrasekhar’) entnommen.) 
Wegen a, = 1 haben wir 


(III, 2) Wi 


geh 
Außerdem ist Y eine gerade Funktion,d.h.: vee 


(—t) = (III, 3) 


Wir werden eine abkiirzende Schreibweise benutzen, indem wir setzen: 


Ta... for (III, 4) 
1 (t) 


Die 

Die Identität (III, 4) soll bedeuten, daß beide Seiten als Integraloperatoren for 
auf eine (weitgehend) willkürliche Funktion von t aufzufassen sind. Mit dieser kle 
Schreibweise ist der Operator 2 bestimmt durch 

DF=F@-( fer, F (ILI, 5) 

0) 

(F beliebig), oder auch, da Integration über t und innere Produktbildung ’) X 
miteinander vertauschbar sind, 

QF=F()- je le-¢lt, F (N) 
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Es besteht eine Vertauschungsrelation®). Sie hat in unserer Bezeichnungs- 
weise die Gestalt: 


(QD—DQ)F=F (0) (III, 6) 
Es erscheint zweckmäßig (namentlich im Hinblick auf die spätere numerische 
Rechnung), neben Y noch eine Funktion Z (s) einzuführen, die mit Y einfach 
zusammenhängt. 


-1 
Es ist noch nachzutragen, daß (im Interesse der Durchführbarkeit der mathe- 
matischen Theorie) an die Koeffizienten a,, in der Definitionsgleichung (III, 1) 
der charakteristischen Funktion Y eine weitere Bedingung (neben a, = 1) 
zu stellen ist. Sie läßt sich an einfachsten mit Hilfe der soeben eingeführten 
Funktion L (s) so ausdrücken: L (s) soll für reelle s? aus dem Intervall von 
—1 bis +1 stets positiv sein. Das ist aber, wie man leicht einsehen kann, 
äquivalent damit, daß die Gleichung Y = 0 keine (reellen) Wurzeln im Inter- 
vall von t hat. Dem entspricht in unserem „‚Musterbeispiel‘‘, daß man Yin der 
Form schreiben kann: 
(t) = (0) - (1 — #2). (III, 8) 
Zwischen s, und der charakteristischen Zahl x besteht dann der Zusammenhang: sr 
und für 7 (0) ist zu setzen: wahl 


un 


Wir setzen hier noch einige Formeln hin, die sich unter anderem aus dem vorher- 
gehenden leicht ergeben: 


2, (e-#4, DF) = 8 f (8) — F (0) 9) 


Q et? = — s2L(s)e: + Sr 


Die zweite ist die bekannte Formel Pi das „Bild“ (d.h. der Laplacetrans- 
formierten) von DF, also der Ableitung von F. Wir wählten, wie üblich, 
kleine entsprechende Buchstaben für die Bilder, z. B. hätten wir zu setzen: 

G (2)) = g (8) = L{@}. (III, 10) 
Aus der dritten Formel bekommt man unmittelbar Ausdrücke für 2 (1) 
und Q(z), wenn man beide Seiten nach Potenzen von z entwickelt und die 
Koeffizienten vergleicht: 


Q(1) = 


(t) 


Um 


8) Der Beweis ist einfach. 


hr 
bei 
Ju - 
zu 
ch 
di- | 
uf- 
‘ile 
le, 
et. 
it- 
ine 
1) 
n.) | 
2) 
NZ 
3) = 
en 
er 
5% 
7) e” zt 


Wir betrachten die Gleichung ler ‚ui 
E F, = H, (2) | a (III, 12) 


und fragen zunächst, welche allgemeinen Aussagen man über die eventuellen 
Lösungen F, dieser Gleichung machen kann, wenn an gr Funktion H, die 
Forderung gestellt ist, daß das innere Produkt (e-®?, H, (2)) für wenigstens 
alle nichtnegativen Werte von s existiert. 

Wir stellen drei Fundamentalsätze, die sich aus der allgemeinen Theorie 
herleiten ließen, an die Spitze: 

SatzI. Eine nicht notwendig beschränkte Lösung dieser Gleichung ist 
eindeutig, wenn für einen Wert von z (z. B. z = 0) der zugehörige Wert von 
F, vorgegeben ist. H, darf dazu im allgemeinen nicht identisch Null sein. 

Satz II. Wird von der Lösung F, Beschränktheit verlangt, so ist sie 
ohnehin eindeutig und strebt für große z einem konstanten Wert zu. Wenn 
aber H, = 0 ist, so existiert überhaupt keine Lösung der angenommenen Art. 

Satz III. Eine (für große 2) asymptotisch verschwindende Lösung gibt 
es dann (und nur dann), wenn die (positive) Lösung F der zugehörigen homo- 
genen Gleichung orthogonal zu H, ist, d.h., wenn 


QF=0; (fF, 


Die soeben formulierten Sätze genügen, um mit Hilfe der in (III, 9) an- 
gegebenen Formeln die Theorie zu entwickeln, soweit sie fiir unser spezielles 


Problem notwendig ist. 


(III, 13) 
gilt. 


Erläuterungen zu den drei Sätzen a ey 


Habe ich eine Partikularlösung im Sinne des Satzes II gefunden und be- 
zeichne ich diese mit Fir, so lautet die allgemeine Lösung (im Sinne des 
Satzes I): 

F (2) = Fy (2) + const - F (2). 


Die const bestimmt sich aus der Anfangsbedingung. F als (etwa bereits be- 
kannte) Lösung der homogenen Gleichung 

QF=0 
soll im folgenden immer so normiert sein, daß für z = 0 ist: 

F (0) = 1. (III, 14) 


Auf Grund des Satzes II können wir zeigen, daß F für große z einen linearen 
Verlauf annimmt. Zu diesem Zwecke führen wir eine Hilfsfunktion @ (z, t) 
ein, die von dem positiven Parameter t abhängt und die nach Satz II ein- 
deutige Lösung der speziellen Gleichung 


i 
sein soll. Wir kénnen dann schreiben: 


G (2, t) = Q> ei, 
Fir F ‚erhalten wir die ra (s. (III, 6) mit F (0) = 1) 
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woraus wegen F'= 0 folgt: 
L 1 DF =Q-1 e-zt — Qe-zt, 

stens 
woraus man die Behauptung direkt ablesen kann. 
vr: Indem wir jetzt Satz III hinzuziehen, können wir noch weitere Aussagen 

über den asymptotischen Verlauf von F (z) machen. Wir spezialisieren wiederum 
g ist die Funktion H, in unserer Gleichung, indem wir jetzt setzen 

von 
H, = H,= f ; 
omo- 

Die entsprechenden‘ Lösungen F, und F, legen wir fest durch die Anfangs- 
bedingungen: 
‚13) F,(0)=0; F,(0)=1. 
) an- Mit Hilfe der Vertauschungsrelation findet man dann: 

Wir wissen aber nach ((III, 11)) bereits, daß F, = 1. Somit ist aa 

DF, = 1 + const F. BG" 

| be- Die const ist 0, und wir haben eindeutig F, = z. Wir setzen jetzt weiter “ae 
des H,=H,=ofH, + of Ag. 


Die Koeffizienten (die Bezeichnung wegen des Späteren) af znd «$ sollen 
dabei so bestimmt sein, daß die Lösung der entsprechenden Gleichung im 
; be- Sinne des Satezs III existiert. Wir finden Aue 
0) (t) 
Die so bestimmten Konstanten ergeben den asymptotischen Verlauf von F, 
, 14) denn setzen wir 
F, + of Fy =of2-+ of, 
(2, t) so genügt die Differenz von F, — F, der homogenen Gleichung, und wir haben 
F, — F, = const F. 
Die const ist positiv und ließe sich ebenso wie die Koeffizienten af und af 
berechnen, wenn das genaue Verhalten von F (auch fiir kleine z) bekannt 
wäre. (Wir werden später sehen, daß eine so umfassende Kenntnis dazu gar 
nicht notwendig ist.) 
Die Funktion F (z), namentlich ihr Bild f (s), kann in engen Zusammenhang 
1) gebracht werden mit dem Bild g (s, t) der Funktion @ (z, t). Nach dem frü- 
heren war G (z, t) gemäß Satz II die eindeutige Lösung der Gleichung 
QG (2, t) = (II, 16) 
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die wir jetzt im Auge behalten werden. Der Operator 2 ist natürlich mit t 

und auch mit dt vertauschbar, und wir haben daher: 


Qt-—tQ=0 
[26 = Q [e = fe“. $ 


(t) (t) (t) 


Dies vorausgeschickt, können wir die Vertauschungsrelation in der folgenden 
Gestalt schreiben: 


(z,t) =6 (0,)-2 
(2) 
Die Absicht liegt auf der Hand, denn (D + t) hat, wenn auf die Exponential- 
funktion e”?! angewandt, die Wirkung eines Nulloperators; dadurch reduziert 
sich die linke Seite der letzten Gleichung dermaßen, daß wir das Bestehen der 
folgenden Gleichung behaupten können: 


=A F(a, (III, 17) 
(t’) 


worin A (t) eine zunächst willkürliche Funktion von ¢ ist, die aber identisch 
Null sein muß. Das erkennt man, wenn man das asymptotische Anwachsen 
von F für große z-Werte in Rechnung setzt, dem die linke Seite der letzten 
Gleichung folgen müßte, wenn A nicht identisch Null wäre. Die linke Seite 
bleibt aber auf Grund des Satzes II endlich. Von der nun stark vereinfachten 
Gleichung stellen wir das Bild her und finden 


(0, t) (J +1). (III, 18) 


Eine willkommene Eigenschaft der Bildfunktion g (s, t) ist ihre Symmetrie 
bezüglich der Vertauschung der beiden in ihr vorkommenden unabhängigen 
Variablen s und t. Es gilt also 
g (8, t) = (t, 8). 
Das folgt unmittelbar, wenn wir das Bild g, wie wir es getan haben, in Form 
eines inneren Produktes schreiben: 
g (8, t) = (e-**, G (2, t)). 
Dann hat man (die Variable z in@ lassen wir der Einfachheit wegen fort) 9%): 
(G (8), [2G ())]) — @ [2G (8)]) =g (s, t) —g 8). 

Das identische Verschwinden der linken Seite driickt aber gerade die Her- 


mitizität von Q aus. 
Die festgestellte Symmetrie von g (s, ¢) zieht nach sich: 


G (0, t) =tf (t) 
was wir in (III, 18) eintragen können. Denn RN wir das das Bild von 
DF= G (2, t), d.h. % wil 


9) Diese Schreibweise ist aus der Quantenmechanik geläufig! 
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nit ¢ (wir können das übrigens auch in (III, 18) eintragen), so liefert der Vergleich 
der beiden letzten Gleichungen in Zusammenhang mit der Symmetrie von g 
die Behauptung. 

Nachdem wir die entsprechenden Eintragungen in (III, 18) vorgenommen 
haben, bekommen wir endgültig 


(8 + t) g (s, t) =st f (s) f (t). (III, 19) 


iden Das ist der zu Anfang erwähnte Zusammenhang. Wenn man sich die Frage 
vorlegt, was mit dieser an sich einfachen Beziehung anzufangen ist, so wird 
man leicht dazu gefiihrt, unser so gewonnenes g in die Bilder der uns nun 
schon bekannten Lösungen F, und F, (der entsprechenden inhomogenen 
Gleichungen) einzutragen. F, = z ist zu speziell, dagegen erscheint F, geeig- 


tial- | net. Das Bild von F, = 1 ist offenbar gleich 1/s, und wir bekommen 
nlert 
Io fet III, 20 
der tat (III, 20) 
5 Von dieser Gleichung sollen alle unseren weiteren Überlegungen den Ausgang 
1%) nehmen. 
. Wir wollen diese einzige Gleichung in ein System von vier Gleichungen 
isch überführen, indem wir uns von der eingangs erwähnten Analogie zur Elek- 
u trostatik leiten lassen. Wir deuten das Bild von F, also f (s), als „‚Feldstärke‘“. 
‚ten | £(s) und n(f) seien „Potential“ und „Ladung“, und wir wollen schreiben 
sus: (die „Analogien‘‘ setzen wir in Klammern dahinter): 
W(t) n (t)tdt 1 fodt 
18) 1+ st r ), (III, 21) 


82 & (8) 


, 
sf (8) = (div E= 9). (III, 24) 


Man kann leicht einsehen, daß tatsächlich diese vier Gleichungen äquivalent 
9): der einen Gl. (III, 20) sind (wenn man noch die Definitionsgleichung (III, 5) für 
den Integraloperator 2 beachtet). Die zweite unserer Gleichungen ist zuerst von 
Ambarzumian und dann später von Chandrasekhar untersucht worden). 
[er- Wir wollen sie als „C-Gleichung‘‘ bezeichnen). Um einen bestimmten theore- 
tischen Fall vor Augen zu haben, nehmen wir von ihr an, sie wäre (etwa durch 
ein Iterationsverfahren) bereits gelöst, d.h., wir wollen die „Ladung“ für 
das Intervall0 <t<1 als bekannt voraussetzen. Aus der dritten und vierten 
Gleichung kann man dann f eliminieren und hat 


yon 


was besagt, daß, in dem angenommenen Fall das Potential & (s) für 1 < s << oo 
bei bekannter Ladungsverteilung auch mitbestimmt ist. Für große s ist der 
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x 
Verlauf sowohl von f als auch von £ sehr einfach; beide gehen wie 1/s. Denn 
nach der Voraussetzung ist F (0) = 1, d.h. lim sf (s) = 1, und die dritte 
Gleichung liefert ein entsprechendes Verhalten fiir £ (s). In der ,,Nahzone“ 
s< 1 kann man entwickeln: 


- (III, 26) 


wobei nach (III, 21) die einzelnen „Potentialmomente‘“‘ a* gegeben sind 

durch : 

ar 


(III, 27) 


Die ersten beiden Potentialmomente x und «# sind identisch mit den früher 
gleich bezeichneten Koeffizienten (III, 15), durch die, wie wir festgestellt 
hatten, das asymptotische Verhalten von F (z) für große z bestimmt ist. Die 
einzelnen Potentialmomente a} sind untereinander nicht vollständig unab- 
hängig. Wie wir sehen werden, bestehen noch gewisse Relationen zwischen 
ihnen. Eine andere Kategorie von Momenten sind die „Ladungsmomente“ « 
die definiert seien durch die Gleichung: 


0 


Sie haben, wie wir im Laufe der Untersuchung erkennen werden, eine ganz 
reale physikalische Bedeutung (sie stellen die Winkelmomente der aus der 
Grenzfläche heraustretenden Strahlen dar). Man kann die beiden ersten La- 
dungsmomente einzeln berechnen, wenn man das Potential & (s) in gewissen 
Punkten kennt. Nehmen wir unser ,,Musterbeispiel‘‘, wo 


(t) = (0) (1 — (III, 8) 


In den beiden charakteristischen Punkten + s, kann man leicht den Zusam- 
menhang finden. Setzen wir nämlich in der C-Gleichung (III, 22) s = 5, 
so wird der Nenner des Integranden Teiler des Zählers, und wir haben: 


E (8) = 2 (0) (a, — 8 


Das entsprechende Potential in dem zu s, spiegelbildlich gelegenen Punkt —s, 

erhalten wir, indem wir einfach s, in der letzten Gleichung durch —s, ersetzen. 

Wie wir sehen, stehen uns zwei Möglichkeiten zur Verfügung «, und a, aus- 

a zuwerten; einmal bei bekannter „Ladungsverteilung“ n die auf Grund der 

Definitionsgleichung gegebene und das zweite Mal durch Auswertung des 

„Potentials‘‘ € (s) in den beiden charakteristischen Punkten + s,. Die zweite 

Möglichkeit, die prinzipiell von größerer Bedeutung ist, besteht aber nicht 

mehr, wenn die beiden charakteristischen Punkte in einen zusammenfallen. und 

das ist gerade der Fall, wenn ,,Isotropie‘‘ vorliegt, d.h., wenn s, = 0 oder 

die charakteristische Zahl x = 5 ist. Das Wesentliche unserer Methode kann 

dadurch gerade gekennzeichnet werden,daB der bisher immer als ,,einfacherer“ 

„isotroper‘‘ angesehene Fall «= 5 sich als ein singulärer Grenzfall heraus- 

stellt; das äußert sich auch darin, daß dann kein Unterschied zwischen den 

Br, Potentialmomenten «% und Ladungsmomenten «,, mehr besteht, beide fallen 

zusammen, und die von uns durchgeführte Unterscheidung von Potential 
und Ladung wird wesenlos. 
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Für die numerischen Rechnungen ist es praktisch, für das Potential £ (s) 
eine Beziehung aufzustellen auf Grund der „ersten‘‘ Gleichung (III, 21); jene 
ist gültig, wenn s dem Betrage nach kleiner als Eins ist. Dazu multiplizieren 
wir die beiden „spiegelbildlichen‘‘ Potentiale & (+s) und & (—s) miteinander 
und können dann schreiben: 
fad 
Die rechte Seite wird man umformen zu 


1 fa. 
0 


t + t’ af; - 
Die in den eckigen Klammern stehenden Ausdrücke sind nach der ‚zweiten‘ 
Gl. (IIT, 22) identisch mit Eins. Beachtet man jetzt noch die von der charak- 
teristischen Funktion geforderten Eigenschaften, so erhält man die gesuchte 
Beziehung in der Form 


E (s) § (—s) = L(s), (III, 29) 


wo L(s) die in durch (III, 7) eingeführte Funktion ist. Die bei der Her- 
leitung von (III, 29) vollzogenen Vertauschungen der Integrationsfolgen 
sind erlaubt, weil unter der gemachten Voraussetzung |s| < 1 die Integranden 
stetig sind. Die soeben abgeleitete Beziehung (III, 29) spielt in der Theorie 
nicht die Rolle, die ihr in vielen Untersuchungen beigemessen wird (nach einer 
Substitution von & (s) durch 7; sie stellt im wesentlichen die bekannte ,,Fak- 
torisierung‘‘ dar). Sie an die Spitze der Untersuchung zu stellen, was natürlich 
prinzipiell möglich wäre, ist wenig zweckmäßig). Wir sehen sie als Hilfsmittel 
für die numerische Rechnung an. 

Entwickelt man beide Seiten der soeben erhaltenen Beziehung (III, 29) 
nach Potenzen von s, und vergleicht man die einzelnen Koeffizienten der Ent- 
wicklung auf beiden Seiten, so erhält man 


und = — (0). (III, 30) 


Man kann auf diese Weise eine Anzahl von Relationen zwischen den x% be- 
kommen, wenn man immer höhere Potenzen vergleicht. Für unsere Zwecke 
genügen die hingeschriebenen Relationen. Typisch ist, daß die ungeraden 
Potenzen herausgefallen sind, und es steht uns keine Gleichung zur Ver- 
fügung, die uns gestatten würde, &$ separat zu bestimmen; wir sind daher 
gezwungen, x} durch ein numerisches Verfahren zu berechnen, was auch am 
Ende dieses Abschnittes geschehen soll. Zunächst wollen wir eine Reihe von 
Relationen zwischen den Ladungsmomenten x, und den Potentialmomenten 
1% herleiten nach einer Methode, die — im Falle unseres Musterbeispiels — von 
der C-Gleichung (III, 22) äusgeht und mit Hilfe derer es z. B. auch möglich 
wäre, die soeben erhaltenen Relationen der «* untereinander herzuleiten. Zu 
diesem Zwecke schreiben wir die zweite Gleichung (ITI,22) in abgekürzter 
Form: 
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wobei hier (und im folgenden) die Zeichen f und f’ die Abkürzungen be- 
deuten: 


1 1 ‘ 
f'--- ..; (IIT, 31) dur 

Es ist also z. B. (IIT, 28): 


Wir multiplizieren beide Seiten der obigen Gleichung mit {2” und integrieren wil 
5 GERIORN beide § ‚Seiten über ¢ von 0 bis 1 und erhalten: Au 


Man kann die rechte Seite in zwei Teile aufspalten. Wenn das geeignet ge- 
schehen ist, kann man sie wieder zueinander addieren und sieht dann leicht 
ein, daß folgende Gleichung richtig ist: 


Der Grund, weswegen wir mit einer geraden Potenz von t: #2” multipliziert | ‘™ 
haben und dann integriert, wird offenbar, wenn man bemerkt, daß dann der ist, 
Nenner eines jeden Integranden Teiler des Zählers wird. Führt man die Divi- 
sion durch für n = 0 und n = 1, so erhält man unmittelbar die folgenden Be- 
ziehungen zwischen den ersten Ladungsmomenten: (di 


1 s 


Als letzte Gleichung haben wir noch eine triviale, aus der ersten Gl. (III, 21) : 
unmittelbar folgende, hinzugefügt (wenn man an das Gesetz denkt, nach dem en 
das Potential & (s) für große s wie 1/s verläuft). Die zweite Formel (III, 32) 
hätten wir auch bekommen, wenn wir die Potentiale in den charakteristischen 


Punkten miteinander multipliziert und die Formel (III, 8) benutzt hätten. Bi 
Die Formeln (III, 32) stellen drei Gleichungen dar zwischen den Momenten im 
Og % und a. Da die linken Seiten bekannt sind, könnte man denken, daß sa 
dadurch hier im Falle des Musterbeispiels (wofür die Formeln ja nur gelten) - 

diese drei Momente eindeutig bestimmt sind. Es stellt sich bei näherer Be- ZW 
trachtung aber heraus, daß die drei Gleichungen nicht unabhängig vonein- K 
ander sind. Der Rang der Substitutionsdeterminante hat, wie man sich iiber- ei 
zeugen kann, den Wert 2. Im isotropen Fall x = 5 sinkt er sogar auf 1 herab, ge 


— wieder ein Hinweis darauf, daß der Fall x= 5 qualitativ eine Ausnahme- 
stellung einnimmt —. 

Hätte man das oben angegebene Verfahren weitergeführt für n = 2 W 
und n = 3, so hätte man gefunden, daß sich außer den Beziehungen zwischen 
den Potentialmomenten, die wir vorhin RE hatten, nichts Neues ergibt. 
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IV. 

Die asymptotische Darstellung von F(z), der Lösung der homogenen 
Gl. (III, 13), ist für große z, abgesehen von einem bestimmten Faktor, gegeben 
durch: 

(IV, 1 


Da nach (111, 30) das erste Moment af bekannt ist, so bleibt nur noch eine 
numerische Berechnung von a3 durchzuführen. Zweckmäßigerweise nutzen 
wir die vorhin hergeleitete Zerlegungsformel aus. Es ergibt sich daher als 
Aufgabe, die Funktion Z (s) in das Produkt zweier spiegelbildlicher Poten - 
tiale aufzuspalten: 


E (s) (— s) = L(s), (IIT, 29) 


Die Methode, eine derartige Aufspaltung eindeutig durchzuführen, ist gut 
bekannt und trägt den Namen „Wiener-Hopf-Technik‘“. Wir wollen hier 
nicht auf die mathematischen Einzelheiten eingehen und wählen nur zur Illu- 
stration ein einfaches Beispiel, wo die geforderte Zerlegung sofort angebbar 
ist. Es handele sich um folgende Funktion: 


1 ) 

Di) = = 

(die unserem J (s) in qualitativer Hinsicht sehr ähnlich ist. Qualitativ ent- 

spräche sie etwa einer charakteristischen Funktion Y, die im Intervall von t 
nahezu konstant mit Werten zwischen $ und (3/e) - (4)). Wir setzen 


I. (8) = (s) & (s), 
Wo 


J 


1 st | | + st | 
0 
mit 
a(t) =1—- ror O < 1; a)=-1-5 für —1<t< 0. (1V,2) 
ildet man das Produkt von &* mit &-, so können +7 = 
im Exponenten die beiden Integrale in eines zu- a 2 
sammengefaBt werden, dessen Integration von : 
—1 bis + 1 läuft. Die Funktion a (ft) besteht aus ‘ 
zwei Geradenstiicken (s. Abb. 2, ausgezogene ‘ >t 
Kurve) und ist bei t=0 unstetig; sie macht dort -7 +71 
einen Sprung um zwei Einheiten. Ferner ist sie un- R 
gerade, und wir können substituieren 
(s) = &* (— 8). -1 
Wir bekommen, wie man sich überzeugen kann,in Abb. 2. Zur Berechnung des 
zweiten ,Ladungsmome n- 
L(s) = &* (s) (— s) tes“ (s. Text) 
Ann. Physik. 7. Folge, Bd.1 ty 28 
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die gewiinschte Zerlegung. Wir setzen 

und müssen das ‚richtige‘ a (t) so bestimmen, wie es der vorgegebenen a 
teristischen Funktion WY entspricht. Die zwangsläufige, aber etwas weit- 
schweifige Rechnung ergibt für das a(t) folgende Bestimmungsgleichung: 


t, t’) dt’ 
cot za (t) = d. 
(t\2 Pr) 


a (t, t’) 
Das Integral ist ein unbestimmtes (im Sinne Cauchys) und muß durch Bil- 


dung des entsprechenden Hauptwertes eindeutig gemacht werden. F 
Für den Fall unseres ‚„Musterbeispieles‘‘ haben wir zu setzen: 


t’\2 1 — 
t,t’) =(— 
( t ) 1 — 3 
und bekommen nach einer elementaren Rechnung mit s? = 3 ef 


& 3 +3(2 —5)# "Br 
a (t) = are cot [= (t a 


welcher Ausdruck in 


£ | a(t) dt 
0 
einzusetzen ist. x sind die charakteristischen Zahlen (I, 1) (x = 15 cos? 9). 
Für t-Werte kleiner etwa 4 bekommt man angenähert: 
4t 


- 


4 


a(t) =1-, 


a (t) nimmt (im Gegensatz zum Beispiel Z (s), wo a fiir t = 1 } ist) an den 
Intervallenden den Wert 0 bei vertikalem Einmiinden an. Im Abb. 2 ist der 
qualitative Verlauf gestrichelt gezeichnet, entsprechend dem Werte 5 für x. 
man Potential (s) wie früher, so findet man 


Tabelle 1 Das Verhältnis a*/af stellt sich 
also als der Flächeninhalt unter 


(2/7) der rechten Hälfte der Kurve für 
3 1 0,0372 a (t) heraus (s. Abb. 2: <0,75). 
. 4 1 — 0,0200 In der nebenstehenden Tabelle 1 
- 0 sind die Ergebnisse der Auswertung 
7 +4, +0,0504 naan ‘ er ampsonec en ege 
+3), +0,0855 für die verschiedenen charakte- 
9 +1 +0,1285 ristischen Zahlen x zusammen- 


gestellt. Eingetragen sind die 
Abweichungen des „isotropen‘ Wertes von 0,7104 für —=5, also des Wertes 
0,7104 von den übrigen, wie der entsprechende Kopf von Tab. 2 anzeigt. 
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Folgerungen aus dem Alternativsatz. 
ırak- Wir hatten den Operator 2 durch (II, 14) definiert, wobei wir zur Ver- 
weit- einfachung alle w; 52 Null setzten. Die Gl. (II, 16) läßt sich in die Gestalt 


# denn es gilt 2) (a + bz) =a E, — b E, bei willkürlichen Konstanten aundb. 
Nun muß nach Satz III ((III, 13)) die rechte Seite von (V, 1) orthogonal zur _ 
Lösung F(z) der zu (V, 1) gehörenden homogenen Gleichung sein, damit yp, 
(als Lösung der inhomogenen Gleichung) ihrer physikalischen Bedeutung 
Bil- gemäß für große z asymptotisch gegen Null geht. In (V, 1) sinddie E, = E,(z) | 
die durch (6,4) für n = 2 und n= 3 definierten Exponentialintegrale, die 
man auch durch entsprechende Substitution der Integrationsvariablen als 
Integrale zwischen den Grenzen 0 und 1 schreiben kann. Man hat dann 


1 1 1 1 
0 0 Ö 


u 
In der letzten Gleichung ersetzt man die ,,Feldstirke% f durch die Ladungs- i 
dichte mittels (III, 24) und findet dann auf Grund der Definitionsgleichung 
(111, 27) für die Ladungsmomente «,, daß das innere Produkt zwischen F — 
und E, gleich ist dem Doppelten des n-ten Ladungsmomentes. Man hat jetzt _ 
3 J 
52 9) (F,2y,) =0 = -%); (V, 3) 
also wird 3 
2 — Oy . (V, 4) 
den Andererseits ist nach (I) z, als Integrationskonstanten der zweiten Planck- — 
; der schen Gleichung (II, 4) auch gleich dem zweiten Moment der Winkelverteilung 
ir x am Rande z= 0. Das erste Winkelmoment ist wegen der Stromnormierung ge 
(II, 5) gleich — 1. Allgemein gilt Pa 
1 
sich 
oo (V, 5) läßt sich auch für n + 1,2 streng begründen, worauf wir aber nicht 
nter a 
. für MA A wollen. Die Integration im letzten Integral kann wegen der Rand- _ 
| bedingung (II, 11) natürlich auch zwischen den Grenzen —1 und 0 erstreckt | 
slle1 gedacht werden (da das restliche Intervall nichts zum n Integral beiträgt). 
tung 
kte- Die Hopt- Bronsteinsche Beziehung und Schluß 
Für den „isotropen“ Fall (x = 5) gilt die aid bekannte Formel fir die 
eigt. =y3, 


= 
; 
> rg 
fr 
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die wir, ohne Riicksicht auf die Stromnormierung, schreiben kénnen: 
0 


0 — 
Lem) way. (VI, 2) 


Diese als Hopf-Bronsteinsche Beziehung bezeichnete Formel bildet einen 
Kernpunkt des früher schon erwähnten Buches von Chandrasekhar’). 
Dieser Autor hat nicht weniger als fiinf verschiedene Beweise fiir ihr Bestehen 
geliefert. Ein unmittelbar evidenter Beweis existiert nicht, was nach unserer 
Meinung einen tieferen Grund hat, den wir jetzt aufdecken wollen. Nach 
unserer allgemeinen Theorie lassen sich zwei Gleichungen fiir yp (0) und y, (0) 
herstellen, die die Berechnung dieser Größen einzeln gestatten. Sie lauten, 
wie wir jetzt zeigen wolien, 


Yo (0) + 2 Ya (0) = 3 % > (A) ay (VI, 3) 
Yo (0)? +5 


Durch Auflösung dieses Gleichungssystems erhält man ein 


3 
H+ 5 B+ 2+ 
4 2 2 

Wie man sieht, geht unsere Formel für x > 5 direkt in die oben genannte 
Beziehung über. 

Zum Beweis von (A) und (B) bemerken wir, daß (A) nach (II, 9) selbst- 
verständlich ist; man hat dort ja nur z = 0 zu setzen. Formel (B) läßt sich 
mit (III, 32) und w, = (x — 5)/4 umschreiben in 


(VI, 4) 


Zum Beweis von (B’) schreiben wir die in (III, 32) abgeleiteten drei Bezie- 


hungen zwischen a, %, und a, noch einmal ausführlich hin, wobei wir s? 


durch 3 m ersetzen. Wir haben dann: 
2 
3 a 
| 1 = (1+ $2) 08 — (2000, — 
3 
= (1+ 0% (2) (III, 32’) 
3 
1 = (1+) — 5 (3) 


(Die Formel (2), die uns später nützlich werden wird, gebrauchen wir nicht 
zum Beweis.) Formel (3) schreiben wir so: 

(3 — O49) = — 1. (3) 

Aus (1) und (3) erhalten wir jeweils für 3 w, &,%, zwei Ausdrücke: BER 

(14) +5 wy az —1 und 20% — a2. = 

Setzt man diese beiden Ausdrücke einander gleich, so kommt if 
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Der Vergleich dieser Formel mit der quadrierten Gl. (3’) ergibt  — A 
7 


(3 — — (3 «2 — a) | = 0. 

Setzt man w, + 0 voraus, so folgt unmittelbar die Beziehung (B’), die wir 
aa beweisen wollten. Da kein Grund zu finden ist, warum die Koeffizienten bei 
ha W, = 0 (d. h. x = 5) unstetig sein sollten, so gilt diese Formel auch noch fiir 
wo x=5, womit ein Beweis der Hopf-Bronsteinschen Beziehung geliefert ist. 
er Zur angenäherten numerischen Berechnung von y, (0) und damit auch 
0) q (2) (s. I) gehen wir von der Tatsache aus, daß sich die „Momente“ 

co 
= 2 f 2" We (2) dz 
0 
direkt aus den Planckschen Gleichungen (II, 4) ergeben. Denn nehmen wir 
, 3) z. B. die dritte Plancksche Gleichung, die wir wegen 
1 2 
bi = 5 af? pe. 
ıcn 
so bekommt man 
2 3 a, dow 


es 7 Benutzt man jetzt die Gl. (2) (III, 32’), die man offenbar so schreiben kann: 
= [3 (1 — P, (%))]*. 


32") wo P, wie früher die zweite Kugelfunktion ist, und setzt man dies in die vor- 
letzte Formel ein, so bekommt man einen expliziten Ausdruck für m). Ganz _ 
analoge Substitutionen ergeben die höheren Momente. Wir stellen die Ergebnisse _ 


= 
W, = 
2 as 2 
4 
2 We 2 3 dr 7 1 


Die hingeschriebenen Momente (aber nur diese!) lassen sich mit unseren nu- 
merischen Hilfsmitteln (Tab. 1) auswerten, indem man für 2, den „isotropen“ 
Wert 0,7104 einsetzt (der, wie wir sahen, eine genügende Näherung ist). 

Der Grenzübergang w, > 0 (x > 5) läßt sich einfach nur beim nullten 
Moment durchführen und ergibt eine von mehreren Autoren auf umständliche 
Weise abgeleitete Formel (der Zusammenhang mit P, wurde nicht erkannt!) 

2 9 3 
Py (%) =n-73 
Für die höheren Momente ergeben sich keine einfachen Ausdrücke, mit Hilfe 
derer sie sich etwa mit alleiniger Benutzung der Tab. 1 berechnen ließen. 

Für den allgemeinen Fall (aber x +5) machen wir mit Eddington 
(s. I.) folgende Entwicklung: 

Ye (2) = A E,(2)+ BE, 
wobei wir uns auf die hingeschriebenen Glieder beschränken müssen. Die 
Koeffizienten A, B,C werden so bestimmt, daß die Momente m, (n = 0,1) 
und %, (0) richtig herauskommen, wenn y, durch den oben angegebenen drei- 
gliedrigen Ausdruck dargestellt ist (eine in der Statistik übliche Methode). 

Berücksichtigt man die folgenden Formeln für die Exponentialintegrale, 
die sich leicht aus den Definitionsgleichungen der E, ergeben: 


so erhält man als Bestimmungsgleichungen für die drei Koeffizienten A, Bund C 
A = 6 m_, — 36 m, + 30 m,, 
ee B= —36 m_, + 192 m, — 180 m,, 


30 m, — 180 m, + 180 m. 


iow ah bee 


Die einzelnen Werte für y, (0), m, (= 0, 1, 2) sind in der folgenden Tabelle 
zusammengetragen, wobei zu ihrer Berechnung die Formeln (VI, 6) zugrunde 
gelegt wurden (m_, ist eine sinngemäße Abkürzung für 2%, (0)). 


Tabelle 2 

| We (0) | | Mg | Mm, | Ms Ms 
Tr 3 1,7426 | 0,3838 0,1239 0,0674 | 0,0773 | 0,0778 
—0,25 4 1,7376 0,3912 0,1326 0,0723 | 0,0881 0,0887 

0 5 1,7321*)| 0,3993 0,1429 — 
+0,25 6 1,7260 | 0,4081 0,1550 0,0847 | 0,1194 | 0,1203 
+0,50 7 1,7194 0,4179 0,1697 0,0928 | 0,1428 | 0,1443 
+0,75 8 1,7120 | 0,4289 0,1877 0,1028 | 0,1753 | 0,1760 
+1,00 9 1,7037 | 0,4414 | 0,2107 0,1184 | 0,2226 | 0,2211 

= V8. 


Die Güte der Näherung fiir y, als dem dreigliedrigen Ausdruck kann man 
überprüfen, wenn man die aus der Formel 


422 poe Annalen der Physik. 7. Folge. Band 1. 1958 yA 4 Er 4 


| 


was 
exa 


ye 
sich 
Air" 
= 
2 
> 
Abb 
== ist ¢ 
läng 
2 y, 
zwe 
ent 
mor 
Inte 
a. (I, 2 
a The 
wur 
Zeit 
| 
ue 


dC 


alle 


1an 


} i, D. Lyons: Das Milnesche Problem bei anisotroper Streuung 423 


sich ergebenden Werte mit den richtigen Werten m, (s. Tab. 2) vergleicht, 
was geschehen ist, indem die angenäherten Werte m, in Klammern hinter den 
exakten Werten m, eingetragen worden sind. In Abb. 3 ist der Verlauf von 


nal 


- 


| 


a 6 


Abb.3. Die berechnete ‚„Randverteilung“ [2 y, = 3 (z+ 29) — nach (11, 9)]. (z) 

ist die Dichte als Funktion von z [= Abstand vom Rande in Einheiten der freien Weg- 

länge]. Die angeschriebenen Zahlen x = 3 bis x=9 sind die Werte von 15 cos? #. 
(Mittelung über die elementare Streufunktion); 2, = 0,7104 


2 we, für alle ganzzahligen x-Werte (außer für den Wert x = 5) angegeben; und 
zwar in Abhängigkeit von z zwischen 0,1 und 0,4. Für den isotropen Wert 
x=5 haben wir die entsprechende Kurve der Arbeit von Mark (s. I.) 
entnommen. 


Anmerkung bei der Korrektur: Der numerische Wert für das Ladungs- 
moment aus der Tabelle 1 (S. 418) wurde nachträglich durch Lösung der nichtlinearen 
Integralgleichung (III, 22 (,‚C-Gleichung‘‘)) mittels eines Iterationsverfahrens von neuem 
bestimmt. Als Norm der Güte dieses Lösungsverfahrens galt der Placzeksche Wert 
(I, 2) für z, (8.404). Uber eine diesbezügliche Betrachtung, die anscheinend auch die 
Theorie des ‚kleinen Parameters‘ betrifft, welche bei Turbulenzproblemen entwickelt 
wurde, (s. neue sowjetische Literatur), soll gegebenenfalls in einer mathematischen 
Zeitschrift berichtet werden. 


Zeuthen-Miersdorf, Deutsche Akademie der Wissenschaften zu Berlin. 
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Der Hodograf der Bewegung eines materiellen Punktes 
in der relativistischen Mechanik 


i, 
Inhaltsiibersicht 


Die Hodografengleichung wird aus den relativistischen Bewegungsglei- 
chungen abgeleitet und zwar wird nur die zweidimensionale Bewegung be- 
trachtet. Diese Gleichung ist vom Ricattischen Typ. Falls die Komponenten 
der Kraft konstant sind, stellt der Hodograf eine Aquidistante dar. 


1. Betrachten wir, im Minkowskyschen Raum, die Weltlinie eines 
materiellen Punktes, der sich im physischen Raum bewegt. In jedem Punkt 
dieser Linie kann man den Vierervektor der Geschwindigkeit mit den Kom- 
ponenten 


u; = dt 2, a 4 | 
definieren wo 
dr = — dx — dad — dz? 


das Element der Eigenzeit!) darstellt. pid erg 


Aus (1) geht hervor: 
(2) 


Die Komponenten der Geschwindigkeit des materiellen Punktes in Rich- 
tung der Raumachsen Ox,, Ox,, Ox, sind 
dx u, 
k = 1,2, 3. 
Wir betrachten diese Größen als die carthesischen Koordinaten eines 
Punktes im dreidimensionalen Geschwindigkeitsraum. 
a Aus (2) und (3) folgt, daß die Metrik in diesem Raum hyperbolisch ist?). 
Der elementare Abstand im Geschwindigkeitsraum ist durch 


do? = du? + du? + du — du? (4) 
definiert. 
Das Quadrat des absoluten Wertes der Viererbeschleunigung beträgt 


du, du,\2 du, 
Yu 1) Die Lichtgeschwindigkeit c sei = 1 gesetzt. 
2) A.Sommerfeld, Atombau und Spektrallinien, Bd. 1, F. Vieweg & Sohn, 
Braunschweig 1931, S. 709; E. Borel, Introduction géométrique 4 quelques théories 


physiques. Gauthier-Villars 1914, pag. 51—55; V. A. Fock, Die Theorie des Raumes, der 
Zeit und der Gravitation. Moskau 1966, S. 68. 
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Aus den letzten zwei Beziehungen folgt aaa 
|B| dr. (6) 
Diese Gleichung stellt die Beziehung zwischen dem elementaren Abstand 
im Geschwindigkeitsraum und dem elementaren Abstand im Minkowsky- 
schen Raum dar. 
2. Künftig beschränken wir unsere Betrachtungen auf den zweidimen- 
sionalen Geschwindigkeitsraum. Dann ist 
do? = du? + du — du. (7) 
Wenn wir niehthomogene Koordinaten v,, v, verwenden, erhalten wir 
(1 — v3) + 2 v, v, dv, dv, + (1 — - 
(1 vi — Us 2)? 


do? — (8) 


Man kann zeigen’), daß wenn eine Variabelnsubstitution durchgeführt 


wird. und zwar 25 a 


wo p?= 3 + v8 bedeutet, dann erhält der elementare Abstand die Rie- 
mannsche Form 
(10) 


P+ + ¢ 
(1- ify my 


Aus (9) folgt: 


(11) 


4 q 
3. Nun gehen wir auf die Ableitung der Hodografengleichung über. Wir 
beschränken uns auf den zweidimensionalen Fall, alse u, = 0. 
Die EEE lauten in diesem Falle®) 


Min Fe = Fy + Fy ug, (14) 


wo F, (i = 1,2) die Komponenten der Newtonschen Kraft und m, die Rih- 
masse des materiellen Punktes darstellen. 

Wenn wir beachten, daß 


8) F. Klein, Nichteuklidische Geometrie, J. Springer, Berlin 1928, S. 295. 

4) E. R. Neumann, Vorlesungen zur Einführung in die Relativitätstheorie. G. 
as Jena 1922, S. 113. 
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gehen (12) und (13) über in aul 


d ( v ) = F 


Durch Differentiation erhalten wir 
ar (1 — + (I — — v3) 


dv j 


Auflésung nach und 4 ; liefert die Bewegungsgleicbungen in der Geschwin- 


digkeitsebene : 
d 
dvy,_ F, 1 v2 he ¥ 
| = 8) — 
ey Die Hodografengleichung erhält man durch Division der beiden Gl. (18): 


Die Gl. (18) erhalten eine symmetrischere Form in der Ebene n=Pp +ig. 
I : & Wir gehen von den Gin. (11) aus. Durch Differentiation erhält man: 


dp (4 — + q*) — 2 pq- dq 


Einsetzen von (11) und (20) in (18) liefert: nodes 


d 
+99) 2 = [4 — + 16 — 216 pg 


dp F, (21) 
s[4— — 16 - B16 pg. 
Mg Mg 
Auflösung nach ergibt: 4 
d F F 
{bi} 
(22) 
2 


® Diese beiden ere lassen sich in einer einzigen zusammenfassen: 
d 
My, = 7 i P,) + (F, 4-1 F,), (23) 


wo n=p-+ ig gesetzt ist. 
Wenn man die Größe 
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(16) q 


einführt, kann man die Gl. (23) auch folgendermaßen ausdrücken: Be ze 
myst = — 4=F,—iF,. (25) 
Die Bewegungsgleichung in der Geschwindigkeitsebene ist also vom 

Ricattischen Typ. 

Ein Beispiel: 
Wenn wir F, und F, als konstant annehmen, ist (19) eine Jakobische a 
Gleichung. Die allgemeine Lösung dieser Gleichung ist ) 


(F, — F3)(— — + Vr? + F3) = K (F, —F, v4)? 
(26) 


oder 
— (Ff + K F3) — (F3 + K 3 4+ 20 VF, Fy (K — 1) + F$ = 0. (27) 


Hier ist K eine beliebige Konstante. Nach (27) ist der Hodograf ein Kegel 
schnitt. Seine Invarianten sind §): 


=A+C=— (Fi + ( 14K 


2 
BD 

d4=|BCE=K (Fi + 2. 

DEF 


Daraus folgt, daß für X > 0, ö> 0 ist, der Kegelschnitt also eine reelle 
Ellipse darstellt. 

Fir K < 0, ist 6 < 0; also ist der Kegelschnitt eine Hyperbel. 

Fiir K = 0 artet der Kegelschnitt in ein Paar von Paralellen aus, welche 
folgenden Gleichungen genügen: 

ER. D= F, v. +F, D) 4+ F?— F3=0 Wi 
D,=Fy + Fy, -YR-R=0. 

Diese beiden Geraden stehen senkrecht auf einer dritten Geraden D,, 

welehe zur Gleichung gehört: 
D, = Fy», — 0. (30) 

Man kann a daß die Geraden D,, D, durch die Punkte M und N 

laufen, wo 


= VFi+ Fi 

NI- F, PF, 

Diese Punkte befinden sich auf dem Fundamentalkreis mit dem Radius 1; 
mit den eingeführten Bezeichnungen läßt sich der Kegelschnitt folgendermaßen 
schreiben: 

D, D; = K Dj. (32) 
5) W. W. Stepanov, Vorlesungen über Differentialgleichungen. Moskau 1952, S. 41. 
Die Gleichung des Kegelschnittes kann A vi +2 Bu 


= 0 geschrieben werden. 
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Der Kegelschnitt hat also in den Punkten M, N die Geraden D,, D, zu 
Tangenten’”). 

Aus der Parallelität der Geraden D, und D, folgt, daß D, eine der Achsen 
des Kegelschnitts darstellt. 

Aus Abb. 1 ist der Kegelschnitt zu ersehen. Für K > 0, erhält man eine 
Ellipse, die sich zwischen den Geraden D, und D, befindet. Beiabnehmendem K 
bleibt die eine Halbachse der Ellipse 
OM konstant, während die andere 
Halbachse ständig zunimmt. Für 
K=0 erhält man die beiden Ge- 


'. raden D, und D,. Im Falle K < 0 
erhält man Hyperbeln. 

Da nach der Relativitätstheorie 
stets v < 1, kann allein den Ellipsen 
ein physikalischer Sinn zugeschrieben 
werden. Wie bekannt®), stellt eine 

. solche Ellipse den geometrischen Ort 
der (im nichteuklidischen Sinn) gleich- 
entfernten Punkte zur Geraden D, dar. 

Abb. 1. Die Ellipse, die sich im Inneren des } 
absoluten Kreises befindet, stellt den Hodo- 
grafen dar Der Hodograf eines materiellen 
Punktes, der sich in einer Ebene unter 
der Wirkung einer Kraft mit den konstanten Komponenten F,, F, bewegt, 
ist eine Äquidistante zu einer Geraden durch den Koordinatenanfang des 
Geschwindigkeitsraumes (s. oben). Diese Gerade bildet mit der Richtung 
der v,-Achse O v, einen Winkel x, der durch die Beziehung tg x = == ore ist. 
7) T. Lalescu, Tratat de Geometrie Analiticä. Bd. 2, Bukarest 1944, S. 48. 
8) F. Schilling, Pseudosphärische, hyperbolisch-sphärische und elliptisch-sphäri- 
sche Geometrie. B. G. Teubner, Leipzig 1937, S. 169. 
Timisoarä (Rumänien), Facultatea de Matematica si Fisica. 
a : Bei der Redaktion eingegangen am 11. November 1957. A 277 
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DEUTSCHE GLIMMLAMPEN-GES. PRESSLER | JOHANN AMBROSIUS BARTH / VERLAG /LEIPZIG 
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Finführung in die Theorie der Elektronenoptik 


Von Dr. phil. JOHANNES PICHT 
o. Prof. für theoretische Physik und Optik, Pädagogische Hochschule Potsdam-Sanssouci 
> (früher Techn. Hochschule Berlin-Charlottenburg) 
2., erweiterte Auflage 
VIII, 274 Seiten mit 70 Abbildungen im Text. 1967. Gr. SE 
DM 29.20, Leinen DM 31.— 
Die neue Auflage ist auf den neuesten Stand der Wissenschaft gebracht worden. wodurch 
sich eine Erweiterung des Umfanges ergab. Der Verfasser war bemüht, den Charakter 
einer Einführung in das Gebiet zu wahren, insbesondere die mathematischen und ge- 


danklichen Überlegungen soweit durchzuführen, wie dies zum ausreichenden Verständnis 
der theoretischen Zusammenhänge erforderlich schien. 


Klassische Arbeiten deutscher Physiker 


Herausgegeben von der Physikalischen Gesellschaft in der 
Deutschen Demokratischen Republik 


Heft 1: W. C. RONTGEN, Grundlegende Abhandlungen über die X-Strahlen 
44 Seiten und 1 Porträt. 1954. Gr. 8°. DM 2.70 


Die Technik: Es ist sehr zu begrüßen, daß die Originalarbeiten Röntgens über seine 
Entdeckung der X-Strahlen neu herausgegeben worden ‘sind. Das Studium der kleinen 
Schrift führt uns die Größe der Entdeckung Röntgens vor Augen und zeigt, daß er dank 
seiner Gründlichkeit bereits die wichtigsten Eigenschaften der X-Strahlen feststellen 
konnte. Das Studium der Abhandlungen kann sehr empfohlen werden. Freyer 


Heft 2: W. FRIEDRICH, P. KNIPPING, M.v. LAUE, Interferenzerscheinungen 
bei Röntgenstrahlen 


36 Seiten mit 3 Abb. im Text und 4 Tafeln. 1955. Gr.8°. DM 3.— 


Röntgen- und Laboratoriumspraxis: Nachdem im ersten Heft die grundlegenden 
Arbeiten Röntgens abgedruckt worden waren, folgt jetzt im zweiten Heft der Serie die 
originalgetreue Wiedergabe der Veröffentlichung von Laue, Friedrich und Knipping aus 
dem Jahrgang 1913 der „Annalen der Physik“ über Interferenzerscheinungen bei Röntgen- 
strahlen. Es ist außerordentlich reizvoll, den grundlegenden theoretischen Ausführungen 
zu folgen, die Laue über die vermutliche Wellennatur der Röntgenstrahlen und die Mög- 
lichkeit, sie mit Hilfe von Interferenzerscheinungen an dem regelmäßigen Raumgitter von 
Kristallen zu beweisen, in dieser Arbeit gemacht hat. Gj. 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 
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